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ABSTRACT 


In recent years due to requirements of modern technology and advances in 
computing, there is a considerable interest in a wide variety of inverse 
problems. These are basically IBVP incomplete in the normal sense and 
ill-posed from that point of view. These problems are solved with the help of 
additional observed data. 

The thesis addresses itself to two such problems. It consists of five 
chapters. In the first chapter we give a brief introduction to literature on 
the subject. The basic content of the thesis consists of two blocks of two 
chapters each. In chapters two and three the inverse problem called the 
source problem for coupled Reaction-Diffusion system is considered, whereas in 
the second block i. e. , chapters four and five, space 
variable-dependent-source-problem for the Heat Conduction problem is analyzed. 

In chapter two, the coupled Reaction-Diffusion system in a semi-infinite 
spatial domain is studied. The basic results consist an iterative algorithm 
for construction of the reaction function and showing the existence of such a 
function in Holder function space. Towards that end the problem is 
reformulated as a fixed point problem of a certain operator. Using the 
celebrated SCHAUDER’ s fixed point theorem, the existence of the original 
problem is established. 

In chapter three, under additional hypotheses, we have established the 
uniqueness of the above solution by the method of center contraction mapping 
property of the operator and there by simultaneously showing the convergence 
of the iterative scheme. 

The second block of results consist of determination of a source function 
depending on the spatial variable in addition to the unknown solution variable 
in a bounded spatial domain. 

In chapter four we have analyzed the existence of the source function and 
in chapter five the uniqueness and the convergence property of the iteration 
scheme are established. 
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ANALYSIS AND DETERMINATION OF THE SOURCE TERM IN COUPLED 
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KANPUR-208016, MAY, 1994 


In recent years due to requirements of modern technology and advances i 
computing, there is a considerable interest in a wide variety of invers 
problems. Most of these arise in situations where certain basic structure o 
physical models are assumed and observed data are required to be used t 
concretize the models. It is usually done either by determining unknow: 
parameters involved in the model or determining unknown dependence, expressei 
as functions of certain physical entities such as temperature, density or thi 
rate of reaction in a chemical process. The source of inverse problems coul< 
be as diverse as medicine, geophysics or space technology. 

There are many types of inverse problems involving PDE. These ar< 
basically IBVP with incomplete information such as "lack of information aboul 
certain parameters" or "the absence of knowledge about certain inpul 
functions". Determination of source term or the input function, is called th« 
source problem. While solving inverse problems one faces the lack oi 
information which has to be compensated by additional supply of observed oi 
observable data. As one would expect, these problems are generally ill-posec 

and difficult from the view of analysis, because of generally nonlinear nature 

! 

of the problems. 
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Several people have worked on source problems. Indeed, assuming the 
monotonicity of certain class of overposed data, DuChateau and Rundell [1] 
have determined the unknown reaction term in a n-dimensional single R-D 
equation. An analogous result in R 1 using a different approach was obtained by 
Muzylev in [2]. Although in [1] and [2] the uniqueness of the solution of the 
above source problems have been shown in the case when the forcing function is 
assumed to be analytic, the question of the existence was not answered. 
However, Pilant and Rundell in [3] have shown the existence and uniqueness to 
the above problem for a certain class of functions in a small time interval. 
In [4] Pilant and Rundell and in [5] Anger have studied the same problem by 
reformulating it as a fixed point problem. The present work possibly the first 
work for a coupled system of two equations dealing with reaction diffusion 
system derives results analogous to those by Pilant and Rundell obtained for 
the one dimensional heat conduction equation. Later on, source problem for 
heat conduction equation, with the source function depending on space variable 
is considered, where the more complex case of bounded spatial region is the 
underlying domain. 

The thesis consists of five chapters. In the first chapter we give a 
brief introduction to literature on the subject along with an outline of the 
contents of different chapters. The basic content of the thesis consists of 
two blocks of two chapters each. In chapters two and three the source problem 
for coupled Reaction-Diffusion problem is considered, whereas in the second 
block 1. e. , chapters four and five, space variable-dependent-source-problem 
for the Heat Conduction problem is analyzed. 

Now we summarize the basic results contained in the first block i. e. , 
chapters II and III. 

It is well known that Reaction-Diffusion equations arise in most 
chemical processes including the problems of flame propagation. One of the 



major problems is to assess the nature of the reaction. For example, 
problems of one space dimension, the equations generally take the form 


U t ~ u xx = f(u) 2 (v) ! v t ~ v xx = " f(u) 8 (v) ’ 

where f and g are to be determined. From physical considerations one of 
functions "g" is known analytically. Therefore, the source problem consists 
the determination of f from experimental data. Quite often, in practice 1 
structure of f is assumed involving certain parameters. Then using some adl 
procedures like assignment of values to these parameters and comparing wi 
experimental data, one is led to a set of parametric values for which i 
experimental results are very close. However, in the present work we ha 
given an analytical and constructive way to determine such an f withe 
assuming explicit structure of f in terms of some parameters. 

In chapter two, we have considered the following mathematical mpdel in 
semi-infinite domain (0,«) and t > 0. 


u - u = f(u)g(v) 

t XX 

(A) u(x,0) = u q (x) 

u (0, t) = g C t ) 

X 1 


V - V = - f(u)g(v) 

t XX 

v(x,0) = v q (x) 
v^(0, t) = g 2 (t) 


with superposed known data ; u(0,t) = 0(t). 

Then our quest is to determine f assuming g to be a known function of 
such that the overposed condition is satisfied. The basic results of thi 
chapter are as follows : 


BASIC RESULTS : 

(1) Presenting an iterative algorithm for construction of f. 

(2) Showing the existence of such an "f" in a certain function class. 



iv 

The basic algorithm consists of a Picard type function evaluation scheme. 
If f <n> is a given function then determine u (n) , v (n) from the direct problem 
(A) which will not in general satisfy the overposed Dirichlet data at the 
boundary. However, at the boundary x = 0, we need to satisfy the following 
relation 

9' - u = f(0) g(v), 

XX 

which is not valid for arbitrary choice of f. Thus the following updating 
technique is used i.e., knowing f (n) , u (n) , v (n) , one obtains a new function 
f (n+1) satisfying 

9 ' - u (n) = f (n+1> (9) g(v (n) ). 
xx 

This inspires the introduction of the following operator 

0' (t) - u (O.tjf.g) 

T [f 3 Ct ) - — . 

g(v(0, t;f ,g) ) 

Then the iteration scheme is given by 

f (n+1) (0(t)) = T [f Cn) ] (t) . 

U 

It is shown that the existence of the solution of the inverse problem 
is equivalent to the existence of a fixed point of T in C°\ 0 < « < 1. The 
self mapping property of T^ is established in a bounded ball in C** and using 
the compact imbedding in the scales of C a the existence is established. 

Since T 0 involves u^ , it is difficult to obtain different properties 
of the operator, therefore, it is expressed in an integral form with respect 
to u which measures the nonlinearity of the R. H. S. of (A). By using the 
Green’s function, the solution u is expressed in terms of integral expression 
which is much easier to handle. The properties of Green’s function are used to 
give the required estimates in C* for 0 < « < l. The actual bounds are 
calculated explicitly and using the property of the Identity map to be a 



compact imbedding from C* — » C f for p < «, the existence of a fixed point c 
the operator is shown. 

In chapter three, we have studied the uniqueness of the above solutio 
and have also shown the convergence of the iteration map. Extra smoothnes 
property on both f and g ensures certain estimates which lead to the require 
conclusions. The whole process involves a long list of estimates. 

The second block of results consist of determination of a source functioi 
depending on the spatial variable in addition to the unknown solutioi 
variable. The domain considered here is a bounded domain. Pilant and Rundel] 
have considered the heat equation of the following type in a semi-infinitt 
region. 

u - u = F (u), 

t XX 1 

where the source term is a function of u only. Clearly, more general form of 
involving spatial variable is of the type F^u.x). In the thesis we have 
considered F^u.x) where variation of F with respect to x is allowed up to 
linear terms in "x". Then the equation can be written as 

u - u = F (u) + x F (u), 

t xx 1 2 

which can be rewritten as 

u - u = (1 - x) f (u) + x f (u) + y(x,t) 

t xx 1 2 

in a domain (0,1) x R + . Besides, we assume the following Initial Boundary 
conditions. 

u(0, t) = u(x) 

u^O.t) = g^ ( t ) , u x (l,t) = g 2 (t). 

Since there are two unknown functions to be determined, we require two 
overposed Dirichlet data which we give at the two boundary points i. e. , j 


u(0,t) = 0(t) and u(l,t) = q(t). 



vi 

In chapter four we have analyzed the existence of the solution (u,^,^) 
and in chapter five the uniqueness and the convergence property of the 
iteration scheme. Here the Green s function is an infinite series and the 
underlying function space is C with 0 < « < 1/2. 
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CHAPTER I 


INTRODUCTION 


1.1 MOTIVATION : 

The modern inverse problems will date back to N. H. ABEL, who 
pioneering work in 1923 led to an integral equation of the first kind (KELLE 
1976 [76], ANGER, 1990 [7]). His original problem was to determine the sha 
of a hill from the time a particle taking to return when it slides up on 
frictionless hill with certain initial energy. Questions regarding uniquene: 
of the solution to this problem were answered by KELLER (1976) [76] and ANG1 
(1987) [7]. In course of time more several inverse problems were posed in th: 
connection as well as in connection with investigation of the interior of tl 
obstacle. Some of them turn out to be very important from the point of view < 
modern technology. A large number of inverse problems are motivated by tl 
desire to know the past states of physical system, or determination of unknov 
parameters involved in a physical model or to determine certain parameters 3 
order to influence a system to behave in a desired manner in future. j 

1.2 INVERSE PROBLEMS : j 

Physical problems in different fields of Science, Technology and Medicir 
often give rise to mathematical models of inverse problems. Severe 
practically important models involve differential equations with insufficient 



informations. These deal with determination of unknown coefficients, 
boundary or initial conditions or source terms from observed or experimental 
overdetermined data. The Cauchy problems for a parabolic one in backward time 
direction and Cauchy problems for hyperbolic equation with data on a time like 
manifold (JOHN, 1955 [73]) are some classical examples. 

1.3 ASSOCIATED DIFFICULTIES : 

The first question that arises in the inverse problem is the physical 
wholesomeness of the proposed model. If it is not so, then it is one of the 
many finite models that may satisfy the observed data. The reasons for this 
can be traced to (1) inherent nonuniqueness; for example potential fields are 
inherently nonunique (2) uncertainty in data; caused by experimental error and 
insufficient model specification or the mixture of these two. For example the 
information provided by the finite number of sensors inside a heat conducting 
body are at discrete locations and measured at discrete time interval results 
in incomplete data of spatial temperature variation. 

Another difficulty associated with the heat conduction problems in 
particular is, the temperature response of an internal point is quite 
different from that of a point at the surface. These are some of the demanding 
aspects encountered in inverse problems. 

1.4 SCOPE OF INVERSE PROBLEMS : 

It is easily observed how inverse problems arise in the following areas 
of scientific activity. 

(1) EARTH SCIENCE : 

The principal task of geophysics is to study the internal structure of 
the earth based on surface or sub surface observations. The shortage of the 
information on gravitational field inside the earth results in uncertainty and 



nonuniqueness of solutions. But the first uniqueness results (for special mi 
distribution) were obtained by G. HERGLOTZ (1914) and P. S. NOVIKOV (193$ 
The spherical symmetrical model of the earth played a major role in know! 
the internal structure such as crust, coating and nucleus. The coupling 
different fields in geophysics i.e., coupling of gravitational and seisn 
fields (A. S. ALEKSEEV, B. A. BUBNOV (1981,1984)) as well as coupling 
gravitational and magnetic fields (D. ZIDAROV (1968-1984), M. S. ZHDAN 
(1981, 1988)) are of particular interest. In addition, satellites introdu 
the possibility of measuring various kinds of radiation. Associated wi 
gravitational and magnetic surveying is one more inverse problem dealing wi 
the theory of potential. It is used in geophysics for locating the deposi: 
of mineral resources via gravity anomalies. The first theorem of uniquene: 
for an inverse problem of the Newton’s potential was formulated by P. i 
NOVIKOV in 1938. 

(2) SCATTERING THEORY : 

The inverse problem of scattering theory originates in the theoretica 
physics. The problem here is of determination of the parameters of sources 
diffusor or propagation media by means of the radiation received by a detecto 
(ANGER, 1990 [7]). The first uniqueness result in inverse scattering theor 
related to Sturm-Liouville operator was formulated and investigated by V. A 
AMBARTSUMIAN in 1929. Since 1964, various form of inverse problems have bee 
studied by G. BERG, N. LEVINSON, V. A. MARCHENKO, I. M. GELFAND and B. M 
LEVITAN etc.. FADDEYEV (1959) [45] determined a medium for which a! 
electromagnetic wave is reflected knowing the reflection coefficient. Similar 
questions are also asked in quantum mechanics. The paper translated by SECKLE1 
(FADDEYEV, 1959) [45] gives a nice review on the inverse scattering theory* 

I 

The early book of AGRANOVICH and MARCHENKO (1963) [1] discusses the variants 



if inverse problem arising in connection with the quantum theory of scattering 
,hich is apparently the most interesting from the stand point of application. 

3 ) MEDICINE : (COMPUTER-AIDED TOMOGRAPHY) 

Computer — aided tomography in medicine consists of the reconstruction of a 
'unction from their projection along lines (NATTERER, 1986) [95). "Tomography" 
lS derived from the greek word meaning slice. In the early 1970s CT was 
introduced in diagnostic radiology. Complete bibliography with excellent 
survey on the applications of CT is cited in (DEANS, 1983) [37]. 

1.5 ILL-POSEDNESS : 

Mathematical physics frequently give rise to ABEL’s integral equation 
rtiich is of first kind. The small structural' difference between the first and 
second kind of Fredholm integral equations changes the whole theory 
Irastically. The Integral equations of second kind are wellposed whereas the 
squations of first kind are inherently unstable. The smoother the kernel, the 
greater the risk to handle such equations. 

In the beginning of the twentieth century J. HADAMARD (1923) defined 
'well-posedness" of a problem in a classical sense. A problem is called 
/ellposed if the operator equation corresponding to it has a unique 
:ontinuously dependent solution. If there is any deficiency as regards any 
>f the above three conditions, the problem is called illposed or incorrectly 
x>sed in the sense of HADAMARD. 

A closer view of HADAMARD’ s definition ensures that the wellposedness of 
-he operator equation is intimately connected not only with the operator but 
also with the spaces and the topology they carry i. e. , it is the property of 
che triplet, space, operator and topology. 

For the problems which are not correct in the classical sense TIKHONOV 
suggested a new notion of correctness (TIKHONOV, 1963) [118]. His 



regularization method is based on the celebrated Riemann-Lebesgue lemma, 
role of the regularizing parameter is proposed to damp out the oscillations 
the solution. 

Almost all works on Inverse Problems discuss about the illposed nature 
the inverse problems. ROMANOV (1987) [106] and DONTCHEV (1993) [39] comp; 
Tikhonov Regularization with the original classical notion of HADAMARI 
wellposedness. They give a nice distinction between the two. Further tl 
illustrate that by modifying the definition of stability and assuming apric 
the existence of the solution, one can overcome the rigidity of the classic 
wellposedness. 

The mathematical tool developed by TIKHONOV’s conventional correctne 
has by now a wide variety of applications. DIMRI (1992) [38] discuss 

different aspects of this method which are proved to be applicable f 
geophysical problems. In another classification the regularization procedure 
have been classified into three categories. (1) Tikhonov regularization (: 
Constrained least squares regularization and (3) Wiener Filter Regularizatioi 
From computational point of view regularization methods are grouped togethi 
as (1) Sequential regularization method (BECK 1985) [18] (2) Trial functi< 

regularization method (3) Zeroth order regularization method (4) Generalize 
Sequential function specification regularization method . For more informatic 
see DEMRI (1992) [38] and the literature cited there. 

1.6 METHODS : 

I 

Every inverse problem requires a particular approach since no sing! 

| 

approach is suitable for all types of problems. The methods are general! 
designed to solve a particular category of problem, each category beip 
defined by the relationship between the perturbation of the model and it 
effects on the observations. In several engineering context, it is sometimej 



necessary to estimate the surface temperature in a body from a measured 
history at a fixed location inside the body. The corresponding inverse thermal 
problem may be called internal in LAVRENT’EV’S terminology (KLIVANOV, 1985) 
[80]. 

Various solution methods have been applied to the inverse problems 
including (1) Integral equation solution (2) Series solution (3) Transform 
solutions and (4) Function minimization technique. 

One of the earliest papers on the IHCP (Inverse Heat Conduction Problem) 
was published by STOLZ in 1960 [113], where he calculated the heat transfer 
rates during quenching of simple finite shapes. He obtained a linear solution 
by numerical inversion of the integral solution of the direct problem. His 
result was found to be unstable for small time steps. 

FICHERA (1961) [46] succeeded in extending the single layer theory to 

Dirichlet problems for higher order strongly elliptic differential equations 
m two independent variables. The corresponding integral equations derived 
from the boundary conditions have logarithmic Kernel, rather than Cauchy 
cernels which were used in the classical approach. This replacement has the 
idvantage of allowing simpler numerical computations and obtaining optimal 
-ate of convergence (HSIAO-WENDLAND, 1977) [66]. 

BURGGRAFF (1964) [21] devised a series solution to linear inverse problem 
rhich is exact only for continuous output data. An important drawback of this 
•orm is that they contain derivatives of arbitrary orders with respect to time 
>f the experimentally determined temperature and its gradient at points inside 
*e body. IMs immediately limits the applicability of series constructed in 
*is way because of the unreliability of the values of the derivatives. In 
•TSIREL MAN, 1985) [121] the author has demonstrated the construction of an 
alytical solution using Laplace transform which is free from this difficulty 
correct within the class of analytical functions. IMBER (1974) [67] has 



developed a transform solution for two dimensional bodies of arbitrary sha 
Details can be found in (WOO, 1981) [127]. 

BECK (1968) [16] utilized a least squares technique to generate solutii 
for a much smaller time step using the integral approach similar to STOLZ. 
(ARLEDGE et al. 1977) [8] the authors used an integral solution approach whl 
is valid for constant thermal properties. 

The finite difference and finite element methods have been t 

predominant numerical techniques for solution of the direct problem of he 
conduction and they have been applied to non-linear inverse formulation by 0 
(1975) [96], BECK (1970) [17], BASS (1980) [15] and YOSHIMURA (1985) [128]. 

In (TALENT I -VESSELLA, 1982) [117] the authors have presented t. 

stability estimates of the heat equation. The problem satisfying some grow 
condition restores stability of the solutions. MIKHAILOV (1983, 1985) [89 

[90] in an external formulation investigates the convergence of iterati< 
methods. These methods are based on the search for boundary functions 1 
starting from the requirement of minimization of a certain functional 
Iteration approximations to the desired function are constructed by tl 

conjugate gradient method. 

KEROV (1983) [77] solved the heat conduction problem in a cylindrica 
coordinate system where as in 1985, ALNAJEM and OZISIK [5] solved it for thre 
dimensional problem using the least square technique to compute the unknow 
parameters associated with solution. 

BEN-HA I M and ELIAS (1987) [20] measured the surface temperature and hea 
flux by optimal design using convexity analysis which is developed by BEN-HAIj 
in 1985 [19]. 

ALIFANOV solved the IHCP by iterative methods in 1974 [2]. But he founi 
that there is no uniform convergence in the calculation of the gradient an< 
the accuracy of the solution of IHCP depends to a considerable extent on th< 



selection of the initial approximation. But MIKHAILOV (1983) [89] modified 
the procedure slightly and got a better result. The Iteration scheme based on 
searching for the surface temperature is not only uniformly convergent but 
also assures the greatest accuracy in restoring the heat flux density. A 
number of methods for solving inverse heat conduction problems are analyzed in 
(ALIFANOV, 1983) [3] from the point of view of their practical use. 

Based on the use of gradient of the error functional and the functional 
deficiency gradient for the iterative solution of inverse parabolic problem 
various formulations of nonlinear inverse problems of generalized heat 
conduction were discussed in (ALIFANOV, 1987) [4]. 

Parametric optimization methods are analyzed in (VIGAK, 1986) [123], 
(MATSEVITY, 1988) [88]. In (VIGAK , 1986) [123] the authors have elaborated 
the result of VIGAK (1983) [122] and proposed a method to solve a control 
problem for a one-dimensional temperature system. 

MURIO (1988) [93], introduces a new automatic algorithm to uniquely 
letermine the radius of mollification as a function of the amount of noise in 
:he data. This parameter selection criterion is a very important practical 
letail when attempting to solve a real problem and leads naturally to a simple 
ind powerful computational technique. 

BANKS AND WADE (1989) [14], [124] have introduced weak Tau method for 
arameter determination problems. 

INGHAM (1992) [68] used BEM to solve the inverse nonlinear problem. 

For internal measurement please refer to (CARASSO, 1982) [33], (KLIBANOV, 
985) [80], (MUZYLEV, 1985) [94] , (WEBER, 1981) [125], (IMBER, 1974) [67], 
HORE, 1977) [63], (KEROV, 1983) [77]. 



1.7 COEFFICIENT PROBLEMS : 


(1) DISTRIBUTED PARAMETER SYSTEM : 

The measurement problem in distributed system is posed as 
observability question. Sensor location and the information content of 
resulting signals relative to a partial differential equation model are 
primary questions. New definition of observability has been introduced for ' 
purpose by the authors GOODSON and KLEIN (1970) [54]. Indeed the abc 

definition states observability as a uniqueness question and does not consic 
the subsequent estimation problem. However, the definition is independent of 
particular analytical solution form. Therefore, the full range of classic 
and modern mathematical techniques may be employed to answer observabili 
questions. 

(2) PARABOLIC PROBLEMS : 

The coefficient problem for heat operator was formulated by JONES Jr. 
1962 [74], He determined the conductivity of a medium if it was known aprioi 
to be a function of time only. It was shown that the problem is equivalent I 
a certain nonlinear integral equation for the coefficient and both posses 
unique solution simultaneously. The constructive existence proof suggested t 
KELLER is utilized to find a numerical approximation to the solution (DOUGLAS 
1962) [40]. 

In (JONES Jr. , 1963) [75] the author solves five problems using th 

integral method equation developed in 1962. The unknown coefficient is assume 

to be positive and continuous to establish the results. For special casi 

1 

explicit formula has given. 

CANNON (1963) [22] considered linear transformation with respect to dat< 
and the integral transformation for the coefficient to yield a simple analysis 

I 

of the existence and uniqueness of JONES problem with some changes . Some 



continuity conditions on the data ensures the solution of the problem. (The 
stability of the problem is also studied and it was found that [a - a^| = 
0(N -4 ) for 0 s t s T. ) 

CANNON (1967) [24] solves a nonlinear coefficient problem. The 
overdetermined data is considered to posses a continuous derivative with 
respect to arc length in order to solve the problem. In (CANNON, 1973) [28] 
the authors determined the positive function (a(u), b(u), u(x,t)> from the 

equation a(u)u fc = (b(u)u ) . They have also shown that for a constant k > 0 
such that b(u) = k a(u) there exists a unique solution provided that the 
additional data is specified where k is unknown. Both for semi infinite and 
finite intervals, using variational formulation of nonlinear inverse problem 
for the determination of coefficient has been solved in (CANNON-DUCHATEAU, 
1980) [29]. The authors have used maximum principle to solve the problem. 
Earlier a related problem was studied by ISKENDROV (1975) [70] as well. G. 

CHAVENT and P. LEMONNIER (1974), used an optimal control theory using interior 
measured data. But from the point of view of applications the problem 
involving interior measurements cause a lot of difficulties. DUCHATEAU (1981) 
[41] showed that the mapping from the coefficient to overposed data is 
strictly monotone (isotonic) and as a consequence the solution of the inverse 
problem, was shown to be unique. Under appropriate conditions, the spatially 
varying coefficient problem can be treated as a first order hyperbolic 
equation in the unknown coefficient (RICHTER, 1981) [105]. If the forcing 

function is positive and HOlder continuous in the domain and u = 0 on the 
boundary then they obtain that the solution for the diffusion coefficient does 
not need Cauchy data to obtain an apriori bound on its stability in terms of 
relevant properties of the domain, f and the unknown diffusion coefficient. 

In (RUNDELL, 1987) [108] the author points out the demerits of the 

earlier methods chosen by PIERCE (1979) [100] and SUZUKI (1983) [115] using 
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the technique of inverse Sturm-Llouville operators developed by GELFAND an< 
LEVITAN . Although they have shown the existence and uniqueness of the solutioi 
pair for spatially dependent coefficient problem from overdetermined date 
still the problem is not completely resolved. Since Gelfand-Levitan approacl 
is limited to Ordinary Differential Equations therefore attempts to recover 
a(x) (from u t - A u + a(x) u = 0 ) for x € R n with n > 1 has not proved 
successful. RUNDELL’s integral operator method is limited to one space 
variable. He also discusses the nature of the prescribed initial-boundary 
conditions to get ill or well posed ness of the parabolic equation. An example 
has been furnished to show that one can not in general obtain continuous 
dependence on the boundary data unless the class of coefficients considered is 
suitably restricted. 

The computational algorithm for inhomogeneous quasi linear heat 
conduction equation is studied in (ARTYUKH I N-NENAROKOMOV, 1987) [9]. 

A bibliography on inverse problems for parabolic equations can be found 
in CANNON’ s book (1984) [26]. 

18 SOURCE PROBLEM : 

In (SOLOVEV, 1989) [110] the author determines the source term for a 
parabolic equation where the source F(x,t) = f(x) h(x,t) + g(x,t) is given in 
this fashion. The only unknown quantity here is f(x). The overdetermined data; 
is given on the upper base of the cylinder. A good survey is presented 
connected with this problem. It is found that the problem is Fredholm solvable 
and with some sufficient conditions the solution is unique. The solution is 
determined in the function space [u e H 2+<X,1+ ° C/2 (Q T ) and f e H (fl)]. 

In (DUCHATEAU-RUNDELL, 1985) [42], the authors seek to determine an 

unknown source term in a reaction-diffusion equation. The analysis is based on ; 

! 

the observation that the over specified data depends monotonically on the 



unknown source term in the equation. This monotonicity is used to establish 
the unicity result for the inverse problem. Methods which apply to a general 
class of problems are very few. Monotonicity methods, however, apply to a 
variety of inverse problems involving partial differential equations of 
parabolic or elliptic type. The results presented there do not rely on 
special representations and hence valid for n-dimensional regions and for 
general parabolic operators, although for the sake of simplicity it has been 
restricted to n-dimensional heat operator. If the unknown function is 
restricted to analytic class then uniqueness for all choices of initial data 
follows. 

CANNON (1968) [25] proposed to determine the spatial source in an 
n-dimensional setting considering the eigen values and normalized eigen 
functions of difference of two solutions of the source problem. He establishes 
the uniqueness, and the procedure itself gives the essence of a numerical 
procedure. In order to specify the class of function where the continuous 
dependence holds, it is necessary that the source should possess a first 
derivative. 

ISAKOV (1991) [69] has used the monotonicity condition to study the 
inverse problems when the right side or a coefficient of a parabolic equation 
is unknown. He proves sharp estimates of solutions to the inverse problem in 
Holder spaces and reduces this problem to a fredholm equation. One important 
counter example emphasizes the importance of the monotonicity condition for 
the uniqueness theorem. 

In (ANGER, 1990) [7] and (PILANT-RUNDELL, (1987, 1988)) [103], [104] the 
authors have solved the inverse source problem for a heat conduction equation 
by reformulating it as a fixed point problem. This motivated us to look at 
other problems from the point of view. 



1.9 PREVIEW OF THE THESIS : 


In chapter two a system of two reaction-diffusion equations ar 
considered. As in standard forms of reaction diffusion equations, the sourc 
term is assumed to be the product of two functions. The function involving th 
reaction is to be determined from the available data measured at the boundary 
The problem is reformulated as a fixed point problem. The basic results are 
presenting an iterative algorithm for construction of f and showing th 
existence of such an f under certain conditions. In chapter three, thi 
uniqueness of the inverse source problem and the convergence of the iterativ* 
procedure are established. 

In chapter four we have considered a nonlinear source function problem 
for a heat conduction equation in a spatially finite region via the ovei 
determined data at the boundaries. The problem is again reformulated as a 
fixed point problem on a function space. Existence of the source term is 
established in chapter four and in chapter five, uniqueness as well as the 
convergence of iterative method has been proved. 




CHAPTER II 


EXISTENCE OF A SOLUTION OF THE SOURCE PROBLEM 
IN A REACTION-DIFFUSION SYSTEM 


2.1 INTRODUCTION : 

The Reaction-Diffusion equations are the simultaneous equations arising 
in almost all chemical reactions. The basic equations governing these 
reactions are derived from the conservation of mass and heat. Chemical 
reactions are classified into homogeneous and heterogeneous processes, 
depending on whether they take place in a single phase or at an interface. In 
both cases, the rate of reaction is a function of the temperature and of the 
concentrations of the reactants (FIFE, FRANK-KANETSKII, STACKGOLD) [47], [48], 
[ 112 ]. 

Often, in practice, the gradients corresponding to both concentration 
and temperature are present in the same system and heat transfer and diffusion 
occur simultaneously. The process then becomes more complex and entirely new 
phenomena, known as thermal diffusion and diffusion thermo effect, arise. This 
is essentially due to the fact that the heat flow depends not only on the 
temperature gradient but also on the concentration gradient and the diffusion 
flow depends not only on the concentration but also on the temperature 


gradient. 
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The rest of the chapter consists of six sections. The inverse problem is 
formulated in section 2. Some basic notations and assumptions are stated in 
section 3. In section 4 an iteration scheme is developed and in section 5 it 
is shown that the original problem is equivalent to a fixed point problem. 
Existence of the fixed point is shown in section 6. The chapter ends with a 
brief conclusion in section 7. 


2.2 FORMULATION OF INVERSE PROBLEM : 


In this chapter we consider the inverse problem dealing with the 
determination of source term in a system of two reaction diffusion equations 
where the source term is assumed to be the product of two functions fCu) g(v) 
where g is assumed to be known. Our basic problem consists of an initial 
boundary value problem for a system of two coupled equations given in a 
semi-infinite domain. 


and 


d u - 3 u = f(u) g(v) , x > 0, t > 0 ; (2.2.1a) 

t XX . 

u(x,0) = u q (x) , 0 — x < oo ; (2.2.1b) 

u x (0,t) = gi (t) , t > 0 ; (2.2.1c) 

3 t v - S^v = - f(u) g(v) , x > 0 , t > 0 ; (2.2.2a) 

v(x,0) = v q (x) , 0 i x < o) ; (2.2.2b) 

v x (0,t) = g 2 (t) , t > 0 ; (2.2.2c) 


where f is an unknown function of u and g is a known function of v. Since the 
data available is insufficient so as to solve the inverse problem (that is to 
determine the structure of f) one needs some additional informations such as 
the values of the solution at the boundary. Assume that the measured value at 
x = 0 is given to be 
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u(0,t) = © ( t ) , t > 0 ; (2.2.3) 

Now the overdetermined data is imposed on equation (2.2.1a) at x = 0. Hence it 
will become 

au(o,t) - a u(o, t) = f (u(o, t ) ) g (v(o, t) ) 

t XX 

that is, 6'(t) - U xx (0,t) = f(6(t)) g(v(0, t) ) (2.2.4) 

2.3 NOTATIONS AND ASSUMPTIONS : 


We shall use the following notations throughout the rest of the work. 

n = | (x, t ) | 0 < x < tn, 0<t<l| = R + x (0,T) (2.3.1) 

The following notations for norms and semi-norms are defined for 
functions defined on 0 or on appropriate domains : 

|u(x,t)| : (x.t) e £) 1 (2.3.2) 



I u( . , t ) | 


a 



|u(£,t) - u(r), t) 


le - -n\ 


a 


: € , V € R 


(2.3.3) 


||f|| w = sup | |f(£)l : € e dom (f)| 

f |f(€)-f(u)| 

If I = sup -j , £ , t) e dom (f) • 

“ Z*T> [ l€ - 7 , 1 “ J 


(2.3.4) 


(2.3.5) 


+ Ifl 

00 oc 


(2.3.6) 


C and C denote respectively the space of the continuous functions and 


the space of Lipschitz continuous functions of order « on appropriate 
domains. 



The role of the time variable is one of the compl icating factors, so the 
usual semi-norm for solutions of parabolic equations is not considered here. 
The following assumptions are made on the data. 

Al: u and v e C 2+ “ [0,®) and g = g = 0. 

0 0 * * 

A2: The over posed data 6 is a monotone function, 0" e C !0,») 

and inf |0'(t)| 2t 6 > 0 for some 5. 

taO 

A3: f , g € B £ s | || . | a s E J- for some constant E. 

2.4 ITERATION SCHEME : 

In order to solve the inverse problem ( 2 . 2 . 1 ) and ( 2 . 2 . 2 ) numerically, an 
iteration scheme is developed. If f(u) were known then ( 2 . 2 . 1 )- ( 2 . 2 . 2 ) would 
define a well posed problem for the functions (u,v). MORGAN, HOLLIS, MARTIN, 
PIERRE [92,61,60,62] and others have studied the direct problem dealing with 
global existence and the well-posedness of the reaction diffusion systems 
under various boundary conditions. If one assumes f = f*°* satisfying 

desirable properties, we shall have the solutions u = u <0) = u(x, t ; f <0> ,g) and 
v = v <0) = v(x,t;f (0> ,g) . 

However, u <0) , so obtained, will not. In general satisfy the superposed 
B. C. (2.2.3). So the function f <0) is updated with the help of the equation 
(2.2.4) to a new function f (1) defined as follows. 

f (1) (0(t) ) g(v <0) (0, t) ) = 0'(t) - u (0,t;f (0 !g) (2.4.1) 

XX 

The basic idea behind the updating scheme is the following. If "f" 
proper choice of function, then 


is the 
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u t (°,t) - u CO, t) = f (u(0, t) ) gCvfO.t); 


will be satisfied. Since u(0,t) = e(t). 


f (0(t) ) g(v(0, t) ) = e'(t) - U CO. t) 

XX 

The above expression motivates the updating procedure. 

This procedure is continued and the general iteration scheme is given by 


f (k+1) (e(t)) g(v(o,t;f (k \ g )) = e'(t) - u xx (o,t : f (k) ,g) 
for k = 0 , 1 , 2 ,.... 


(2.4.2) 


In view of the equation (2.2.1), we have the following equivalent form 


for the iteration 


f Oc+D ( 0 ( t ) ) g(v(0,t;f (k \g)) 


= 0'(t) - u" (0,t;f lK \g) + f (k) (u (k) (0, t) ) g(v (k) (0, t) ) 


= f <k> (e(t ) ) g(v (k) (o, t) ) 


.00, (k) 


+ j f (u (0, t) ) - f lKJ (0(t)) g(v (k) (0,t)) + 0' (t) - U< k) (0,t;f (k) ,g) 


= f (e(t)) g(v (k) (o, t) ) + f [e(t) - u (k) (o, t ) ] 


(2.4.3) 


Where F is a nonlinear map. 


2.5 EQUIVALENCE RELATION : 


e'(t)-u (0, t;f ,g) 
T 0 [f](t) = ” 

g(v(0, t; f ,g) ) 


(2.5.1) 



denote a nonlinear map in the function space B £ . Then f is called a ©-fixed 

point of T fl , if f(9(t)) = T 0 [f](t). 

In view of (2.5.1) the iteration scheme can be written as 

f Ck+1) (©(t)) = T e lf (k) 3(t) (2.5.2) 

The following lemma establishes the equivalence of the original inverse 
problem and determination of a fixed point of the map T 0 . 

LEMMA 1 : 

Let f satisfies the Lipschitz condition. Then the pair (u,f> Is a 
solution of the Reaction-Diffusion equations if and only if f is a ©-fixed 
point of T 0 . 

PROOF : 

Let us assume that {u,f> is a solution of R-D equations. Our aim is to 
show that f is a fixed point of T . 

Since {u,f} is a solution, it satisfies the over posed data, 

u(0,t) = e(t) 

Therefore, 

f(©(t)) « f (u(0, t) ) 

u (0 , t ) - u (0, t;f ,g) 

t XX 

g(v(0, t;f ,g) ) 

©' (t)-u (0, t;f ,g) 

= XX 

g(v(0, t;f ,g) ) 

= T 0 [fHt) 

Hence f(0(t)) = T 0 [f](t). This implies that f is a ©-fixed point of T . 

U 
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CONVERSELY: 

Suppose f satisfies the Lipschitz condition and f is a 0-fixed point of 
T . It is sufficient to show u(0,t;f,g) = e(t). 

Indeed, from (2.2.1), for x = 0, we get 

f(u(0,t)) g(v(0, t) ) = u (0, t) - u (0, t; f ,g) 

t XX 

<=> u xx (0,t;f,g) = u t (0,t) “ g(v(0, t) ) 

Since f is a 0-fixed point of T , 

0 

6 ' (t) - u (0,t;f,g) 

f(e(t)) = t [f ) (t) = ™ 

g(v(0, t) ) 

=> f (0(t)) g(v(0,t)) = 0' (t) - u (0, t;f ,g) 

XX 

= 0'(t) - u t (0,t) + f (u(0, t) ) g(v(0, t) ) 

* e'(t) -ujo.t) = |f(e(t)) -f(u(0,t))] g(v(0,t)) (2.5.3) 

Let a(t) = 0(t) - u(0, t) (2.5.4) 

If f has Lipschitz constant C and 

|g| = sup |g (v (x, t ) | 

CO 

V 

Then (2.5.3) can be estimated as 

l«'(t)| | [f (0(t) ) - f (u(0, t) ) ] | |g(v(0, t) ) | 

s C |g| |0(t) - u(0, t) | 

00 

s Cjaft)! (2.5.5) 

where C = C Igl 
l 6, « 

We claim a(t) ■ 0 for all t. 
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Note that the compatibility condition ensures that e(t) must agree with the 
initial condition m q M at x = 0. So for t = 0, cc(0) = 6(0) - u(0,0) = 0. 

Now suppose a(t) * 0. Represent a(t) as 


Then 

1 <x ( t ) I s f I a' (t) I dT £ C f k(T)| dT 
•»o J o 

Thus by Gronwall’s lemma cc(t) =0, is a contradiction. So our claim 
implies, 

u(0,t;f ,g) = e(t) 

In other words (u.f) is a solution of the R-D equations. This completes 
the proof of the Lemma. | 


2.6 EXISTENCE OF FIXED POINT : 

The main aim of the chapter is to show the existence of a fixed point of 
the R-D source problem. Therefore, in this section it has shown that for 0 < « 

< 1, T : C** » C“ and from SCHAUDER’ s Theorem one can get the required 

result ( see FRIEDMANN [49]). 

Before proceeding to prove the main result of this section we perform 
some preliminary preparations for subsequent use. Consider the comparison 
function which satisfies the differential equations, 

8 Or(x,t) ~ 5 0 (x,t) = 0 , x > 0 , t > 0 ; 

t XX ’ 

V»(x,0) = U Q (x) , x > 0 ; (2.6.1) 

tf x ( 0 ,t) = gi (t) , t > 0 ; 

If one writes r = u - 0, then r will satisfy 
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9 t r - d^r = f(u) g(v) , x > 0, t > 0 ; 

r(x,0) =0 , x > 0 ; (2.6.2) 

r (0, t ) =0 , t > 0 ; 

X ’ 


Hence the solution (2.6.2) can be represented as (see CANNON [26]) 

r=u -0= f f K(x,y, t-r) f(u(y,x)) g(v(y,x)) dy dx (2.6.3) 

^0 ^ 0 

where K is the Neumann function for the semi-infinite domain of the 
homogeneous equations (2.6.1). The expression for K is as follows (STAKG0LD 
[ 111 ]). 


K = K(x,y,t) 


_1 

ZVn~ 


2 2 
“(x-y) / 4 t -(x+y) / 4 t 

* + 

vr vr 


(2.6.4) 


This Neumann function satisfies 


.00 00 
I K (x,y,t-x) dy = 0 = | K (x,y,t-x) dy 

j n X J r. Xx 


(2.6.5) 


Hence for any function heC, 0 < « < l, the following is true 

J J K xx (x,y,t-x) h(x,x) dy dx = J h(x,x)j^ J K^fx.y.t-x) dy j dx 


= 0 


Therefore, we have the following identity property. 


( 2 . 6 . 6 ) 


u: 


K (x,y,t- 

XX 


x) h(y,x) dy dx 


= j f K„(x.y.t-x)[l 


h(y,x) - h(x,x) dy dx 


(2.6.7) 

During the course of the proof of the main result, the above identity is 
repeatedly used and this is the basic tool for proving the theorem. 
Differentiating (2.6.3) twice with respect to x, we get 
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t, CO 

r =u - 0 » f f K (x,y,t-x) f{u(y,x)) g(v(y,x)) dy dx (2.6.8] 

xx xx *xx J 0 J 0 xx 

t 00 

* u (x,t;f ,g) = 0 (x,t) + f f K (x,y,t-x) f(u(y,x)} g(v(y,x)) dy dx 

XX xx xx 


Evaluating at x = 0, we get the following relation 


u (0,t;f,g) = 0 (0,t) + K [f(u) g(v)] 

XX XX 


(2.6.9) 


= 0 fc + K [f(u) g(v); 


( 2 . 6 . 10 ) 


where 


t 00 

K [f(u) g(v)] = f f K (0,y, t-x) f (u(y,x) ) g(v(y,x)) dy dx (2.6.11) 

J ^ XX 


In view of the relation (2.2.4) at x = 0, (2.6.10) reduces to 


T 0 [f)(t) = f (8(t) ) 


6'(t)-u (0, t;f ,g) 

XX 

g(v(0, t;f ,g) ) 

e'(t) - 0 (0,t) + 0 (0, t) - u (0, t; f ,g) 

g(v(0, t;f ,g) ) 


e' (t) - 0 t (O,t) - K [f (u) g ( v ) ] 
g(v{0. t;f ,g) ) 


( 2 . 6 . 12 ) 


f(6(t)) = T 0 [f ] (t) 


0 (t) - 0 t (O,t) J J K xx (0,y,t-T)f(u(y f r))g(v(y,T)) dy dt 


g(v(0, t; f ,g) ) 
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Hence , 


f (© ( t ) ) g(v(0, t;f ,g) ) = 


0'(t) - * (O.t) 



KJO.y.t-x) f (u(y,T) ) g(v(y,T)) dy dr 


(2.6.13) 


We now state the following Theorem. 


THEOREM 1 : 

If the assumptions (A1)-(A3) are true, then for sufficiently small T, 

T maps C** [0,T] into itself. 

8 

PROOF : 

oc ~ 

Given a function f e C , define the function f by 
f (e(t)) = T 0 [f](t). 

For a monotone function "6", f is a single valued function. 

To show T q : C** » C**, it is enough to show j| f is bounded by || f fl^ 

for 0 < ot < 1 . Equation (2.6.13) can be written in an equivalent form, by 
applying the identity property (2.6.7) at x = 0. Omitting the dependence of v 
on f, g in the notation, the following equation is obtained. 


f (e(t)) g(v(o, t) ) = e'(t) - ^ t (o,t) 


-n 


K (0,y, t-r) 

XX 


f (u(y,r) ) g(v(y,r) ) - f(u(0,x)) g(v(0,x)) dy dr 


(2.6.14) 


By adding and subtracting some suitable terms, the bracketed expression 
appearing within the integral sign becomes 
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f(u(y,T)) g(v(y,x)) - f(u(0,x)) g(v(0,x)) 

= (f(u(y,T)) - f (u(0,x) 5 ] g(v(y,T ) } 

+ f (u(0,x) ) [g(v(y,T)J - g(vCO,x))3 (2.6.15) 

Assuming f , g e C°\ (2.6.14) can be estimated as 


I f (e(t))| s - [ le'(t) - * (o,t)| 

|g(v(0, t)) I L 



|K (0,y,t-x) | |f(u(y,r)) g(v(y,x)) - f (u(0,t) ) g(v(0,x))| 

XX 


dy dxj 
(2.6.16) 


Consider the second term of (2.6.16), which is denoted by 1^ and 
substitute equation (2.6.15) into I . Since f and g are in C** CO.oo), then it 
becomes 


I 

l 



(0,y, t-x) | 


(0,y,t-x) | 


If (u(y,x) ) - f (u(0,x)) | lg(v(y,x)) | 
+ |f (u(0,x) ) | Jg(v(y,x)) - g(v(0,x) ) | 

[ l u( y- T ^ ~ u(0,t)|° c 


dy dx 


+ ,f, oo l g i<x I v( y> T ^ “ v(0,x)|“ 


dy dx 


Now, expanding u and v at the point y = 0 and assuming u (0,t) = g (t) = 0 and 

x X 

vW.t) = g 2 (t) = 0, we get 

5 J j ,K xx (0 , y' t_T)l [ |f| « •sl eo (ly| 20< /2 0< ) (|u | )“ 

J 0 J 0 L<xco XX CO 


+ \f\ a Igl ( ly l 2oe /2 oe ) (|v | )**] dy dx 

w oc XX CO 
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Since (see Appendix A , (A6)) 


t -00 

I I |K (0,y, t-r) | |y| a dy dx £ -i- C («) t°‘ /2 

K J„ xx 2<x 3 


0 •'0 


one gets, 


I, s 2"'“* 2> c 3 (2 ° ,) lgl„ m„Uvj)“] 

Hence, the sup-norm of f can be given by 


|f | £ (inflgl) 


' 1 [ le ' - <\i. 


+ 2 (<x+2) c^(2<x) (tVoc) x 


I f l Igl (iu I ) + Igl If I Civ 

oc CO xx 00 CX 00 XX 


'-1 ] 


(2.6.17) 


(2.6. 18) 


In order to compute the a-norm of f, which is defined as || f = | f | 

+ I f I , we now set out to compute the «-semi norm of f. Without any loss of 
generality it is assumed that t^ > t^. We calculate the value of f (0 ( t^ ) ) - 


f(0(t )) first. 
2 


f (eft ) ) - f ( 0 (t ) ) = 

1 2 


t 00 f «■ 

e'(t ) - 0 (0,t ) - f 1 f K (0,y, t -t) f (u(y,x) ) g(v(y,x)) dy dr 

1 1 1 J 0 J o xx 1 L J 


gfvfO.t^f.g)) 


t oo r* 

o ' (t ) - 0 (0,t ) - f 2 f K (0,y,t -x) f(u(y,x)) g(v(y,T)) dy dx 

2 t 2 J J xx 2 

O 0 *- 


g(v(0,t 2 ;f,g)) 


Therefore, 
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1 


f(e(t )) - f(e(t )) 

1 2 gCvCO.tj) ) g(v(0.t 2 )) 


J g(v(0,t 2 )) |e' (t^ - 0 t (o,t i )| - gCvfo.t^) je' cy - yo.yj 


. A 00 


- g(v(0,t 2 )) J 1 | K xx (0,y,t r x) f (u(y,x) ) g(v(y,x)) dy dx 


'O "0 


*t A C 0 


g(v(0, t )) I 2 j K (0,y,t -x) f (u(y,x) } g(v£y,x}) dy dx 
1 xx 2 


0 0 


At this point, only the following terms are considered. 


A co 


g(v(0, t ) ) f 1 f K (0,y,t -x) f (u(y,x) ) g(v(y,T}) dy dT 
2 Jr > xx 1 


'o v o 


*t „co 


- g(v(0, t ) ) f 2 f K (0,y,t -x) f(u(y,x)) g(v(y,T)) dy dx 

1 J rk Jo xx 2 


f 0 *'0 


Adding and subtracting 


A oo 


g(v(0,t )) f 1 f K (0,y,t -x) f(u(y,x)) g(v(y,x)) dy dx 
1 J J ^ XX 1 


0 "0 


from the above expression, one gets 


rS r 00 

+ g(v(0,t 2 )) j J K xx (0,y,t i -x) f(u(y,x)) g(v(y,x)) dy dx 

r\ r 00 

- gCvCO.t^)) J j K^CO.y.^-x) fCu(y.x)) g(v(y,r)) dy dx 


JX > 


g(v(0,t i )) J Kyo.y.^-r) f (u(y,x) ) g(v(y,T>) dy dx 


rS f W 

- gtvfO,^)) ^(O.y.^-x) f(u(y,x)) g(v(y,x)) dy dr 


(2.6.19) 
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P i .t j) 

= g(v(0,t ) - g(v(0,t )) j j K (0,y,t -x) f(u(y,r)) g(v(y,x)) dy dx 

L “ J J o J o 


r t. oo 

gCvCO.t^)^ J 1 J ^(O.y.tj-x) f (u(y,x) ) g(v(y,x) ) dy dx 


t. 00 _ 

- f 2 f K (0,y,t -x) f(u(y,x)) g(v(y,x)) dy dx 
J o J o xx 


( 2 . 6 . 20 ) 


Similarly, the other two terms of (2.6.19) can be written as 

g(v(o,t 2 )) |e' (tj) - ^(o.t^j - g (v(o,t i )) je'(t 2 ) - ^ t (o,t 2 )| 
= je' (t^ - ^(o.y). x |g(v(o,t 2 ) - g(v(o, t j )) j- 
+ g (v(o,t )) je'(t ) - if) (o,t ) - e'(t ) + 0 (o,t jl 

111 tl 2 t 2 I 


Combining (2.6.20) and (2.6.21), (2.6.19) can be represented as 


( 2 . 6 . 21 ) 


f(0(t )) - f (6(t )) 
1 2 


g(v(0, ) ) g(v(0, 2 


T^[f 


g(v(0, t 2 )) - g(v(0 


•V’} 


- { r* r 

V J 0 J 0 


K (0,y,t -x) f(u(y,x)) g(v(y,x)) dy dx i 

XX 1 J 


|g(v(o,t 2 )) - g(v(o,t i ))l x 

gfvfo.tj)) |e' (tj - ^(o,^) - e' (t 2 ) + ^ t (o,t 2 )| 


+ g(v( 0 , t 


,» [{( 1 J K^fO.y.t^r) f (u(y,T) ) g(v(y,x)) dy dxj- 
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a t JO 

I. 


K (0,y,t -t) f(u(y,T)) g(v(y,xj) dy dx 

xx 2 


( 2 . 6 . 22 ) 


Now 


all the terras of (2.6.22) are simplified one after the other. 


Consider I first, where 


t on 

I = f 1 f K (O.y.t -t) f (u(y,T>) g(v(y,x)) dy dr 
2 J J xx 1 

J 0 J 0 


t 00 

- f 2 f K (O.y.t -t) f(u(y,T)) g(v(y,x) ) dy dx 
J o ” 2 


Since t < t , therefore, 
2 1 


■fT{[ 

J 0 J 0 v L 


K (0,y, t -x) - K (0,y, t -x) 

xx 1 xx 2 


f (u(y,x) } g£v{y,x) } \ dy dx 


} 


. JO 


J 1 J K xx (0,y,t r T) f (u(y,r) ) g(v(y,r) ) dy dx 


t ' o 
2 


I = I + I 

2 2,1 2,2 


(2.6.23) 


These two terms are separately estimated. Using the mean value theorem, 

I . can be written as 

2,1 

l z,\ = 1 2 J { [ K xx (0 ' y ’ t r T) ~ K xx (0,y,t 2 ~ T) ] f(u( y* T)) g(v(y. T )) j- dy dx 

pt JO t / r 

~ J 2 J J 1 | K xxt (0,y,1 “ T) [ f(u( y> T ^ g(v(y,rj) 

0 0 t 2 


- f (u(0,x) ) g(v(0,x) ) \ dl dy dx 
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*t *«> *t / P , 

= J 2 J o Jt* \ K xxt (0 ' y,1-T) [ g (v (y. T ^ | f(u(y,T)) - f(u(0,x)) | 

+ f (u(0,x) ) | g(v(y,x) ) - g(v(0,x))j j | dl dy dx 


Taking the estimate on both the sides 


H 2fl l “ J 2 J J 1 { lK xxt (0, y ,i-x)l J lg(v(y,x) ) | |f (u(y,x) ) - f(u(0,x))| 

y 0 1 2 

+ If (u(0,x)) | |g(v(y,x)) - g(v(0,x)) | j- dl dy dx 


s |' 2 J” j‘‘ { "w' 0 -*- 1 -*’ 1 [ |f| « I*'. flyl 2 “/2“) 


i2 "“ [ l f i. 


+ If I Igl (|v I )* ( |y| 2 ° c /2° c ) 1 1 dl dy dx 

oo oc XX 00 J j J 

Igl ( |u | )* + |f | Igl ( | V | )“ ] X 

00 XX 00 co ° oc XX 00 J 


„00 


f 1 f 2 f |K (0,y, 1-x) | |y| 2oc dy dx dl 

j X j A J A 


t •'O *'0 

2 


Making use of (A7) from Appendix A one can derive 


s 2 


[ lfl “ 


igl du I ) + If I Igl 

00 XX 00 00 oc 


( I v I )“ ] C (2«)f 1 f t2 |l-x| <x ‘ 2 dx dl 

xx® J 4 J t J 0 


Therefore, I can be weighted as 

x 


'V * 


2 K C (2<x) 


«(i 


— — if i igi dv \ j* + ifi igi du i j* it -t r 

x 00 OC XX 00 OC 00 XX 00 I 1 2 

-oc) L J 


(2.6.24) 



quatlon (2.6.23) is estimated as follows. 


Now I , • ‘I* ° ther te ™ ° f 6 
2 > 2 


- f* 1 f°° K (O.y.t -t) f(u(y,x)) g(v(y.x)) dy dT 

8.2 "J t J 0 XX 1 


= f 1 f K xx ( 0 ,y,t r r)[f(u(y.x)) g(v(y,x)) - f(u(0,x)) g(v(0,x))J dy dx 

J t J 0 xx 

2 

. p f ^(O.y.yx) [ g(v(y,s)) { f(u(y.T)) - flulo.x))} 

J t J o 

2 

+ fCu(O.T)) | g(v(y,x) ) - g(v(0,x) ) | j dy dx 


taking the modulus on both the sides, we have 


- P 1 «.» 

J* J n 


t 

2 


> y ( t^ — T ) | lg(v(y,x))l |f (u(y,x) } - f (u(0,x) ) 1 

+ |f (u(0,x) ) I lgCv(y.x) ) - g(v(0,x) ) 1 dy dx 


T „ , t .<x , . , 2 cx a, 

. * ik (° > y ,t i"' )l [ lfl “ lEl " - * y 

J t J o L 

2 

+ If! Igl (Iv | )“ ( iy| Z0£ /2° t ) ] dy dx 

00 <X XX 00 

*{ I8l. 'Iv 1 - 5 " * lfl » lgl - Uv - I » , “} 




IK (O.y.t -x)i (lyl 2 ^) dy dx j 

XX 1 J 


C (2«) 
- (oc+2) 3 


— f ifi. i«i. uuj.r * |f| » |gl « hr 1 . 1 

OC 


(2.6.25) 



Therefore, (2.6.23) becomes. 


g * igj + ig 2 i 


C(«) 


v ^ r i 

7T— T~ If 1 oo ,g! oc C |V I + IfL Igl UU | )“ It -t I' 

l 1 "<x J CO OC XX oo OC co XX co 12 


(2.6.26) 


Where, 


C(«) = 2 -<X C (2a) + 2 ~ (cx+2) (1-oc) C (2a) 

4 3 


Now represent the first term of (2.6.22) by I, 


x 3 = | g(v(0,t 2 )) - g(v(0,t i )) | 


-t . -CO 


f 1 f K (0,y,t -x) f(u(y,x)) g(v(y,x)) dy dx 
J J /-% XX 1 


0 " 0 


- 'V 


£ | |g(v(0,t 2 )) - gfvCO.tj))! J- 


JX> 


J 1 J l K xx ^°,y, tj-x) I / |f (u(y,x) ) g(v(y,x) ) - f(u(0,x)) g(v(0,x))| dy dx j- 


o 'o 


One can expand v(0,t ) considering the centre at t and the radius as 

1 2 

(t - t ). since g e C* we have from (A6) 

12 

his 2 -(oc+2) | g | i v (o , t ) - V (o, t ) r C3(2) t* 

3 oc 2 1 1 


Igl I f I Civ I ) + Igl Ifl (lu I ) 

OC 00 XX CO 00 OC XX oo 


1 


C (2«) 

* 2-*~ 2 > igi. iv r it r t r_i — t“ 
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Denoting 

i 4 = |g(v(o,t z )) - g(v(o,t x )j je'(t x ) - * t t0 'V} 

the second term of the equation (2.6.22) is estimated to give 

|I| < |g| |t -t I” ( I v | )** |0'- ip I (2.6.28) 

1 4 1 & cc 1 1 2 too t £0 

Expressing the third term of (2.6.22) as I g and taking modulus one gets 

|i 5 | i |g| w | 0' (t x ) - ^ t (0 ( t x ) - e'(t 2 ) + 0 t (O.t 2 ) ] (2.6.29) 


Now combining (2.6.26), (2.6.27), (2.6.28) and (2.6.29), we obtain the 

following estimate from (2.6.22) : 


| f(e(t x )) - f(0(t 2 )) I 


1 [|g| |8' 

|g(v(0,t x ))| |g(v(0,t 2 ))| L 


(t ) - * to.t ) - e' (t ) + * (o,t )} 

1 t 1 2 t 2 


♦ igUVt/ (Iv t i w ) a 


+ 2 is'oc 'Vv (|v tU s 


C (2tx) 


jigL lf L 

I <x 00 


Uv JJ + Igl If I (lu I ) 

XX co 00 oc xx 00 


■} 


C(cc) f . - 

* |8| » ‘VV“ { Ifl „ I*I« - m. Igl. ciuj.)“ } ] 


(2.6.30) 

So the «-semi norm of f is obtained by dividing both the sides of 
(2.6.30) by lett^ - 0(t Jl* and taking the supremum. 
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Thus, 


I f (e(t i )) - f(e(t 2 ))| 

I e ( t ) - eft) 

1 2 


s |inf |g| 2 j 


.-2 It - t | 
2 12 


left ) - eft ) r 

1 2 


Igl. 


10' (t i ) - 0 t (O, t 4 ) - 0'(t 2 ) + i/» t (0,t 2 )| 

it - 1 r 
1 2 


+ (|v tL )0C l 0 ' -* t L + igL it.i* ( |v t L )CC 

octco too CX 1 too 


Cfot) 


ocfl-cc) 


{■ 


X i\g\ Ifl (|V I ) + Igl Ifl flu I ) 

CX CO XX oo 00 <X XX CO 




CM 


OC (l-oc) 


igi {ifi igi ( i v i r + ifi igi flu i ) 

1 ° 1 00 1 00 OC XX CO OC ° CD XX 00 


■}] 


So this will imply, 


Ifl 


(inf Igl 2 ) -1 
(inf |0' I ) c 


[ 1*1. 10' 1*1. < ie ' 'KK 


U<XJ s \ 

- - - 1*1. + 1*1. 'S'“ 

OC (l-oc) v 

{ifl Igl ( I V I )" + Ifl Igl flu | )* 1 1 

1 00 OC XX CO OC CO XX 00 J J 


(2.6.31) 


Here, equation (2.6.18) will be repeated for the sake of completeness. 
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Ifl 


00 


2S 


inf Igl 



IS' - * t l 


CO 


+ 


- (cx+2) 

2 


C (2«) 

3 

<x 


{l f l« Is'. * lfl « lsl 


( i V 

XX 



or, 


i?i 


00 


+ 


[inf Igl ] [ )e' - * t \„ 

t“ {in. Igl. CI..I.)- - Ifl. 1st. } 

« Cl-oc) V 

(2.6.32) 


Therefore, adding (2.6.31) and (2.6.32) the oc-norm of f can be found out. 

II f 1 = Ifl + I f I 

II II OC 00 oc 

s [inf Igl ] ‘ [le' - (SJ. 

C(«) f \ ■ 

+ t“ 4 Ifl Igl (lu | )“ + Ifl igl (|v | } 

. . . I OC ° 00 XX CD CO « XX « 

OC (1-oc) v j J 


(inf Igl 2 ) -1 
(inf |e'|)“ 


Igl le' - * I * Igl (|v I ) W le' - * I 

00 t OC OC 1 00 tv 


C(cc) 


oc (1-oc) 


( 'V“ ti\i.)“ lei.) 


{ If| . Igl « <lv J» ) ° < * Ifl. I el. fl^l„>“ } 


[ 


£ I (inf Igl)" 1 + 


(inf Igl 2 )" 1 

— Igl 

(inf |8' |) 60 
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+ (inf Igl 2 )" 1 |g| 


(Ivl ) 


t to 

f Ifl' I ■ 


(inf |0' |) 


X | 0' - 0 


t"« 


c(°o t 

« (l-c<) l 


t~ (inf Igl)" 1 + 


(inf |g| 2 ) -1 


(inf |0' | ) 


Igl. 


+ |g| I 1 , I 

OC 1 


« (l \ n „ i 

- x (inf |g| ) J. 


(inf |0' |) 


*{ 


igi du i ) + igi dv i r\ | f n 

CO XX CO <x XX 00 [ 11 ,l oc 


'1 


i M, I 6' - III + M„ I f 1 

1 » *t « 2 11 11 « 


(2.6.33) 


(2.6.33) implies that || f is bounded above by || f || k . So if || f | £ E, we 

can have || f £ E provided M.^ and are small quantities. Here one needs to 
show, (|u | ) <x , ( | v | )** and (|v | / (inf |0'|)) <x are small so that M. , 


XX to 


XX to 


are small enough quantities. 

To estimate u , consider the relation 

XX 


„t JO 


u = 0 + f f K(x,y,t-x) f(u(y,x)) g(v(y,x)) dy dx 

■^o ^ 0 


,t „oo 


u = 0 + f [ K (x,y,t-x) f (u(y,x) ) g(v(y,x)) dy dx 

XX XX J„ J„ XX 


0 •'o 


= 0 + f f K (x,y, t-x) [ f (u(y,x) ) g(v(y,x)) - 

xx ^ 0 ^ 0 xx L 

f(u(x,x)) g(v(x,x))j dy dx 

= 0^ + J J K^tx.y.t-x) £ g(v(y,x) ) j f(u(y,x)) - f(u(x,x))| 
+ f (u(x,x ) ) | g(v(y,x) ) - g(v(x,x))j- j dy dx 
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Now taking the estimates on both the sides and expanding u and v with 
respect to x and from (A5), it becomes 


lu 


I * f f IK (x,y,t-r)i { lgl w lfl w ly-xl“ CiuJJ* 
XX J o J 0 ' 


♦ i f L I*I« 


dy dx 




,t „00 


[ f IK (x,y,t-x) I ly-xl* dy dr 


* I lgl w if i„ (lu.U a + IgL l f L flv |j a 


\ c («) 


OC X 00 


OC 00 X oo 


3 ' t <x/2 


In veiw of Appendix A (A3) we can similarly estimate Ju | 

X CO 


(2.6.34) 


|u - \p | s f f IK (x,y, t-x) | |f (u(y,t) ) | |g(v(y,x) ) | dy dx 
x x J 0 J o x 


- Ifl w lgl w J J lK x (x,y,t-x) | dy dx 


s C CO) VT |f| Igl 

c 00 00 


(2.6.35) 


Substituting (2.6.35) into (2.6.34), we therefore have 


\n - \p \ 

XX XX 


* I 1 *'. If i« 'iVJ” + I*I„ i f L 

' oc 


l oc/2 


{igl. Ifl. (i* x i * c 2 (0) vrin. Igl.]” * Ifl. Igi ra (Iv.i.r } 


i CM 

oc I 3 


(2.6.36) 
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Let <f> be another comparison function corresponding to v, which satisfies 
the differential equation in a semi-infinite domain 


d 4> - d <f> - 0 

t XX 


, X > 0, t > 0; 


0(x,O) = v q (x) , x > 0; 

= g 2 (t) , t > 0 ; 


(2.6.37) 


So, if q is represented by v - <f> , then q becomes 


q = v - $ = J J K(x,y,t-x) | - f(u(y,x)) g(v(y,x))j dy dx (2.6.38) 


Differentiating (2.6.38) with respect to x will give rise 


v x = * x + J J K x (x,y,t-x) | - f(u(y,x)) g(v(y,x)) | dy dx (2.6.39) 


Now taking modulus on both the sides, 


|v - <f> | s | f I Igl f f IK (x,y, t-x) | dy dx 

x x oo oo J o J q X 


~ |f L l«L C , (0) 

oo oo 2 


(2.6.40) 


putting (2.6.40) into (2.6.39), can be written as 


|u - 0 | a | |g| If I f 10 I + C COD VT If I Igl ] 

'xx xxlcooc^x 2 CO 00 J 


* l f l„ Igl. ( l*,l - c,(0) VT Igl. Ifl.) } 


« \ c(«) 

_3 t a/2 


(2.6.41) 


This implies that 
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u - \f> — » 0 

XX XX 


like Oft 0 ' 72 } 


if I f l« 5 E and UL s E 


=» 


lu | a i^& j 


XX co 


XX CO 


k 


( 2 . 6 . 42 ) 


Similarly, v can be estimated. 

XX 

Since, 


v = d> + 

XX XX 



K (x,y,t-x) 

XX 


- f (u(y,x) ) g(v£y , x ) ) 

+ f (u(x,x) ) g(v(x,x) ) 


j- dy dx 


So 

|v xx " ^xx 1 * J j lK xx (x,y,t-x){ j lg(v(y,x))Hf(u(x,x)) - f(u(y,x))| 

J 0 J 0 V 

+ If (u(x,x)) | ig(v(x,x)) - g(v(y,x) ) | 1 dy dx 


s { lgl . lf| « lf L i*i, dv.ij”} 

f f lK(x,y,t-x)| |y-x|“ dy dx 
J 0 J 0 

S { If '« - lfl„ l8l„ (lv x l„)“ } 

* oc 

s { lgl » lf '« ( i#,i * <y 0 ) vr ifi„ i s i„ )“ 
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|f L '*'« ( * c 2 <°> igi„)“} 


<* A C («) 

t a/2 (2.6.43) 


So v - <p — > 0 like 0(t°' /2 ) if Iff] < E 

XX XX I II oc 


* |v I s II I 2 |v" I 

XX 00 XX 00 0 CD 


(2.6.44) 


Now we turn to estimate v 


Since 


= <P + J J K(x,y,t-x) | - f (u(y,x) ) g(v(y,x)) j dy dx 


=* v = d> + -- - - - 

t v t at 


£ 


F(t,x) dx 


= <t>. + F(t, t) + [ 

J n 


aF(t,x) 

at 


dx 


(2.6.45) 


Where 


00 f v 

F(t,x) = J K(x,y,t-x) | - f(u(y,x)) g(v(y,x)) j dy dx 

Using the Causal Property of the Fundamental Solution K (STAKGOLD [111]) 
we say F(t,t) = 0. Since K(x,y,t) = 0 for t £ 0, or in other words, 
K(x,y,t-x) = 0 for t s x, differentiating F(t,x). 

F(t,x) = J K(x,y,t-x) | - f(u(y,x)) g(v(y,x)) j- dy dx 


We get 


3F ( t , x ) 

at 


= J K t (x,y,t-x) | - f(u(y,x)) g(v(y,x)) | dy dx 


Using the identity property (2.6.7), the above expression can again be 


written as 
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SF(t,x) 

St 


r 

* r\ 


Taking modulus on 


K (x,y,t-x) 


g(v(y,x) ) 


{ 


+ f £ u(x,t) ) 


both the sides, one gets 


f (u{x,x) ) 

g(v(x,x) ) 


f (u(y,x) } 


g( v (y, x ) ) 



dy dt 


SF(t,x) 

St 




Is!. 


If I (lu | ) 

oc X oo 


If 1 Igl ( iv 1 ) 

oo oc x CO 


\ r °C 

“I 

J c 


iK t (x,y,t-x)! ly-xl dy 


Now, integrating both the sides w.r.t. to x and considering the limits 
from 0 to t, the right hand side will be given the value of the domain 
constant as calculated in Appendix A, (A9). 


I 


SF 

St 


(t,x) 


dx 


s { 's'. 


Ifl (lu I ) + If I Igl <|v 

oc X 00 00 oc 


I )* ) X 

X CO j 


C («) 

s t «/2 


As t — » 0, since all the other terms are bounded above, the L. H. S. goes to 
zero. 

=> l v t ~ ^ t l — * 0 as t — » 0. 

Consequently, 


IV 1 1 

oc 

r !♦ 1 1 

OC 

|v"| ] 

t 

< 

XX 

££ 

0 

inf |6'(t)| 

. inf |G' (t) | 

_ inf |G' (0)! _ 


Now, in view of the inequalities (2.6.42), (2.6.44) and (2.6.46) we can 
say that these estimates can be made very small by choosing flat initial data 
fu^'l c 1 and |v' ' | « l. 
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Thus by taking t small enough and initial data flat enough , M 1 and of 
(2.6.33) can be made sufficiently small. 

I f L - M 1 I 6 ' - + »2 I f l« 


- M i I e ' - * 


CM (inf | g | 2 ) _1 


« (!-«) (inf \Q ' | ) c 


Igl 


{ 


Igl (lu I ) + |g 

CO XX 00 


I (|V | )* \ 

<x XX 00 J 


C(«) , , 

a «! I®' - # t l« + - , { 1 * 1 . "V l.>“ * 1 * 1 . } E 

oe (1 — oe J V J 


£ E 


(2.6.47) 


In view of the smallness of , a choice of E can be found so that 
£ E if flffl^ £ E. 

Hence T : B c c“ > B c C* . I 

e e e ■ 

Finally, we note that T is a continuous map from C * to C°\ This follows 
from the continuous dependence of the solution "u" of the direct problem on 
the right hand side function f. 


LEMMA 2 : The set U = j f : || f fl w £ E, | f | a £ E J- is closed set in 


C* 


forO<0<oc<l. 


PROOF : Clearly, U is a closed set in C . 


Let 


M 


be a sequence in U such that f — » f in for 0 < « 


We claim that f € U. 


Since f 


f in 





* if - f| — - > 0 

" n *® 


so I f I. » I { X * l f „ - X 


Hl.^ ( I g„ * l f „ - f L ) 3 E 


Further , 


jf(x) - f(y)J 

|x - yr 


11m (f (x) - f (y ) 

n n 


x - y 


s lim 


|f (x) - f Cy } 

f n n 


* f s E 

oc 


5 lim E = E 


Hence U is closed in C^. g 


THEOREM 2 (EXISTENCE) : Assume that T Q is a seif mapping from U into 
itself, where U = | f : JfJ^^E, J f s E j- . Then has a fixed point. 

PROOF : The set U is a closed, convex set in C°\ By the above lemma U Is 
also closed and convex set in . Further, since U is a bounded set in c“, 

T e U is P reco ®pact in Hence T g is a completely continuous operator taking 
: U c cP >UcC^ . Therefore, by SCHAUDER’ s fixed point theorem 
(FRIEDMANN, p. 189 [49]), T g has a fixed point in C^. g 

REMARK : 

However, the above theorem does not prove the existence of a limit of the 
iterated sequence jf j obtained above. In the next chapter we shall prove a 

— tract ' ion ma P pi ~ ng Pr °Perty of the T g map under additional hypotheses 



which will show the convergence of the iterates and there by establish 

the convergence of the iterative scheme outlined above. 

2.7 CONCLUSIONS : 

The basic problem dealt with in this chapter is to prove the existence of 
a non-linear source function of a system of two Reaction-Diffusion equations 
in a semi-infinite spatial domain. At the out set the original problem is 
reformulated as a fixed point problem. An iterative scheme is outlined via a 
self mapping operator in a bounded set in the Holder space. Using the 
properties of Green’s function many domain estimates have been calculated 
explicitly. Finally making use of the SCHAUDER’ s fixed point theorem, the 
existence of a fixed point of the solution of the original source problem is 


established. 




CHAPTER - til 


UNIQUENESS OF THE SOLUTION OF THE REACTION-DIFFUSION SOURCE 
PROBLEM AND CONVERGENCE OF AN ITERATIVE SCHEME 


3.1 INTRODUCTION : 

The existence of the solution of the Reaction-Diffusion source problem 
with superposed Dirichlet data is established in the previous chapter. In this 
chapter, the uniqueness of the solution and the convergence of the proposed 
iterative scheme outlined in the previous chapter is established. We follow 
the notations of chapter II and assume the validity of the conditions stated 
in that chapter. 

The existence theorem ensures existence of a 0-fixed point of the 

operator T . Therefore, there is an "f" such that 
0 

f(e(t)) = T 0 [f](t) (3.1.1) 

In order to prove the uniqueness, we proceed as follows : 

Denote p (0 ( t ) ) =T [p](t) (3.1.2) 

0 

We show that 

|f - P|„ * A |f - Pl. with A* (0.1) 

Consequently, there can not be any other 0-fixed point "p M of , for in that 




case 


Jf - p|^ s A |f - p{ o which is a clear contradiction. 


Below we show that Tg , indeed satisfies the following contraction 
property i. e., there exists A satisfying 0 a A < 1 such that 

Jf(e) - p(9)| a = - Tg(p) | w 


5 A |f - pf a (3.1.3) 

The main result of this chapter is stated in section 2. The sup and 
semi-norm are computed in sections 3 and 4 respectively to establish the 
central result. Finally a brief summary with possible future studies are posed 

in section 5. 

3.2 THE THEOREM : 


It is assumed that the overposed boundary function © belongs to the 
class of Lipschtz functions with some additional smoothness assumptions to be 
stated precisely later. 

For any uniformly Lipschltz function f and any g € C* , f o g € C*. 
In fact 


if o 





(3.2.1) 


We saw that for any Lipschitz class of functions on a bounded domain, 
preserves C regularity under composition. This is not sufficient to show 
uniqueness. We say that a function f e Lip has “property S" if 

PROPERTY S : Given function u, v e C*. for 0 < « < 1 , the mapping u — » f(u) 

oc 

is C , that is there exist C < « and 

s 

|f(u) - f (v) | < c jju - v | (3.2.2) 

THEOREM . If T 0 operator for the the overposed boundary value problem 
possesses a ©-fixed point f and both f , g of R-D equations satisfy property S, 



then f is the unique fixed point. Furthermore, the operator T is contractive 

0 

in d* for « < 1 about f. i. e. , 

i T e [f] " T 0 [p] l« s A ll f - P|« with 0 s A < l. 

3.3 SUP-NORM ESTIMATION : 


Let u = u(x,t;f,g); be the solution of the differential equation 



U t ~ u xx = f(u) g(v) ’ 

x > 0, 

t > 0; 

(3.3.1a) 


u (0,t) = 0 , 

X 

t > 0; 


(3.3.1b) 


u(x,0) = u q (x) , 

x > 0; 


(3.3.1c) 

satisfying 

the overposed boundary condition, 





u(0,t) = e(t) , 

t > 0; 


(3.3.2) 

where v is 

the solution of the differential equation 




v - v = - f(u) g(v) , 

t XX 

x > 0, 

t > 0; 

(3.3.3a) 


v (0, t) = 0 , 

X 

t > 0; 


(3.3.3b) 


v(x, 0) = v q (x) , 

x > 0; 


(3.3.3c) 

Let p 

correspond to another set of solutions of 

the over-posed 

problem 

with w = w(x,t;p,g) and v (x,t;p,g) satisfying 

the following : 



w - w = p(w) g(v) , 

t XX 

x > 0, 

t > 0; 

(3.3.4a) 


w (0, t ) = 0 , 

X 

t > 0; 


(3.3.4b) 


w(x,0) = u q (x) , 

x > 0; 


(3.3.4c) 


w(0,t) = 0(t) , 

t > 0; 


(3.3. 4d) 

and 

v - v = - p(w) g(v) , 

t XX 

X > 0, 

t > 0; 

(3.3.5a) 


v (0, t) = 0 , 

X 

t > 0; 


(3.3.5b) 


v(x,0) = v q (x) , 

X > 0; 


(3.3.5c) 



The expression for the map T 0 from (2.6. 13) is given explicitly by 


f(e(t)> = T 0 lf](t) 


f t r ® 

K (0,y,t-tJ f (u(y,t 1 ) g(v{y.x) ) dy dt 


©' (t) - 0 t (O.t) 


'0 0 


g(v{0, t, ) 5 


(3.3.6) 


and analogously. 


pCe(t) ) = T 0 [ P ] (t) 


t 05 


e' (t) - \}t (0, t) - K (0,y, t-T) p(w(y,T ) ) g(v(y,T) ) dy dr 

V A V ft 


'0 "0 


g(v(0, 1 3 ) 


Here \ji is the comparison function (repeated here for convenience). 


(3.3.7) 


* CO, t ) = 0 , . 

X 

^(x.O) = u q (x) , 


X > 0, t > 0; 

t > 0; (3.3.8) 

x > 0; 


First of all we set out to estimate If - pf , for which we use the 
expressions in (3.3.6) and (3.3.7). In order to avoid repeat itive expressions 
we divide the expression for I below into several parts and estimate each 
part separately. 


I Q = f(9(t)) g(v(0, t) ) - p(9(t)) g(v(0, t) ) 


*t JX> 


~ I o J o K xx^ 0,y,t_T ^ P (w( y* T) ^ g(v(y,x)) - f (u(y,x) ) g(v(y,x}) 


dy dx 


Hence, 


( iC0(t)) - p(e(t)) | g(v(0, t) ) = p(e(t)) | g(v(0, t) ) - g(v(0, t) ) J 



J J K^CO.y.t-c) p(w(y,x) ) - f(w(y,x)) j g(v(y,x)) dy dx 
J J K xx (0,y,t-x) f (w(y,x) ) - f(u(y,x)) j g(v(y,x)) dy dx 
J J K xx (0,y,t-x) T g(v(y,x) ) - g(v(y,x)) 1 f(u(y,x)) dy dx 


= I + I + I + I 

12 3 4 


(3.3.9) 


In order to find an estimate for L. H. S. of equation (3.3.9), 1^ , for i 
= 1, 2, 3 and 4 are computed separately. In the expressions I ,1 and I , f 

13 4 

- p does not enter explicitly. But the arguments w, u, v and v can be 
expressed in terms of f and p as they are solutions of the corresponding 
non-homogeneous differential equations. To calculate the max-norm of L. H. S. 
the idea is to express I , I and I in the form of |w - u | and |v - v |. 

13 4 xx xx 

Since the latter quantities involve first space derivative of the Neumann 
function, the domain estimates (Appendix A) will not blow up. Later |u x - W x l 
and |v - v | are expressed through lower order terms and their weighted 

X X 

values are calculated. These modifications are needed because in I and I the 

3 4 

kernel is the second derivative of Neumann function. 


Since u and w are the solutions of the non-homogeneous differential 
equations (3.3.1) with (3.3.2) and (3.3.4) respectively, so the difference is 


given by (see CANNON [26]) 


w(x,t) - u(x,t) = 


CENTR4L LIBRARY 

I I. T.. KANPUR 


4m Afe A 


J J K(x,y,t-x) £p(w(y,x)) g(v(y,x)) - f(u(y,x)) g(v(y,x)) j dy dx (3.3.10) 


Similarly, the difference between the solutions of (3.3.3) and (3.3.5) is 



v(x,t) - v(x,t) 


C f K(x,y,t-x) r f(u(y,T)} g(v(y,x)) - ptwjy.x)) gfv(y.x)} 

■*0 '• 


dy dx (3.3.11) 


From (3.3.10) and (3.3.11) it Is evident that ufx, t) - w(x.t) * ^(x,t) - 
v(x, t) and hence (w - u| = |v - v| and the sane is true for their higher 

order derivatives. 

Since the domain integrals (Appendix A) have a form, where kernel is 
multiplied by |x - y|°\ we require to use the identity property (2.6.7) 
defined in chapter II, to bring the expressions in the proper form. 

The computation of if - pi involves several integral estimates. 

CO 

Therefore, it is subdivided into several steps. 


STEP 1 : ESTIMATION OF 

Let us compute 1^, 
it can be written as 

- 1 ] f a 

r t .00 

+ I I K xx (0,y,t-T) 
J 0 J 0 


I CO < « < I). 

0 

first. Assume f, g € C°\ Using the identity property, 


p(w(y,x) ) - f (w(y,x) ) 

- p(w(0,x) ) + f (w(0,x5) 
£ p(w(0,x) ) - f (w(0,x)) 
[ g(v(y,x)) - g(v(0,x) ) 


g(v(y,x) ) dy dx 


x 


dy dx 


Hence , 

ft CO 

£ J J * K xx^ 0 ,y, t~x) I IP “ f I lw(y,x) - wCO.x)!* jgi dy dx 

0 0 ** co 

+ I 0 J 0 |p - fJ w jv(y,x) - v(0,x)| oc jgj^ dy dx 
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Expand w and v by Taylor’s expansion, considering centre at 0 and radius 

as y. From the assumptions that g = 0 = g and Appendix A (A6), |I | becomes, 

1 2 2 

|I 2 I £ ISL fP ~ f l« (|W X x l a, ) °'J J |K xx (0 ’ y ’ t ' T) ll(y 2 /2)| OC dy dT 

J 0 J 0 

+ Ip - f l ro J o J o | K xx ( 0 , y ,t-T)||(y 2/2) r dy dr 

Hence , 

III s IlgL IP ~ fl a (lw I . )“ + 1 g I Ip - f| C I v | ) K 1 C (2<x) 2" (a+2) — 

2 I co ocxxoo oc oo xx oo I 3 oc 


Noting that o<~semi norm and sup norm are individually less than oc-norm, 

we can write estimate I as 

2 


* \ ** Ip - f L 

Where A = I Igl (|w | )“ + |g| |v | a 1 - 

1 I 03 XX oo CX XX 03 j 

chapter II (in the box). For the estimate of I 3> 


C 3 (2o0 

CX 


p - f 


(3.3.12) 

C(oc) 

, C(.) is defined in 

«(1 - oc) 

we use the "property S" . 


|u(£,t) - u(r) , t ) | 

Since, by definition |u(.,t)| = sup 

K i? - v r 

Put £ = y, 7 ) = 0 and u(.,t) = f(w(.,t)) - f (u( . , t ) ) 

then | f (w( . , t ) ) - f(u(.,t))l t = 

sup |f(w(y,t)) - f(u(y,t)) - f(w(0,t)) + f(u(0,t))| / |y| 
y*o 

Since f € Lip class and satisfies property S, so making use of the 
identity property, I 3 can be expressed in the following form. 
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u K^O.y.t-x) g(v(0,x)) - g(v(0,x)) J x 

^ f (u(y,x) ) - f (u(O.r) ) j dy dx 


= I + I 

4,1 4,2 


Taking into account |w(x,t) - u(x,t)| = |v(x,t) - v(x,t)|, the 
w(0,t) - u(0,t) = 0 implies that I =0. Thus, 

4, 2 


ill - III £ -4- D |f| C (2) t |v - v | 

4 4,1 4s c»3 xx xx i 


s A t |u - w | 

3 xx xx co 


For g e C or Lip class the above argument is also true for I 


I s o 

l 


Thus "step 1" leads to the following estimate in view of 
(3.3.14), (3.3 16) and (3.3.17). 


If (0(t) ) - p(0(t))| |g(v(0, t) ) | £ III + |I | + |I | + |I | 

12 3 4 


or 


I f (e ( t ) ) - pCeCt ) ) | ^ inf |g| x 


| inf |g| j 


t lip - f II +At|u -w| + At 

11 oc 2 xx XX CO 3 


h 

If (©(t) ) - p(e(t))| s 
[ A, t“ Ip - fl„ * A t t lu^ - ] 


inf |g| x 


|u - w | ] 

XX XX 00 J 


Where A = A + A 

4 2 3 


STEP 2 : ESTIMATION OF w - u : 

XX XX 


(3.3.15) 

fact that 

(3.3. 16) 

. Hence 

(3.3.17) 

(3.3.12), 


(3.3.18) 


Differentiating (3.3.10) twice with respect to x, one gets 
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w (x,t) - u (x,t) 

XX xx 

= f f K (x,y, t-x) p(w(y,x)) g(v(y,x)) - f(u(y.xj) g(v{y,x)) 

= r K (x,y,t-x) [ p(w(y,x)) - f(w(y,x}) g(v(y , x 3 3 dy dx 

^ 0 ^ o xx 

t oo r 1 ^ 

+ f f K (x,y, t-x) fCw(y.T)) - f(uCy,x)} g(v(y,x)) dy dr 

J o J o xx 

+ J J K xx (x,y,t-x) J g(v(y,x) ) - g(v(y,x)3 J f (u(y,x)3 dy dx 


dy dx 


= I + I + I 

5 6 7 


(3.3.19) 


We note that, these expressions are similar to I I and 1^ respectively 
(except that x stands in place of zero in K ). Now we provide estimates for 
these expressions in terms of the first derivative of u and w. Applying the 
identity property (2.6.7) with the assumptions f , g « C** and Appendix A (A5). 

t 00 

|I 5 I £ J J IK joe (x,y,t-r)| | p(w(y^T) ) - f(w(y,T)) 

- p(w(x , t ) ) + f (w(x,t) ) | lg(v(y,x ) ) | dy dr 
f f IK (x,y, t-x) | | p(w(x,x)) - f(w(x,x))| x 

J o J ^ XX 

i g(v(y,x)) - g(v(x,x))i dy dx 


0 ”0 


C {«) 

s Igl Ip - fl (lv I )“ t ^ 2 


« X 00 


C (oc) 

+ Igl Ip - f I Civ I )“ — t K/2 


00 x 00 


\ IP - t 0 " 2 ♦ A 


l oc/2 


* A 7 Ip - f| a t 


ck/2 


( 3 . 3 . 20 ) 


Since f e Lip class and satisfies the property S, I is estimated as follows : 
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t, 00 _ 

1 6 = J o J K xx (x *y> t_x) [ f ( w (y.T)) - f(u(y,T)) 

- f (w(x, r) ) + f(u(x,T)) ] g(v(y,T)) dy dT 

+ J J K xx (X,y,t-T) [ f(w(x > T ^ - f (u(x, t) ) J x 

£ g(v(y,T)) - g(v(x,T)> J d y dT 


= I + I 
6,1 6,2 

As f satisfies the property S, I is estimated as (3.3.13). Indeed, 

6, 1 

J*> 


(3.3.21) 


U 6>1 I s J J lK xx (x,y,t-x)| |f(w) -f(u)l i Igl^ lx - yl d y dT 

- c f f IK (x,y, t-r) | |w(.,t) - u(.,t)I l x ~ y * * g L dy dx 

J 0 J o xx 

t 00 

- C J g L l u - w I f f IK (X.y.t-T)l |x ' yl d y dT 

s 00 X xooJ J XX 


A |u - w | t 

8 x x co 


1/2 


(3.3.22) 


Further, since f, g € Lip class with constants C and D (say) respectively, we 
have 


- J J | |K xx (x, y ,t- T )| l f | lip |(w(x,r)) - (u(x,t))I lgl Hp 

|x - yl H^IJ | dy dr 

- C |w - u| D ( | v | ) f f |K (x,y,t-r)l l x " y l dy dT 

" * “ J „ J „ xx 


1/2 


s A |w - u| t 
9 oo 


Hence from (3.3.22) and (3.3.23) 


(3.3.23) 


III (A u - w I + A |w - u| ) t 

6 8 x x oo 9 oo 


1/2 


(3.3.24) 
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I is estimated exactly like I & . So we can write 

1 1 | s A iu - w 1 t 1/2 + A |w - u I t 1 

*7 io x n w 


(3.3.25) 


In view of (3.3.20), (3. 3. 24) and (3.3.25) Step 2 leads to the estimate 


I w - u | a III + III + III 

XX XX ® 5 6 1 


0 .«/2 j, , , *1/2 ^ . , , .1/2 

t + A !u - W 1 t + A I W - U j t 

“<x 12 x x co 13 e& 


(3.3.26) 


Where A = A + A and A = A + A 

12 8 10 13 9 11 


STEP 3 : ESTIMATION OF w - n 

X X 


This step is devoted to calculate |u ■ w | . Since 

X X 00 

w (x, t) - u (x,t) 

X X 

= | J K x (x,y,t-x) ^ p(w(y,x) ) -f(w(y,x)) g(v(y,r)) dy dr 

+ I 0 I K x (x,y,t ' T) [ f(w(y,x)) - f (u(y,x) ) g(v(y,T)) dy dt 

+ -L I 0 K x (x,y,t ~ T) [ g(v(y>x)) “ 8 (v£y ' T)) f (u(y,x) ) dy dr 


= 1+1 + I 
8 9 10 


(3.3.27) 


From (A3) of Appendix A, we get 


f t r® 

'V s Jo J 0 |K x (x,y,t " T)| lp(w(y,x)) - f(w(y,T))j |g(v(y,T))|dy dr 
S ,p " f Ll g Lj J lK x (x,y,t-T)| dy dr 


0 *0 


1 A „ Ip - f L t’ 


a u - c a (oi 


'*1.) 


(3.3.28) 
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For f e Lip class 

,t -CO 


UJ 5 f f IK (x,y,t-x)| |f(w(y,r)) - f(u(y,r))| lg(v(y,x) ) | dy dx 
9 J o J o 

C f f IK (x,y, t-x) | | w(y,x) 

j A ^ A ^ 


-u(y,x)| |g(v(y,x))| dy dx 


o 'o 


C |g| C (0) |w - u| t 

co 2 co 


1/2 


2 A |w - u| t 

15 co 


1/2 


(3.3.29) 


Now for g e Lip class, I is estimated similar to I and the result is 

10 9 


given by 


I | ^ D |f 1 C (0) | v “ v | t 

10 oo 2 co 


1/2 


^ A |w - u| t 

16 oo 


1/2 


(3.3.30) 


Therefore , in view of (3.3.28), (3.3.29) and (3.3.30) "step 3" gives the 


following estimate for w - u in terms of w - u : 

X X 


Iw - u | s |I| +|I | + |I| 

x x oo 8 9 10 


£ A lip ~ f || t 1/2 + A |w - u| 

14 11 11 ot 17 co 


l 1/2 


(3.3.31) 


Where A = A + A 

17 15 16 


STEP 4 : ESTIMATION OF |w - u| . 

1 1 00 

Finally, we estimate |u - w| , where the kernel is the Neumann function. 
Since we have calculated in Appendix A (Al) that the domain integral 
corresponding to K(x,y,t-T), is a finite quantity, the estimation procedure 
will be easier. Expressing w - u in its integral form from (3.3.10) 
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w(x, t) - u(x,t) 

.t r ® 

= K(x,y,t-x) 

Q " 0 

t © 

+ f f K(x,y, t-x } 

•'o ^ 0 

t © 

+ [ f K(x,y, t-T ) 

Jn J ft 


'0 "0 


= I + I + I 

11 12 13 


p(w(y , t ) 5 - ffwiy.r)) 
f (w(y, x ) ) - f (u(y, T ) ) 
g{v(y,x) ) - g(v(y,T)) 


g ( v ( y , x ) ) dy dx 


g ( v { y , t ) ) dy dr 


f (uty.x) ) dy dr 


(3,3.32) 


Let f, g € C , then 


|I I s f f jK(x,y,t-T) j |p(w(y,T) ) - f (wfy.x) 5 j g(vfy.x) ) dy dx 

** 0 J Q 

- |P ~ f L |g| ffi f f |K(x,y, t-x) j dy dx 


s i*i. c,«» l p - f L 1 


=* |I | s A 
11 18 


For f e Lip with constant C, I is estimated as 


(3.3.33) 


ft CO 

1 I i 2 l - lK(x,y,t-r) | |f (w(y,T) ) - f(u(y,T))l lg(v(y,r))| dy dr 

J o J o 


* C |g| C (0) |w - u| t 

tD 1 CO 

s A w |w - ul^ t (3.3.34) 


and finally, for g e Lip class with constant D, we have 

-t -co 

^la 1 S J J lK(x,y, t-r) | lg(v(y, T )) - g(v(y, T ))| |f(u(y,x))| dy dx 
0 v 0 

S D |f I C (0) |v - v| t 
« 1 00 

s |w - u| t 

20 oo 


(3.3.35) 
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So the following estimate follows from "step 4". 


|u - w| s |I | + |I | + |I | 

' oo 11 12 13 


* A ,a l p - f l« 1 * A 2 i ‘ 


Where A = A + A 

21 19 20 


=> |w - uj ^ 
00 


18 


1 - A t 
21 


IP-f| K t 


(3.3.36) 


Where A 


22 


18 


1 - A t 
21 


STEP 5 : ESTIMATION OF THE SUP-NORM OF f - p 


Now substitute the inequality (3.3.36) in (3.3.31) to get 


|w - u | s A |p - f J t 1/2 + A A | p — f II t 3/2 

x x 00 14 “ “<X 17 22 "(X 


s (a 


+ A A t Ip - f I t 

14 17 22 11 ^ 11 « 


1/2 


s A 23 Ip - t 


1/2 


(3.3.37) 


Using the estimates (3.3.36) and (3.3.37), (3.3.26) can be written as 


w - u | s A | p — f II t a/2 + A |w - u | t 1 '* + A |w - u| t 

xx xx oo 7 H 11 oc 12 x x 13 oo 


1/2 


L 1/2 


* \ |p - ^ + A U A 23 I p - f l« * 4 a .3 A 22 Ip ' f L 1 


3/2 




, a a a.I-Ox/ 2) . A A .(3-o<)/2 ) ,, .«/2 

4 A - A - 1 4 A 13 A 22 ‘ J IP f V ‘ 


12 23 


* A 2t Ip ' f L 1 


oc/2 


(3.3.38) 


Finally substituting (3.3.38) in (3.3.18), we get 
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If (e(t> ) - p(e(t))i 


s ; inf Igl 


A Ip - f[ t" + A I u - w | t 

1 * * » 4 *3sr %% m 


s inf |gi 


-5 


A , Ip - f ! a 1 * A 4 Ip ~ f l t 




a inf fgi 


A + A A t 

1 4 24 


I - (a/2 1 


b - ff, t* 


s A 2 s Ip " f L t0 


£3.3.39) 


Where = | inf fgj j 


-1 r 


A + A A t 

1 4 24 


i- to«/2) 


3.4 ESTIMATION OF HOLDER SEMINORM 


For completion of «-norm, «-semi norm of £f - p) is to be estimated. 
Consider t > t Then )f - p| is written as 

12 oc 


l f “ Pl^ - sup 

0 ( t ) *0 ( t ) 
1 2 


f( 0 (t t )) - p(e(t 1 3 ) 


Hi 


f(9(tj) - p(G(t )) 
2 2 


left ) - ©uj I 

l 2 


(3.4.1) 

Since f is a e-fixed point. So in view of (2.6.13) it is expanded 
corresponding to t j and t g explicitly as 

f(e(t )) g(v(0,t )} « e'(t ) - t/r (0 t ) 

1 1 t 1 

J o J o V* t ^ f£u(y,x)) g(v(y,T) ) dy dT (3.4.2) 

f (e(t 2 )) 8(v(o,t a )) = e'(t ) - ^ ( 0 ,t ) 

2 t 2 
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“ J o 2 J o K X x (0,y ’ V x) f(u( y« T) ) g(v(y,x) ) dy dx (3.4.3) 
Similarly we can express for p 
p(e ( t i ) ) g(v(o, ) ) = e / (t i ) - co , ) 

" K v X (0>y ’ t r Tj P (w (y. T )) g(v(y,x)) dy dx (3.4.4) 

J o J o 

p(e(t )) g(v(o, t ) ) = 0 ' (t ) - i/i (o,t ) 

2 2 2 t 2 

K xx (0 ’ y,t 2 _T) P (w( y> T) ^ g(v(y,x) ) dy dx (3.4.5) 

Subtracting equations (3.4.3) and (3.4.4) from (3.4.2) and then adding 
equation (3.4.5). The L. H. S expression of the equation becomes 



f(e(t i )) g(v(0,t i )) - p(e(t a )) g(v(0,t )) 


-[ 


f(e(t )) g(v(o, t ) ) - p(e(t )) g(v(o, t )) 

2 2*2 


] 


-[ 


g(v(o,t i )) + pfeu^) 

g(v(0,t i )) - g(v(0,t i )) j 

g(v(0,t )) - p (0 ( t ) ) 

2 2 

g(v(0,t 2 )) - g(v(0,t 2 )) J 


([ 


f(e(t )) - p(e(t 




f(0(tj) - p(0(tj) j j g(v(0,t,)) 


+ pfOt^)) £ g(v(0,t,)) - g(v(0 


•V 1 ] 


- g(v(0,t 2 )) - gWO^)) j J f(0(t 2 )) - i(0(t 2 )) j 


“ P(0(t o )) £ g(v(0, t^) ) - g(v(0,tj) 


] 


(3.4.6) 


Similarly the R. H. S. becomes 
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f'J 'W 0 ^V‘ ) 

J o J o 

t -CO 

2 K (0,y, t -x) 

J xx 2 

0 ^0 


Hence from (3.4.6) and 



p(wfy.x) ) g(v(y, t ) ) 
p(w£y, x ) 5 g(v(y.r ) 5 

(3.4.7) one gets 


- f ( u l y , x ) ) g ( V ( y , X ) ) 

- f (u(y,t ) ) g ( v ( y , t ) ) 


dy dt 

dy dx 
(3.4.7) 


( [ f(e(t . 


)) - p(e(t )) 


f (e(t )) - p(8(t )) 

4 


g(v(0. t ^ ) ) 


+ p(6(t )) g(v(0, t ) ) - g(v(0, t ) ) 


= [ g(v(0,t 2 » - gtvto.tj)) 


f (Q( t ) ) - p(0(t )) 
2 2 


+ p(e(t 2 )) g(v(o,t 2 )) - g(v(o,t 2 )) 


..t ,co 


f 1 [ K (0,y,t -x) p(w(y,x) ) g(v(y.x)) - f(u(y,x)) g(v(y,x}} 

Jo in * X 1 


'0 "0 


- J 2 | ^(O.y.tjj-x) £ p(w(y,x) ) g(v(y ,x ) ) - f (u(y,x) ) g(v(y,x) ) 


dy dx 

dy dx 
(3.4.8) 


Since, |u(x, t) - w(x,t)| = |v(x,t) - v(x,t)i and at x * 0, u(0,t) « w(0,t) 
6(t). We have for all t, |v(0,t) - v(0,t)| « 0. Therefore, for g «= c“. 


and also 


pCeft^)) j^gtvto.tj)) -gtvto,^)) 


p(e(t 2 )) [ g(v(o,t 2 )) - g (v(o,t 2 )) J = o 


(3.4.9) 


(3.4.10) 


Therefore, (3.4.8) simplifies to 

( [ f(e(t i )) ~ P (0(t i )} ] ■ [ f (ecy) - p(e(t 2 )) J | g(v(o,t j )) 
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f(e(t )) - p(e(t )) 
2 2 


= J g(v(0,t 2 )) - gCvfO.tj)) J 

, f t f 00 r 

+ [ J J K xx ( 0 ,y,t r T) [P (w{ y> T)) g(v(y,x)) - f (u(y,x) ) g(v(y,x))J dy dx 

- J 2 J K xx (0 -y-V T) [pfwCy.x)) g(v(y,x) ) - f(u(y,x)) g(v(y,x))J dy dx j 


= I + I 

14 IS 


( 3 . 4 . 11 ) 


STEP 6 : DERIVATION OF ESTIMATE FOR I + I . 

14 15 


Since > t 2 , I is expressed equivalently as 


- J 2 J ^(O.y.^-t) - K xx ( 0 ,y,t 2 ~x) £ p(w(y,x)) g(v(y,x)) 


15 


o '0 


- f (u(y,x) ) g(v(y,x) ) J dy dx 

t r oo r 

J 1 J K xx ( 0 ,y, t r x) j| p(w(y,x) ) g(v(y,x)) - f(u(y,x)) g(v(y,x)) J dy dx 


t o 
2 


= I + I 

15,1 15,2 


( 3 . 4 . 12 ) 


I and I are estimated following earlier arguments. The very structure 

15,1 15,2 


Of I 

15,2 

that "t " 

i 

gives rise to three expressions analogous to 

( i = 1, 2) stands in place of "0" and “ t“ in 

1 , 1 and I except 

2 3 4 

the integral. 

I 

15,2 

f 1 f K (0,y, t -x) 

J t J o xx 1 

2 

£ p(w(y,x)) - f (w(y,x)) J 

g(v(y,x) ) dy dx 

+ 

f 1 f K (0,y,t -x) 

J t J o ** 1 

2 

£ f (w(y,x) ) - f (u(y,x)) 

g(v(y,x) ) dy dx 

+ 

pt pCO 

J 1 J K (0,y,t -x) 
J t J Q xx 1 

£ g(v(y,x)) - g(v(y,x) ) 

f(u(y,x)) dy dr 


i + i + i 

15,2,1 15,2,2 15,2,3 


( 3 . 4 . 13 ) 
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Let f, g € C°\ for 0 < « < 1- Then, 


,t J» 


T = f 1 f K (O.y.t -t) 

15,2,1 J t J 0 *x 5 


♦ I s f ,t „ (0 ' y -V Tl 

J t J o 
2 


p(w(y,T) ) - f (w(y,t ) 3 

- p(w(0,x 3 ) ♦ HvIO.t)) 
p(w(0, x 3 3 - fltf(O.t)) 

g(v(y , t ) ) - g( v(0 , t 3 3 


F. f v ! y , t 3 3 dy dx 


dy dT 


Therefore, 


5 , 2 , 1 ’ £ { 


U 15 , J £ i lgl m ip - n w ♦ lgi w Hv^lJ Ip - f I, 


C (2» 3 

-i— it, - 1 r 

CT 1 2 


*{ >«'. 


Clw | 3 + igl Civ | ) 

XX 00 OC XX 00 


% ~1<X+Z) 3 


CJ2<*} 


Ip - f|„ It, - t 2 l 


< A 


8 It - t | 

1 2 


(3.4.14) 


Let f e Lip class and satisfy the property S. Then I is estimated 

i n m* f 4m 

analogous to I 3 . The estimate can be seen to be 


II , ,1 s A, |u - w | |t - t | 

15,2,2 2 xx xx co 1 2 


{3.4. 15) 


Similarly for g e Lip class and satisfying the property S, I has the 

estimate 


IL, J 1 A |w - u 1 It - t j 

15,2,3 3 xx xx co 1 2 


{3.4.16) 


Therefore from (3.4.14), (3.4.15) and (3.4.16) we get 

l I ,»l 5 MP - f L It, - t 2 l“ * A, - u I It - t, (3.4.17) 
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Before estimating I , we note that I can be rewritten as 

I i 15,1 


t CO r- 

[ = f 2 f i 

15 > 1 J 0 J o L 


K (0,y, t -x) - K (0, 

XX 1 XX 


.y,t 2 -x)j p(w(y,r)) g(v(y,x) ) 

- f (u(y,x) ) g(v(y,x) ) j dy dx 

r t jx> r t r 

= j 2 j J 1 K xxt (°*y> 1_ x) [ p(w(y,x) )g(v(y,x) ) - f (u(y,x) )g(v(y,x)) 1 dy dx 
° ° 1 2 

t jx> f t j- - 

= J 2 j j 1 K xxt (0 ’ y,1_T) [ P (w(y > x)) - f (w(y, x) ) g(v(y,x)) dl dy dx 

° ° t 2 

+ J 2 J J 1 K^JO.y.l-x) [ f(w(y,x)) - f (u(y,x) ) J g(v(y,x)) dl dy dx 


f t -CO t r -j 

J 2 J J 1 K xxt (0 ’ y,1 ~ T) [ g(v(y,x)) - g(v(y,x) ) J f(u(y,x)) dl dy dx 


= I + I + I 

15,1,1 15,1,2 15,1,3 


(3.4.18) 


Below we estimate these three terms separately. Making use of the identity 
property and noting that f, p and g are in C** for 0 < « < 1, we estimate 
I as follows. 

15,1,1 

J i 5 i,i = J 2 J J 1 K xxt (0,y ’ 1-t) [ P (w(y > T ^ " f(w(y,x)) 

0 0 t 2 

- p(w(0,x) ) + f(w(0,x)) J g(v(y,x) ) dl dy dx 

-t J00 t r -| 

+ J 2 J J 1 K^CO.y, 1-x) f (w(0,x) ) - p(w(0,x)) x 

° 0 t 2 

g(v(y,x)) - g(v(0,x)) j dl dy dx 


II 


o o t 


2 

| lgl w dl dy dx 

+ f 2 f f 1 |K (0,y,l-x)| |p - f| lg(v(y,x)) - g(v(0,x))| dl dy dx 

J 0 J 0 J t 

2 
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Hence , 


II, _ .1 * IP - f l« lgl » £lw *x ! » ) 

15,1,1 


t r t „ f “ 1 ^ 

IK (n.y.l-T)l | 


XX t 


r r j 

J t J o ~ 

2 

* .p - fi. m?jj* i.i. f; CC ^ 10 '— 1 " 

2 

Which by Appendix A (A8), gives 




dy dt dl 


dy dt dl 


C (2«) a 

u i — ip ' f *<* (|w J® 5 lg '« U i 

15 > 1 ' 1 «(!-«) 


+ 2 


C (2<x) 


t 


<x (1 - «) 

s A 2« Ip “ f l. !t , - S' 


(3.4.19) 


Where 


A =2 
26 


-ot 4 


C (2«) 


«(!-«) 


( I U JJ“ >«*» * lgl . 


For f € Lip class satisfying the property S. Following the steps of 


I , I . can be estimated as below ; 
3 15,1,2 


U 15jl2 l £ | 2 J j 1 IK xxt (0.y , 1 -t ) | |f(w(y,T)) - f (u(y.T)) | 

° o t 2 

lg{v(y,T) ) i dl dy dt 

3 c . Iu » - “J. i«i. J* a f j ' 1 " : xxl < 0 -y- ! - r >i ■? dl dy dT 

2 

S C S lg L C ,* 2i l U - w | It - t | I log 1 1 - t I I 

s oq 4 XX XX CD 1 2 I * * 1 2 * 


s A„ |W - U | It - t | 

27 xx xx a> 1 i 2 


(3.4.20) 
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Where 


A = C lgl C (2) [log |t - t I | 

For g e Lip class and satisfying the property S, on similar line I. 
can be estimated as : 


15,1,3 


II I s D |v - V | |f | f 2 [ 

15,1,3 s XX XX CO J o J o J 


I K xxt C 0 , y , 1 — t ) | y dl dy dx 


S D . C 4 (2) iU ,x " "J. If| «, U , ' V l 10 ^ I 4 , - S' I 

s A 2B lu « - U xx'» |l . - ‘a 1 (3.4.21 ) 

Since from (3.3.38) 

|u - w | £ A |p - f| t“ /2 (3.4.22) 

xx xx oo 24 H H oo 1 

Hence substituting (3.4.22) in (3.4.20), (3.4.21) and (3.4.17) one gets 
1 1 I - A |u - w | It - t 1 . 

15,1,2 27 xx xx oo 1 2 

£ A A Jp - ffl t K/Z |t - t | 

27 24 " ~ Hex 1 1 2 


Where 


£ a lip - fj it - t i 

29 11 11 cx l 2 


A = A A t 0 " 72 

29 27 24 1 


(3.4.23) 


and similarly 


II I i A A 
15,1,3 28 24 


I t * 72 it - t i 

“ oc 1 1 2 


* So Ip - f US ' V 

l'u.,1 s ( A 1 + A , S, *r IS - S'” ) l p - f L 'S - S 1 ' 
s A. I p - f L 'S - S'“ 


( 3 . 4 . 24 ) 


( 3 . 4 . 25 ) 
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, , . . m <g) I'l A ?'i) and 13.4 24 1 wr rnn< luic 

Similarly collecting (3.4.19), 

S 'W .! 1 ' ^ 

|p-f|„ It, - t/ * |p- f|„ It, - t 2 l * V lr - >l„ It,- t, 

‘(V'* 8 'VV M *VW '1 |p - fI « ll . - V 


* 32 Ip - f L «\ - 


(3.4.26) 


In view of (3.4.25) and (3.4.26), the following estimate of l j& is clear : 


it I = f* 1 [ W K (0,y,t -t) |p(w(y,T))g(v(y,T)) - f (u(y, t ) )g(v{y.x) ) 

1 15 ' J 0 J o xx 11 


dy dx 




’ o * o 


- f 2 f K (0,y, t -t) p(w(y, x ) )g(v(y, t 3 ) - f (u(y, x 3 )g ( v (y , t ) ) 

L J n xx 2 L 


dy dt 


< |I | + |I | 
1 15, l 1 15,2* 


£ A + A 

32 31 


) b - f l„ 
* A * b - f L i‘, - V" 


it, - t s i 


(3.4.27) 


Now coming to the estimate of I , given by 

14 


I = f g(v(0,t )) - g(v(0, t ) ) 
14 L 2 i 


f (e(t )) - p( 0 { t ) 3 
2 2 


We have, for g e C 




n u \ s | g i a |v(o,y - vfo.yT |f(e(t 2 )) - p(e(t 2 ))| 

s '*'« |f - pi. i‘, - y“ 

* '*•. tivj.)- if - ?|„ it, - 1 r 


S A 34 l f - Pi. I 1 , - ‘,1' 


(3.4.28) 
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STEP 7 : ESTIMATION OF «-SEMINORM OF f - p : 



fCeft,)) - pce( t i ) )J - f(e(t 2 )) - pce(t 2 )) ] J |g(v(o, ) ) I 


ft p® r - 

s | K xx^°’ y ’ t i _T) Pfwfy.T) )g(v(y,T) ) - f (u(y,x) )g(v(y,x)) dy dx 

J o J o xx J 


^ 00 r 

- f 2 f K xx (0 ' y,t 2 _T) p(w(y,x) )g(v(y,x) ) - f CuCy.x) )g(v(y,x) ) dy dx 
•’ o J o xx *- 


+ [g(v(0,t 2 )J - g(v(0,t t ))J f(e(t 2 )) - p(0(t 2 )) 


| | |^f (e(t i ) ) - p(e(t i )) J - £ f(o(t 2 )) - p{0(t 2 )) 

: A Up - f I It - t r + A If - pi It - t r 

33 11 « 1 2 34 11 K|l oc 1 2 



Igl 


(3.4.29) 


Now dividing both the sides of inequality (3.4.29) by |0(t j ) - © C | , 
the «-semi norm of | f — pi will be obtained 


If - pi 


sup 

0(t )*0(t ) 

1 2 


f(0(t i )) - pfettj)) | - J f(0(t 2 )) - p(0(t 2 )) 

|e(t i ) - 0 (t 2 )T 


S | inf |g| j 


-1 


sup 

0 ( t )*0(t ) 
1 2 


It - t I 
1 2 


33 i0(t ) - o(t )T 

1 4S 


»p - fl 


+ A 


•S - ‘a 1 


34 


|0(t ) - act ) I 

1 2 


if - pi. 


Hence , 


if - pI„ = [ inf Igl ) ’ | |p - f| K + a 31 |f - p|„ } ( f»f I®' I ) 


s A 3s Ip - f L + A „ l f 


(3.4.30) 


CX 36 
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We recall the expression (3.3.39) and note that 


and 


• f " Pi„ * A 25 l P " f L 


!f - * A 3S l P ' f L 4 A 3 . I f ~ Pi, 


Now adding up the two estimates above, we obtain f-nmm. of Jf - £| 


Hence , 


If - p| a = if - * if - ?i, 


3 A 2s I* 1 - f l„ t“ * |p - f |, • a jo | f - 


( A * *“ + A 3s ) I p - f h * A 3 * K 


pi, 


Therefore, 

l f - pL : 


f At” + A 1 

25 35 


1 - A 


Ip - f|. 


36 > 


s A 37 Ip - f l« 


£3.4.31) 


Thus we have got an inequality where <*~norm of (f - p) is bounded by the 
« norm of (f p). Therefore, we reach a contradiction to the assumptions of 
the existence of two distinct solutions provided 0 5 A < 1 . 


37 


STEP 8 : ESTIMATE OF A : 

37 

Define A = A t” + A + A 

38 25 35 3^ 

If we said that < 1 then t“ * A_ < 1 - A 


35 


36 


provided A < 1* 


36 


So 


A t + A 

A = L 5 35 

37 


1 ~ A 


36 


is equivalent to A < 1 . 

38 


For convenience A . i = t -5-7 

1 .... ,37 are rewritten to estimate A 

38 
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A = [igl ( |w i )“ + |g| C j v | )“ ] 

1 L XX 00 CX XX 00 J 


C(cx) 


‘(1 - oc) 


A = —j— C Igl C (2) 

2 4 S CO 3 


A 3 - T D , If l« C 3 (2) 


A = A + A 

4 2 3 


1 

4 


c 3 (2) f C s l«L + D if! ] 
3 L 8 00 S 00 j 


C (oc) 

A = Igl ( iw | )* — 

5 co x oo c < 


C (oc) 

A = Igl ( I v | )“ 3 


« X oo oc 


C (oc) 

A = A + A = — 

7 5 6 oc 


r igi dw i )“ + i g i civ i ) w i 

|_ 00 X oo oc X oo J 


A = C |g| C (1) 

8 s oo 3 

AO' 0 , 

A l2 = A 8 + A 10 = C 3 (1) [ C . I*'. +D , l f L] 

A 13- A , * A „ =CDC , (1) [ + 

A ,S * C I*'. C 2 ' 0) 

A it = A .s + \e - [ C 1*1. * D l f '„ ] c 3 (0 > 

A » * C I*'. S CO) 


A = C D (|v | ) C (1) 

9 x oo 3 

A = C D (|u | ) C Cl) 

11 x 00 3 


= C 2 ' 0) '*'. 
\e - D lfl . c 2 ‘°> 

A !8 = 'S'. C l (0) 
A 30 - D l f '. C ,' 0) 
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A 2, ' * A » ‘ 


C Igl • D If I, 


r {pi 

5 


18 

A « 

22 1 - A t 

21 


I " 


k t: i p ) 

1 r ' m J 


c jg ♦ r i i 

1 r 85 I 


r ten t 

i 


A = A + A,„ A t 
23 H 17 22 


- c 2 (« Igl, * t 


C Igl, * D If I, 


R f 10) 

® i 


C (01 
2 


C * n If! 


. » ^ « » , 1- (ot/2) , . . . (3-»)/2 

A = A + A A t + A A t 

24 7 12 23 S3 22 


C («) 

3 


r 

— Igl 

OC CO 


<1« IJ + Igl (Ivl ) 

X CD X CO 


C 3 (l) [c. |g|_ + D_ If !_ | X 


£ CO s 00 


C (0) Igl + t [ C Igl + D If | 

45 <0 CO 60 


cyoi 


♦CD C (1)[ t|? | ) + (|u | ) 

3 X oo X CO 


, (3-OCJ/2 


1 * 1 . V 01 

c lei. * d |f|. 


lei c io) 

° 00 1 


C |g| * D If I 

° *00 m 


f Inf |g| 1 1 

[ A + A A 1 

l J 

L 1 4 24 J 


- ( Inf Igl ] | [l8l„ + Igl. (IvJJ- ] 

* T V 2 ’ [ C. Igl. ♦ D, lf|. ] t- IW » } 


C(«) 


<x(l - «) 


c {05 t 


C j £0} t 


(0) t 
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C(cx) 


A 26 = 


oe( 1 - o< ) 


igl ( lw | ) + |g| ( | V 

CO XX CO OC xx 




A = C Igl C (2) (log It, - t. 

27 s <»4 1 ; 


A = D | f I Cl 2) (log |t, - t. 
28 s oo 4 1 1 2 


A = A A t 

29 27 24 1 


cx/2 


C lg L C 4 (2) I lo 8 It - t || t “ /2 A 

s 00 4 1 1 2 1 1 24 


A = h A t,^ 2 = D If | C (2) I log It - t | I t ° C/2 A 

30 28 24 1 s co 4 • 1 2*1 ; 


24 


A = A + A A t** 72 It - t I 1 “ 

31 1 4 24 1 1 2 


C(oc) 


cx(l - oc) 




(lw 1 ) + |g| (|v | 

XX 00 oc XX 




+ -t C 3 (2) C Igl + D Ifi 

3 E oo s co 


t^ 2 It - t A 

1 1 2 24 


A = A + A |t - t 1 1_OC + A |t-t| 1_0C 

32 26 29 1 2 30 1 2 1 


C(oc) 


«(1 - oc) 


[ 


Igl (lw I ) + Igl (|v I ) 

co XX co oc XX co 


] 


* [ c , * d , i f '„] - vi c . 


(2) It - t | 1_<X A 

1 2 24 


A = A + A 

33 31 32 




A 3S - (inf Igl ) 


-1 


33 


(inf 1 6'\r 
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-( 


a 3s - | taf Igl) 


A 34 -1 Uv I )“ 

^ = (inf ig[) 1 |g| 1 - 


(inf |e' | ) 


(inf IS' | )‘ 


A = 

37 


A t + A 

25 35 


1 - A 


36 


Now one can estimate the the following 


A = A + A + A t 

38 36 35 25 


= (inf Igl ) -1 Igl 


civ | ) 


t 00 


“ (infle'j) 0 ' 




(inf |gj )’ * [ C(oc) 

2 


(inf |e' | ) c 


«(1 - oc) 


{ lgl « * "w V" Isl 


* T C 3 (2) [ C , 1*1. ♦ D |f|_ 


,oc /2 


t. It. - t 


, l-a 


I 2 24 


-C 4 (2) [ C> Igl. + D. IflJ | l0 g | ti 




2 24 


(inf |g|) 


-1 , oc 


C(cc) 


a \ r* 

t — |g| (jv j ) w + (J W | )« . . 

«(1 - oc) L “ X* » xx» igf , 


+ ~T C 3 (2) (inf fgD" 1 [ C |g| + D J f I ] 

L B » E CO J 


K t 

24 


l+(«/ 2 ) 


* Unf Igl )" 1 |g| 

« . . 


(Jvj )" 


(inf Je' | )** 


+ ~r c 3 (2) 


[ c . + in.] 


tr it, 


(inf |gj ) -1 

t | A 

(infle'D* 24 
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+ 


(inf lg( ) -1 

2 + (inf Igl)' 1 1“ 

(Inf ie' | ) 


x 


C(oc) 


ot(l - ot) 



( |v I )' 

XX 00 


+ ( |w 

XX 




(inf |g| ) -1 


(inf |0' I ) 


— C 4 (2) | log It, - t 2 l| It, - t 2 \“ 


* (inf Igi)" 1 t 1+(oc/2) -J-C 3 (2) 


C Igl + D If | ] A 

s oo s oo 24 


Let coefficient of A be 5 (t). Now calculate S (t) A 

24 1 1 24 


(3.4.32) 


6 (t) A 

1 24 


C («) 
set) — 

1 <x 


r |g L ( ,w > )<x + igi uv i )“ 

l CO X 00 <X x CO 


! . (t) [ 


+ 3 (t) 1C Igl + D |f | 1 C (1) t 

E co s oo J 3 


l-(oc/2) 


c (0) Jg| + t f C Igl + D |f | ] C CO) 

<2 oo I co oo J 2 


gl C (0) 

CO I 




+ S ( t ) | C I v | ) + (|u | ) CD C (1) t 

4 1 x oo X 00 3 


(3-t>0/2 


»- [c 1*1. * D l f l»] S' 01 * J 

1*1. C J (0) 


* - [ c 1*1. 4 D l f l. ] c , (0) 1 


(3.4.33) 


Represent the second term of R. H. S. Of (3.4.33) as 5 2 (t). ^ ow 
substitute the above expression in (3.4.32) and with some suitable arrangement 
this can be written as 
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A = (inf Igl)' 1 igl 


38 


(lv I ) 

L_? + 5 (t) 

(inf ie' I ) a 2 



- 

C(«) 

lv | ) 

+ (inf Igl r 1 

1 , I * ©< 

2 + (inf j8' | ) t 

Igl 

00 

XX oo 


_ 

oc( 1 - <x ) 

1 inf 10' | y 


+ (inf Igl) 


2 + (inf 1 8' j )** t C 


C{«) t 

Igl 


c(l - oc) 


XX 00 


( inf i8* i j 


|w | ] 

(X 

lv 1 1 

ex 

lu | ) 

X co 

+ 8 (t) 

4 

X OO 

+ 8 (t) 

5 

X CO 

k inf 18' | 

k inf 10' 1 ; 

k inf 10' | ; 






Where 5 (t), 8 (t) and 5 (t) are the coefficients obtained with some suitable 

3 4 5 


arrangements. 

So it is enough to show. 


Ivi ] 

<X 

|w | 1 

OC 

r iv i i 

t oo 


xx oo 


XX co 

k inf. 18' | , 

* 

k inf 10' | i 


k inf 10' i , 


( lw I ( lv I 1* 


X » 


inf |8' | 


X CD 


{ inf 18' | J 


and 


< |U J» 1 


inf |8' 


are bounded. Utilizing the 


representation of the previous chapter 


A *-» -w 

W (x.t) = 0 (x,t) + j K (x,y, t-x) p(w(y,r)) g(v(y,x)) dy dr 

XX XX J o J o xx 

= *xx (x ’ U + J J ^(x.y.t-r) |g(v(y,T)) jp(w(y,x)) - p(w(x,x))j 

+ p (w(x,t)) | g(v(y,x) ) - g(v(x,x) ) j- j dy dx 

* '*= - *J * [ 1*1. Ipl. ♦ Ipi. Isl„ (|5 X IJ‘ 


T C («) 


!_ t^ 2 


(3.4.34) 
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t 00 

Iw x - s J J |K x (x,y, t-x) | lp(w(y,T))| |g(v(y,T))| dy dr 
* o J 0 


s c to) vr ipi igt 

2 co co 


(3.4.35) 


Similarly, v can be represented as 

t 00 r 

v (x , t ) = 0(x,t) +| | K(x,y,t-z) - p(w(y,r) ) g(v(y,x)) J dy dr 


t. 00 

" - IPi m f [ IK (x,y, t-x) | dy dx 

» A * A 


" fP L lgl m C ? {0) ^ 

oo co 2 


(3.4.36) 


and also 


v = <p + f f K (x,y,t-x) - p(w(y, x) ) g(v(y,x)) 

XX XX J J XX [ 


+ p (w (x , x ) ) g(v(y,x) ) dy dx 


|V -Ms Igl Ipl Clw I ) + Ipl Igl (|v | ) X 

XX XX OB OC X00 CD « X 00 


* f Igl Ipl ( |w | )* + Ipl Igl (|v I )“ 

oo oc x oo oo oc x oo 


t 00 

f f |K (x,y, t-x) | |y - xl* dy dx 
•* o o xx 

1 C (<x) 

IpI «d lgl « nvM- ] — — t“ /2 

J oc 


00 oc x oo 


t*' 2 (3.4.37) 


Now substituting the bound for |w | and |v I, (3.4.34) and (3.4.37) will 

XX 

become respectively 

iu ,« - 5 [ if. if. { i*„i + c 2 ( °> ^ ipi. if. f 


CC C (oc) 


( } 1 3 l ; 

+ |p L i«L i i* i + IpL i«L c (o) vr v t 

co «l .x co oo2 J J 


(3.4.38) 


and 
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i; -♦„! * f 1*1. ipi. { 'V. * V 0) ^ lpl » 

XX xx l v 


+ IpL ! «L 


( I* I 

1 X CO 


+ c (o) VT” Ip! Ig 


4 


a C («) 

_J t“ /2 

cx 


(3.4.39) 


From these two inequalities (3.4.38) and (3.4.39), we can conclude that 


w - Ip 

XX XX 


and v - (p 

XX XX 


0, like 0(t* /2 ) 


if ML ~ E and l*L ~ E 


(3.4.40) 


The last part (|v | )** which is left, has been calculated in chapter II 

r t CO 


Now we define the sequence of iterates 


/ (n) l 

t p j 


defined by 


P (n+1) (0(t) = T Ip'"'] (t) 


(n) - 


since f is a fixed point of T. , clearly by the estimate 

U 


|f - p <n * 11 L * A_ If - p'"’’ 


oc 37 


and consequently p (n) - f — » 0 as n — > to. Since C** is complete, so p (n) — > 

f in c“ To see that T has a unique fixed point, we proceed by the method of 

contradiction. Suppose h is another fixed point of T , but we know that 

0 


l h - f|« - I I e lW -T e [f] | k = A 


37 


f fi« < I h - 


This is possible if f = h. 

This completes basic assertions of this chapter. 

3.5 : CONCLUSIONS : 


In this chapter we have shown the uniqueness of a solution which is 
established by the method of contraction* A large number of estimates have 
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been divided to arrive at the results. The constants are derived explicitly. 
With the assumptions on flat initial values the local contractively result is 
true for small time interval. The present work is possibly the first work for 
a coupled system of two equations along these lines. 

There are several directions in which the work remains to supplemented. 
First of all we have assumed the Neumann boundary condition as the identity 
property is not satisfied for Dirichlet boundary conditions. So this theory 
will not be applicable to Dirichlet Boundary Value problem. Hence, possibly a 
different approach needs to be looked into. Recently HOLLIS and MORGAN [60], 
[62] have studied the wellposedness of the direct problem for system of R-D 
equations (n 2) considering mixed boundary conditions on a global boundary. 
Therefore, the inverse problems for more than two equations with mixed 
boundary condition on local boundary (90 = 91^ + dn^) data possibly be studied 


later. 




CHAPTER IV 


EXISTENCE OF A SOLUTION OF THE SOURCE PROBLEM 
IN A HEAT CONDUCTION EQUATION 


4.1 INTRODUCTION : 


In this chapter we have considered the existence of an unknown forcing 

function as a function not only of the unknown solution but also of the 

spatial variable x. Assuming variation to be space variable, the form of the 

function is assumed to be f (u) + x f (u). In view of the two independent 

1 2 

functions entering into the unknown right hand side, boundary conditions are 
assumed at both the boundary points. Earlier PILANT and RUNDELL (1986,1988) 
[102], [104] have studied the problem in a semi-infinite domain considering 

right hand side as a function of "u" only. 

The preview of this chapter is as follows. In section 2 the inverse 
problem is formulated. In section 3, an iteration scheme is outlined. An 
equivalence relation is established in section 4. Section 5 deals with the 
existence of a fixed point of a nonlinear map and a brief conclusion is given 


in section 6. 



84 


4.2 FORMULATION OF THE INVERSE PROBLEM : 

We consider the following Heat equation in a spatially bounded domain, 

u (x.t) - u (x,t) = (1 - x) f (uCx.t)) + x f (u(x,t)) + yCx.t) ; 

0 < X < 1 , t > 0 ; (4.2.1) 

u (0,t) = g (t) ; t > 0 ; 
x 0 

u (l.t) = g (t) ; t > 0 ; (4.2.2) 

X 1 

u(x,0) = u q (x) ; 0 — x — 1 ; (4.2.3) 

Where, y is a known function and the functions f^.) and f 2 (.) are unknown 
and are to be determined with the help of overposed boundary values of u at x 
= 0 and x = 1, given in the following form. 

u(0, t) = 0(t) ; t > 0 ; (4.2.4) 

u(l,t) = q(t) ; t > 0 ; (4.2.5) 

An introduction of a comparison function i />, defined below will reduce the 
complexity of computations. 

(x.t) - $ (x.t) = y(x,t) ; 0<x<l , t>0; (4.2.6) 

t XX 

O.t) = g Q (t) ; t > 0 ; 

0 x (l,t) = g x (t) ; t > 0 ; (4.2.7) 

tfr(x.O) = U q (x) ; 0 £ x £ 1 ; (4.2.8) 

If one substitutes the values of u at x = 0 and x = 1 in the equation (4.2.1), 
then the following two equations result : 

0,(t) ' U xx (0,t;f i ,f 2 ) = V e(t)) + rCO.t) (4.2.9) 

q,(t) ' U xx (1 ’ t;f i >f 2 ) = f 2 (c I (t)) + yCl.t) (4.2.10) 

The above two equations provide the motivation for an iteration scheme to 
determine the functions f and f . 
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4.3 ITERATION SCHEME 


In order to find the solution of the inverse problem (4. 2. i )-(4. 2. 5) the 
following iteration scheme is developed. For any given initial choice 


(f i ),f 2 } ° f (f r f 2 ) the corresponding direct problem is solved. Let 


u(x,t;f i <0) ,f^ 0) ) be the solution of 


ujx.t) - u xx (x,t;f i <0) ,f^ 0) ) = (1 - x) f x (0) (u(x ( t)) 

+ x f^ 0) (u(x,t)) + r(x,t) ; x e (0,1) , t > 0 ; 
V X ’ t} = S 0 (t) * u x (1 ’ t) = «! (t) J t > 0 ; (4.3.1) 

u(x, 0) = u (x) ; X € [0,1] ; 

However, in general u(x, t ; f^ 05 ,f^ 0) ) will not satisfy the Dirichlet Boundary 
data at x = 0 and x = 1, i.e, 

u(0,t;f i (O) ,f^ O) ) * 0(t) and u(l, tjf^’.f * 0) ) * q(t) 

Therefore, one needs a strategy to update (f (0) ,f (0) ). Let (f (1) ,f (1> ) be 

12 12 

the new value generated by via the following updating scheme. 


f ci> (e(t)) = 6' (t) - u (0,t;f (O) ,ff > ) - r(0,t) 

1 xx 1 2 


(4.3.2) 


and 


^“(qft)) = q'(t) - u^d.tif^’.f^) - y(l,t) (4.3.3) 

Proceeding along these lines, one has the following for k e N. 

f (k+1> (e(t)) = e'(t) - U (o,t ; f (k) ,f‘ k) ) - r(o,t) (4.3.4) 

1 XX 1 2 


f‘ k+1, (q(t)) = q'(t) - u (l,t;f (k> ,f* k> ) - y(l,t) 

2 XX X 2 


(4.3.5) 


Before proceeding further, we introduce the following notion : 



t [f ,f Ht) = [ e'(t) - u - y(o,t) , 

0»q 1 2 L 

q' ( t ) - u (l,t;f .f ) - 7(1, t) (4.3.6) 

XX I Z 

Therefore, the iteration scheme can be restated in terms of T Q q operator as 


>.,[ c’-r] 


= [ e'(t) - u (o,t;f (k> ,f‘ k) ) - r(0,t), 

XX 12 

q'(t) - u (l,t;f‘ k, ,f‘ kl ) ~ y(l.t) 

XX 1 ^ 


f (k+1) (e(t)) , fr i; (q(t)) - T 


(4.3.7) 


which on adding and subtracting f ™ ( . ) for j = 1, 2 and replacing y - u^ by 
its expression in the equation (4.3.1) becomes 


= f‘ k) (e(t)) + | f‘ k) (u ( 


(0, t;f i ,f 2 )) - r (©{t ) ) 


00 fA . -( k ) ,.00 •, 

+ 6' (t) - U (0 , t ; f ,f ) •, 

t 1 * J 


f' k> (q(t)) + | - f^ (q(t) ) 

f t \ \ (k) M 4. r 00 r CkK \ ] 

+ q' (t) - u t (1, t;f i ,f 2 ) M 
J fj k) (e(t)) + Fjj e(t) - u (k> (o,t) j- , 


r (q(t)) + Fj q(t) 


u«Cl.t)}] 


Here F and F are in general nonlinear functions. 


4 . 4 . EQUIVALENCE RELATION : 


The following lemma shows the equivalence of the solvability of the 


inverse problem to the existence of a solution of a fixed point problem. 
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DEFINITION : A function will be called (e.q)-fixed point of T if 

0,q 

given 0 = e ( t ) and q = qCt), 


f (0( t) ) , f (q(t>) 
1 2 



f 2 l(t) 


(4.4.1) 


If f t and f 2 are uniformly Lipschitz, then the existence of (0,q)-fixed 
point of the operator T is equivalent to the solution of the inverse 

H 

problem of determining the source function with the help of the superposed 
Dirichlet data. We state it as a lemma. 


LEMMA 1 : If f and f are Lipschitz function, then (u, (f ,f )> is a 
12 12 

solution of (4.2. 1) — (4.2.5) if and only if (f^.f ) is a (0,q)-fixed point of 


PROOF : Suppose (u, (f ,f )} is a solution of (4. 2. 1 ) — (4. 2. 5) . Now we show 

12 

that (f ,f ) is a (0,q)-fixed point of T . Since (u,f ,f ) is a solution 

1 2 D,q 1 2 

triplet, it satisfies the following relation : 

[ f a (0(t)) , f 2 (q(t)) j = £ f t (u(0,t)) , f 2 (u(l,t)) J 

- [ u (0, t) - u (0, t;f , f ) - r(0,t). 

It xx 12 

u (l,t) - u (1 , t ;f , f ) - y(l,t) 1 

t xx 1 2 J 

« [ e'(t) - u (0, t; f ,f ) - 7 ( 0 , t), 

L XX 1 2 

q'(t) - u (l»t;f ,f ) - 7 ( 1 , t) 1 

XX 1 2 J 

= t r f ,f ] (t) 

e.q [ 1 2 J 


So (f ,f ) is a (0,q)-fixed point of T. . 
12 q 
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Conversely , assuming (f,,f 2 > 1= a (e.q)-flxed point of T^. we claim that 

lu,(f ,f )1 is a solution of the inverse source problem. It is sufficient to 

show u(0, t) = e(t) and u(l.t) = q(t) for the I. B. V. problem (4.2. 1 >-(4.2.3) 

with f and f as the source functions. 

1 2 

Since {u,(f ,f )> satisfies the differential equation (4. 2. 1 )- (4. 2. 3) we 
1 2 


have 


at x 


y(0,t) + U^tO.tjf^f.g) = ~ 

0 , and 

y(l,t) + u (l,t;f ,f ) = u (1 . t) - f (u(l,t)) 

xx 1 2 t ^ 


(4.4.2) 


(4.4.3) 


at x = 1. 

If (f ,f ) is (0,q)-f ixed point of , then from the equations (4,4,2) and 

(4,4,3) we have, 

( ye(t» , f 2 (q(t» ) = [yrj (t) 

= J e' (t) - | u^to.tjf^f.,) + r(o,t) j- , 
q' (t) - | u^d.t;^,^) + r(l,t) | j 
= [ 0'(t) - u (0,t;f ,f ) + f (u(0,t)) , 

I t 12 1 

q'(t) - u (l,t;f , f ) + f (u(0,t)) 1 

*> 12 2 J 


Here y + u has been replaced by Its equivalent expression from the equation 

XX 

(4.3.1). Thus the following holds. 

j^e'(t) - u t (0,t;f it f 2 ) , q'(t) - u t (l,t;f lt f 2 ) J 

- ^(ett)) - f x (u(0,t)) . f 2 (q(t)) - f 2 (u(l,t)) 


( 4 . 4 . 4 ) 
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Denoting a and 8 by the expression below, 

oc(t) = 6(0 - u(0, t) 

and 

8(t) = q(t) - u(l,t) 

and noting that f jF f g are assumed to be Lipschitzian with the Lipschitz 
constants D and C respectively, we get 

la' (t) | = |9' (t) - u t (0,t) | 

= 1 f j (8(0) - f 1 (u(0 , t ) ) | (from(4. 4. 4) ) 

£ D |6(t) - u(0, t) | 

- D |a(t) I (4.4.5) 

and similarly, 

10' (t)| s C |0(O| (4.4.6) 

Since the overposed data 6(t) and q(t) must agree with the initial data u q (x) 
at x = 0 and x = 1 respectively, we have 

a(0) = 6(0) - u(0,0) = 0 and 0(0) = q(0) - u(l,0) = 0 

Then, by Gronwall’s lemma a(t) = 0(t) = 0 for t > 0. i.e, the overposed 
boundary values are also satisfied. | 

So the existence of a (0,q)-fixed point of T , ensures that there is a 
solution for the inverse problem under consideration. 

Thus the original problem reduces in establishing the existence of a 
(6,q) -fixed point of the T operator. 

6,q 
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4.5 EXISTENCE OF A FIXED POINT OF T e : 

First of all we turn our attention to the following preliminary 
considerations. Let K be the Neumann function for the bounded region (0 s x s 
1) with homogeneous Neumann data. Then K has the following form 


K = K (x,y,t) 


00 

J 

1 

r 2 2n 

- ( x - y- 2n ) /4 1 ~(x+y~2n) /4t 

e + e 

L 

-CO 

y^hrt 

- 


(4.5.1) 


If v denotes the difference between u and ^ given in (4.2.1) and (4.2.6) 
respectively then, v satisfies the following equation : 


v fc (x,t) - v xx (*.t) = (1 - x) f ^ (u) + x f ^ (u) , x e (0,1), t > 0 ; (4.5.2a) 

v (0,t) = 0 , v (l.t) = 0 ; t > 0 ; (4.5.2b) 

XX 

v(x,0) = 0 ; x > 0 ; (4.5.2c) 

Therefore, the solution of (4.5.2) can be written as 

v = u- ^ = J J K(x,y, t-x) ^ (1 - y) f^ufy.x)) + y f 2 (u(y,x)) J dy dx 

(4.5.3) 

Now differentiating v twice with respect to x, one gets 

V X* * u „ - *„ = f f f (uty.T)) dy dx 

J 0 J 0 

ft (.1 r 

+ J 0 J 0 Kxx(X,y,t “ T) y [ f 2 (u( y* T)) - f x (u(y,x)) J dy dx (4.5.4) 

Hence the second derivative of u at x = 0 and x = 1 are given respectively by 

u xx t0,t;f i ,f 2 ) * ^xx (0,t) + f f K (O.y.t-*) f (u(y,x)) dy dx 

J 0 J 0 1 

+ -L K « (0,y,t_x) y [ f 2 (u( y> T >) - fjCuCy.x)) J dy dx (4.5.5) 
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n.tsW s ^X (1 * U + J o J o K xx £l,y,t-T) f^uCy.T)) dy dr 
+ J J ^ [ fjfufy.t)) - fjfufy.x)) j dy dx 


(4.5.6) 


From the definition of T , we have, 

y, q 


WW (t) - [ e ' (t) - u xx (0 ' t:f 1 ' f 2 ) -rlo.t). 


q'Ct) - u a (l,t;f it f 2 ) - 7 ( 1 , t) 


(4.5.7) 


Now replacing r by ip - $ and u - t p by their expressions in the 

v XX XX XX 

equations (4.5.5) and (4.5.6), the equation (4.5.7), reduces to the form 
T [f ,fj (t) = f e'(t) - u (0, t;f , f ) - y(0, t), 

o,q 12 ^ xx 12 

q' (t) - U^d.tjfj,^) - y(l, t) J 

= j 0'(t) - u (0, t ; f , f ) + >p CO, t) - <M0,t), 

[ xx 1 2 XX t 

q'(t) - u (1, tjf ,f ) + * (l,t) - (l.t) 1 

xx 1 2 xx t J 

r r l r 1 

- e' (t) - ip (0,t) - J j K (0,y,t-x) f^ufy.x)) dy dx 

- f f K (0,y, t-x) y f f (u(y,x)) - f (u(y,x)) dy dx , 

J o J o xx L 2 J 

q' (t) - ip (l.t) - f f K (l,y,t-x) f (u(y,x)) dy dx 

*• J„ J„ xx 1 


'0 0 


- J J K xx (l,y,t-x) y £ f 2 (u(y,x)) - f^uty.x)) J dy dx J 

= [ 0'(t) - ^(O.t) - ^(f^), q'(t) - ip t ( l.t) - Wf-,) ] (4 * 


(4.5.8) 
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We need the following result for subsequent use and express it as a lemma. 

LEMMA 2 : The Neumann function K defined in (4.5.1) satisfies 
l r 1 

f K (x,y, t) dy = I K (x,y,t) dy — 0 for all x € [0,1] 
xx J o x 


PROOF : Indeed, 


CD - r 

K(x.y.t) = V 

n=-co /47lt *- 


2 2 
- (x-y-2n) / 4t , -(x+y-2n) /4t 

* + e 


so, differentiating K, with respect to x gives 

00 i r 

K (x,y,t) = 7 - 

x n=-00 L 


( X - y - 2n) -(x-y-2n)^/4t 

e 


2t 


(x + y - 2n) 


2t 


-(x+y-2n) / 4t 


Another differentiation of K with reference to x, yields 

X 


K (x,y,t) = y - 

xx n=-co VMi 


L-[- 


(x - y - 2n) ,2 

J - (x-y-2n) /4t 

e 


4t 


1 -(x-y-2n) /4t 

+ e 

2 t 

2 

(x + y - 2n) 2 

-(x+y-2n) /4t 

- e 


4t 


+ 1 - (x+y-2n) /4t 


2 t 



r , ,2 f 

(x - y - 2n) 2 V 

1 

1 - (x-y-2n) /4t f 1 

V4rct 

■ L irr 

4 t 2 J 


+ “ (x+y-2n) /4t f 1 

(x + y - 2n) 2 


^ 2 t 

4 t 2 


(4.5.9) 


(4.5.10) 


(4.5.11) 


In order to prove (4.5.9), we proceed step wise. Integration by parts 
plays a major role in establishing this equality. Inserting the expression in 
(4.5.11) for K , we get 
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f K (x,y, t) dy 

L xx 


** * 

CO 

f 

1" -Cx-y-2n) 2 /4t 

i 


(x - y - 2n) 2 ^ 

n= -oo V^Tlt 

J n 

L 

2 t 


4 t 2 J 


+ e 


■ (x+y-2n) / 4t 


or. 


1 2 t 


(x + y - 2n) 2 v , 

Tt 5 ) J dy 


1 CX) 


J K (x.y.M dy - f £ 


1 - (x-y-2n) /4t 

e 


2... (x - y - 2n) 2 


0 n = -oo V^nT 


4 t 


dy 


1 CD 


l J. 


1 - (x+y-2n ) 2 /4t ( x + y - 2n) 


o n= -to V4irt 


4 t‘ 


dy 


Now, 


r l oo 

r f 1 

r 2 

- (x-y-2n) /4t 

o J 

2 i 

- (x+y-2n) / 4t 1 

^ 0 n = -oo 2 t V^Tlt 


J 


dy 

(4.5.12) 


4 » , , ,2 ^ (x - y - 2n) 

1 - (x-y-2n) /4t , 

e dy 


r 1 00 

i l 


0 n=-co v^TTt 


4 t‘ 


4 » 


*4. W 

- J I — 

J o n=-oo VA reF 


, 2 (x - y - 2n) 

- (x-y-2n) /4t , 

e dy x 

2 t 


(x - y - 2n) 


2 t 


dy 


= J _1 ( lX " y 

n=-co Vint l 2 t 


(x - y - 2n) 2 -.1 

J g -(x-y-2n) /4tl 


* J* -J- 

J o 2 t 


- (x-y-2n) /4t , 

e dy 


} 


00 . r 

- I -4- f 

n=-co V47lt *■ 


(x - y - 2n) 2 

- (x-y-2n) /4t 

e 


2 t 


i: 


* i — r— 

n=-cD V4nt~ J o 2 t 


1 * e -(x-y-2n) 2 /4t dy 


(4.5.13) 
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Similarly, integrating the second term on the R. H. S of expressions (4.5.12) 
by parts we have 


09 

l 


(x + y - 2n) _( X+y - 2 n> 2 / 4 t 

— x e 


n=-co V4TTt 
1 


2 t 


t 1 


I 


n=~co 


v4rrt 


f 

J o 2 t 


e -(x+y-2n) /4t dy 


(4.5.14) 


Now substituting the respective terms of equation (4.5.12) by expressions 


(4.5.13) and (4.5.14) we get 


f K (x,y, t) dy = V — / 
J o xx n=-co vTO t 


(x - 1 - 2n) , 2 

- (x-l-2n) /4t 

e 


2 t 

(x + 1 - 2n) 


• (x+l-2n) /it 


2 t 


= 0 


(4.5.15) 


Similarly, integrating by parts one can show K(x,y,t) dy = 0 

* n 


and hence the lemma is proved. | 


COROLLARY 1 : If h is a continuous function, then 


f [ K (x,y,t-t) h(x,r) dy dr = f h(x,x) f K (x,y,t-r) 

J o J o xx Jq J o xx 


dy 


dr 


= 0 


(4.5.16) 


We record the following statement for later use in the form of a corollary : 
COROLLARY 2 : If for h € C* , « > 0, then 

f f K (x,y,t-r) h(x,r) dy dr 
J o J o ** 

~ J J ^(x.y.t-x) 1" h(y,r) - h(x,r) 1 dy dr (4.5.17) 
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In order to prove the main theorem of this section we record the 
following three assumptions : 

(Bl) u € C 2+W [0,1J, u (0, t ) = 0, u (l,t) = 0 

u x x 

(B2) The functions ©(t) and q(t) are monotone functions whose derivatives 



some 5 and 6 . 
i 2 


(B3) f i and f 2 e B E E j I . s E | for some constant E. 

Now we state the following basic theorem : 


THEOREM 1 : If the assumptions (B1)-(B3) are true, then for sufficiently 

small T and 0 < « < 1/2, T maps a ball in maps c“ x c“ into itself. 

tf , q 

PROOF : 


Given ^ and f g belonging to c“ , we define the functions f j and by 


so 


f(©(t)) , f (q(t) ) ] = T [f ,f ](t) 
L 1 2 I ©,q 1 2 


[ yeft)), ? 2 (q(t))J 

: [ 0' (t) - 0 t (O,t) - J J K xx ( 0 ,y,t-r) f^uty.x)) dy dr 
J J K.to.y.t-*) y [ f 2 (u(y,r)) - f^ufy.T)) j dy dr , 

q'ft) - 0 (l,t) - f [ K (l,y,t-T) f (u(y,x)> dy dr 

J o J o xx 

J J K xx (1,y,t-T ^ y [ f 2 (u(y,x)) “ fjMy.x)) J dy dx j 


(4.5.18) 


(4.5.19) 
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First of all we estimate the third and fourth terms in the expression for f 


Assume 


P t r 1 

I = K (0,y, t-x) f (u(y,x)) dy dx 

i J„ K xx 1 


o o 


f f K xx (0,y,t-x) I* f^ufy.x)) - f^ (u(0,x) } dy dx 
^ 0 o ** *• 


(at x = 0, the identity (4.5.17) gives the above result.) 


Using the fact that f e C K and expanding u(y,x) by Taylor’s expansion about 
(0,t), we obtain 


-t .1 

I I £ IK (0,y, t-x) | If | |u (0.x) y + u (0,x) (y z /2) 1“ dy dx 

ljjxx 1 « x XX 


By assumption Bl, u (0,t) is zero. So 


I I £ |f ,L {|u U* f f l K (0,y,t-T) | (y 2 /2)° < dy dx 
1 l«xxooJJxx 


0 "0 


£ 2" (<x+1) If I (|u | )“ -L C (0,0,2*) t* 

1 « xx 00 <x 3 


(4.5.20) 


(By (B3) of Appendix B). 


we can not apply the identity property for estimation of the last term of 
f (B(t)), because y is multiplied to K . However, it can be estimated 

1 XX 

directly as 



1 V s . 

J lK xx (0 ’ y ’ t-T) l ly l [ 

lf 2 '„ + |f ,L ] d y dT 


or, 

1 5 

. If . l »* l f 2 U C 3 (0.0,1) 

t U2 

(4.5.21) 


similarly, the estimate for is evaluated. Let us denote 
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.t 


1 « 
3 


J J o K xx (1,y,t ~ T) f j f uCy.x) ) dy dr 

J„ W-Y-W fjtuCy.T)) - f^ud.r)) J 


O '0 


u(y,x) can be expanded about y = 1 and we get 


1 V s J 0 I o IK x X n>y>t " T)! * u x (1 > t) (y - 1) 


+ ^ x (e ' T) ~ 1) Z /2)1°' dy dr 


XX 


Since by assumption Bl, u (1 ,t) = 0. Therefore we have, 


1 'a 1 s lf , L j' j’ IK„(J.y.t-r) 


1 y - il 


2<x 


s 2 


-(«+!) 


I f .l« C 3 (l,0.2«) 


The last term of f^ which we call 1^ has an estimate analogous 
i.e. , at x = 1 , this is given by 


III * f If I + If I ] 

4 ( 2 oo 1 co J 


+ ' f tL c 3 (1 ’°’ 1) 1 


1/2 


Hence, from the expressions (4.5.19) and the estimate (4.5.20), 
component of (4.5.19) can be estimated to get 


r .l- s I 0 ' - # t l. * 2 If,l« (l»„l„)“ C 3 (0,0,2«) f 


+ f If, L + l f J 1 C (0.0,1) t I/2 

( 1 oo 2 oo I 3 


dy dr 


dy dr 

(4.5.22) 

to (4.5.21) 

(4.5.23) 

the first 

(4.5.24) 


Similarly using (4.5.22) and (4.5.23), the second component of (4.5.19) is 
estimated to give 
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a> 


+ 2 


(«+i) 


|f | C|u I f C (1,0, 2«) t“ 

l <x xx co 3 


•f 


If I + If I 


1 CO 


2 co 


C (1,0,1) t I/2 

3 


(4.5.25) 


The oc-norm is defined as | | a = | S« + I l« • Therefore, in order to calculate 
the oc-semi norm of f j and one needs to estimate l f ^ (e C t^ ) ) - ^ C0(t 2 ) ) | 

and If (q( t i ) ) - ? 2 (q(t 2 ))l for t j * where 6(t) and q(t) are monotone 

functions. Without loss of generality we assume t^ > t 2 - 

The difference between f (.) at two different time measurements is 
expressed as follows. 


f (e(t )) - f^ef^)) = e' (t^ - ^(o.t^ - o' (t 2 ) + 0 t (o,t 2 ) 




- f 1 f K (O.y.t -x) f (u(y,x) ) dy dx 
J 0 J o xx 1 1 

+ f 2 f K (0,y, t -x) f (u(y,x)) dy dx 
Jo Jo xx 2 1 


'o ”0 


„t „i 


- f 1 f K (O.y.t -x) y f (u(y,x)) - f (u(y,x)) [ dy dx 
J rk J r , xx 1 L 2 1 


F o v o 


f 2 [ K (O.y.t -x) y f (u(y,x)) - f (u(y,x)) dy dx 

J 0 J 0 XX 2 2 1 J 


(4.5.26) 


Let I stand for the expression. 


1 = - f 1 f K (O.y.t -x) f (u(y.x)) dy dx 
& J o J 0 ** 1 1 




f 2 f K (O.y.t -x) f (u(y,x)) dy dx 

J r\ J « ** 2 1 


0 * 0 


Jq J ^(O.y.W) - K^O.y.tjj-T) j f^ufy.x)) dy dx 
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- r ^ 


f 1 K (0,y,t -x) f (u(y,x)) dy dx 

{ XX 1 l 

** t ** n 


i ~ Q 
2 


= I + I 
5 1 5» 2 


(4.5.27) 


Making use of the identity property at x * 0, I g 2 is estimated as 


*t . A 


I = [ 1 f K (0 t y, t -t ) f (u(y,x)) dy dT 

5,2 I xx 1 1 


t "0 
2 


A 


f 1 f K (O.y.t -t) 
J. Jn xx 1 


t ' o 
2 


fj (u(y,x) ) - f j ( u(0,t) ) 


dy dr 


Taking modulus on both the sides and estimating the R. H. S. , we have. 


„t 


II I s f 1 f IK (0, y, t — T ) | If | |u (0,t) (y Z /2) | K dy dr 
5,2 J J xx 1 1 cx xx 


t w 0 
2 


. A 


S 2 "“ |f | ( I u I )** f 1 f IK ( 0 ,y, t-x) |y| 2< * dy dx 

1 CX XX 00 J J XX 

t 0 
2 

it - 1 r 

S 2 " (oc+1) If | (lu i ) w C ( 0 , 0 , 2 «) i — 

1 0( XX 00 3 « 


(4.5.28) 


Now the use of mean value theorem modifies I into 

5, 1 


! 5>1 = J* f 1 K xxt ( 0 ,y,s-r) [ f^ufy.x)) - f^uCO.x)) ] ds dy dx 


SO, 


II 


ft A r t 

J s J 2 I f 1 IK (0,y, s-x) | If (u(y,x)) - f (u(0,x))l ds dy dx 
J Q J o J t xxt i i 




s 2 _<x |f | (|u | )“ f 1 f 2 f IK ( 0 ,y,s-x) 1 |y| 2<X dy dx ds 

i « xx CO J J J xxt 


t *0 ~ 0 

2 


t 


s 2 _<x If I (lu I )“ c ( 0 , 0 , 2oc) f 1 f 2 |s - rl *' 2 dx ds (Appendix B 4 ) 

1 « xx 00 4 J n 


r t v O 
2 
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C (0,0,2c*) 

5 Z'* If I (lu |/ — IS - V 

1 « xx 05 «(!-«) 


Let I denote the last two terms of (4.5.26). The estimation 

6 

for I is analogous to I . Indeed, it can be written equivalently as 


6 “ 5 

,t „1 


= f 2 f K (O.y.t -x) - K (O.y.t -x) 

J 0 J 0 L xx 1 xx 2 


y x 


f 2 (u(y,x)) - f j (u(y, x) ) 


dy dx 


.t 


f 1 f K (O.y.t -x) y f (u(y,x)) - f (u(y,x)) 
J J xx 1 2 I 

J t J 0 L 

2 


dy dx 


= I + I 

6,1 6,2 


I is estimated as below 

6,2 




II .I * f 1 \ IK (0,y, t -t) | |y| If (u(y,x))| dy dx 
J* XX 1 2 


t "0 
2 




1 1 J ^xx^’Y’ t i“ x) i iV u{y,x))! d y dT 


t "0 
2 


S ( If 2 '» + S ' 0 ’ 0 ’ 1 ’ “l - V 


1/2 


Ig 1 can ke estimated as follows. Taking modulus on both sides of 




' J 0 2 f o J t x K xxt (0 ’ y ’ s_T) y [ f 2 (u(y>T)) " s (u{y * T)) _ 


we have, 


II 


C (0,0,1) 


. A * f If I + If | 1 — It - t I 

>,1 I 1 CO 2 00 1 1 1 

V J ( 1 / 2 ) (1 - 1 / 2 ) 1 2 


1/2 


* 4 ( I f 2 l„) C , W.O.D It, - t 2 l 


1/2 


(4.5.29) 

technique 


(4.5.30) 


(4.5.31) 


dy dx 


(4.5.32) 
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Collecting the estimates (4.5.28), (4.5.29), (4.5.31) and (4 5 32) 


IMett^) 


y e <y>! a i®'<y - y°-v ~ ® /{t > + * (o,t 


t 2 


)l 


C { 0 , 0 , <x } 

a 

« ( 1 "* CX ) 


if.i* Uu i )“ it - 1 1 

1 xx oo i 2 


+ C (0,0 

b 


(1) f if lL + ly ] It - 1 1 1/2 

V 1 m 2 CO I 1 2 


where C (0,0, <*) = 2~“ C (0,0,2*) + 2' (a+1) (1 - «) C (0,0, 2«) 


and C b (0,0,l) = 2 C^O.0,1) + C 3 (0,0,1) 


Hence, dividing throughout by {©( t^ ) - 6(^)1“ to obtain an estimate 
the «- semi norm of f , we have 


|f i (e(t j )) - f j (e(t 2 ))i 


le(t i > - e(t 2 )| 


I 0, (t ) - iMO.t ) - e' (t ) + * (o,t)| 

i t 1 2 t. 2 

it - t |“ 

1 2 


It - t | 
1 2 


i©(t i ) - 0(t 2 )| 


c ( 0 , 0 , «) 
& 

«(!-«) 


if I (lu i ) 

la xx oo 


+ C (0 

b 


,0,1) | 


lf l» w + ly 


] It - 
<0 j 1 


V 


(1/2) -a 


it, - y 


le(t, ) - e<yr 


or. 


I 0, - ^ 1 

If | s 3LS1 

1 «* ------ .« 


C (0,0, oc) 


(inf le'l) « (1 - a) 


ifj 


1 00 


(lu I J 

XX co 

(inf |0' |)‘ 


c b (o,o,i) 


(inf 


( If | + If | ) It - 

le'i)* 1 lw 2 “i 1 


V 


(1/2) -a 


we write 

(4.5.33) 

(4.5.34) 

(4.5.35) 

for 


+ 


(4.5.36) 
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Towards the estimation of the «-seminorm of f , we have from the equation 


(4.5.19) 


f (q(t) ) = q'(t) - \l> (O.t) - f f K (l,y,t-x) f {u(y,x)J dy dT 
2 1 J„J„ XX 1 


0 “ 0 


f t r i r 1 

K (1 ,y, t-T) y f (u(y,x)) - f (u(y,x)) dy dx 

L K xx 2 1 J 


'0 0 


An analogous procedure provides the following estimate for if. 


iq' - I C (0.0, «) 

F| s + — If I 

2 “ (inf |q' 1)“ «(!-«) 


inf iq' 


C (1,0,1) 

b 

(inf |q'|) C 


,f lL + 


a/2)-« 


(4.5.37) 


In view of (4.5.36) and (4.5.37) we need to take 0 < « < 1/2. Now to find the 
a-norm of f we combine (4.5.24) and (4.5.36). Since (4.5.24) gives 


JL = I 6 ' - * t L + l f J K Hu i )“ C, ( 0 , 0 . 2 c) t 

l"oo u t u oo 1 CX XX 00 3 


+ ( |f iL + lI 


Jcc = I f ii„ + |f i ! « * we have 


(inf |0' | 3’ 


} IS' - 


r c ( 0 , 0 , «) r 

+ . _± 

L «(!-«) L 


t f lu I 

1 + t* (inf |8' 1)“ 2L” 

J I inf 10' | 


‘Vc 


c ( 0 , 0 , 1 ) 

D 


(inf |e'|) 


x j l- 

It 

, n « L i 


V U/21 '“ + tl/2 (lnf ls'l)“ ] ( lf,l„ * lf a l. ) 


(4.5.38) 


Similar expression can also be given for |f j 
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Now, we study the behavior of the coefficient of |fj a in the above estimate, 
as t approaches zero and show that 


C (0,0, «) 

a 

<x (1 - «) 


1 + t (inf |0 




|u I 

XX 00 


inf |0' | J 


and 


C (0,0,1) 

b 


(inf 10' i ) 

for small values of " t" . 


t - t l <,/2> -* - t“ (inf |e 
1 2 


'!)“ ] 


are sufficiently small 


Towards the end we compute a bound for jj u xx I ^ • In view of the solution 
of nonhomogenous equation of (4. 2. 1 )-(4. 2. 3) , we have 


u - ip = f f K (x,y,t-x) f (u(y,x)) dy dx 

XX XX J J XX 1 

0 0 


f [ K (x,y,t-x) y f (u(y,x)) - f (u(y,x)) 1 dy dx (4.5.39) 

Jn in ” 2 1 J 


O "0 


Due to the identity property, this will imply. 


u 


= ^ + J J K xx (x,y,t-x) £ fjufy.x)) - f l( u(x.x)) j dy dx 



r 1 i 

.1 

+ 

j 

K (x,y, t-x) 

XX 


0 J 

O 


r t 

p l 

— 

0 J 

K (x,y, t-r) 

XX 

0 


(4.5.40) 


Calculating the following part only 


y f^Cufy.x)) - x fj(u(x,x)) for j = 1, 2. 


= (y - x) fjUKy.x)) + x (fjMy.x)) - f^uCx.x))) 


(4.5.41) 



104 


Substituting (4.5.41) into equation (4.5.40) and taking the modulus, one gets 


lu | s 1^ | + f f 1 IK (x,y,t-x)| If I My - x) u (e,T)r dy dx 

XX XX J J 0 XX 

r t r l 

+ |K (x,y,t-x)| |y - x| |f 2 (u(y,x))| dy dx 

J o J o xx 

+ [ f |K (x,y, t-x) | I x | |f 2 l w l(y-x) u^fe.xlTdydx 

+ f f |K (x.y, t-x) | |y - x| If (u(y,x))| dy dx 

•*0 o xx 

r t r 1 <x 

+ IK (x,y,t-x)| j x 1 |f ! I (y - x) u (0,x)| dy dx 

J„ J„ xx 1 “ x 


So, this implies, 


|u I £ 1^ I + 

XX XX 


( lf .'« 


+ Ixl If | + |x| If | C (x,0.oc)(|u I ) — 

1 oc 2 « 3 x oo c< 


+ I lf a L + lf ,'») c 3 (x ’°’ 1) 1 ‘ 


(4.5.42) 


Since, in the above estimate |u | is involved, next we obtain an estimate for 

X 00 

I u | in an analogous manner. 

X oo 

u = + f f K (x.y, t-x) f (u(y,x)} dy dx 

x J o J o x 1 

+ J J K^x.y.t-x) y £ f 2 (u(y,x)) - f^ufy.x)) dy dx 

= + | J K x (x,y,t-x) f^uty.x)) - f^ufx.x)) dy dx 

+ f f K (x,y,t-x) x [ f (u(y,x)) - f (u(x,x)) dy dx 
J o J o L 2 2 J 

+ f f K ( x .y. t-x) (y - x) f (u(y,x)) dy dx 
J 0 J 0 2 

J J K x (x.y.t-x) x £ f i (u(y,x)) - f^ufx.x)) j dy dx 
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" J J K x (x,y,t-T) (y - x) f^uCy.T)} dy dx 


* 

’0 ^ 0 

This implies from Appendix B (B4), 


IV 5 1^1 + l^l* C 2 (x, 0 , 0 ) t 1/2 + If^ C 2 (x,0,l) t 


|f iL C 2 (X '°’ 1) 1 + ( |f 2 L + |f iL ) w c 2 (x ’ 0 ’°) tl/2 ( 4 - 5 -43) 


Substituting (4.5.43) into (4.5.42), one has an estimate for |u | 

XX » 

explicitly. Therefore, for small values of "t" and for some constant E 


If I + If | s E and |f | + |f | s E 

l ot 1 1 oo 2 oc 2 00 


then 


u — > tp like 0(t a/2 ). 

XX XX 


Hence 


f lu 1 ] 

XX 00 

<X 

1* I ] 

XX oo 

OC 

g£ 

r lu' ' (x) | ' 

k inf 10' 1 J 


. inf |0' 1 j 


inf |0' (0)| j 


So considering the flat initial value 


C (0,0, «) 


OC (1 - oc) 


1 + t (inf |0 




f lu I 

XX oo 


inf |0' | 


will become 


C (0,0, oc) 
a 

oc (1 - oc) 


' u"(x) 

0 


( inf le' (0) | J 


and 


C (0,0,1) 

D 




(inf |0' | ) 

Therefore, (4.5.38) shows that 


t | <1/2)_0C + t 1/2 (inf 
2 


le' I )“ ] 


JfJ s E and |f 2 i * E. 
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Hence 


T tf ,f ] (t) : B c C* [O.TJ x c" [0.T] — » B c c“ [0,T] x c“ [0,T] I 

Of (3 1 2 E ■ 


Finally, we note that T Q is a continuous map from c“ x c“ into itself 

tf,q 

This follows from the continuous dependence of the solution "u" of the direct 
problem on the right hand side function f. 


LEMMA 2 : The set 


V = 1 (f 1 ’ f 2 > : fi f Jco - E - I f Jcx 5 E - f ° r 1 = 1. 2 


is closed in C^xC* for 0 < £ < c* < 1. 


PROOF : Clearly, V is a closed set in c” x c'*. 


Let for i = 1, 2, 

| f " | be a sequence in V such that f” — » f ^ in x cP for £ < 


We claim that f e V. 

l 

Since f” — > in cP x cP 


* « f ; - f Jeo 


so I f 


i L - 1 f X + K - f J M 


f 1 1. - lim 


f iL + l f " • f iL 


Further , 


l f ( tx) - fjCy)| | lim (f“(x) - f“(y)) 


x - y 


x - y 
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- f"(y) 

lira ! 1 

x - y 


s lim E = E 



Hence V is 


closed in C* x C* 


I 


THEOREM 2 (EXISTENCE) : Assume that T is a self mapping from V into itself, 

C7, q " 

where V = | (f jt f g ) : | fj^ £ E, | fj^ £ E, for i = 1 , 2 | . Then T 0 ^ 

has a fixed point. 


PROOF : The set V is a closed, convex set in C a x C°\ By the above lemma V 

is also closed and convex set in (P x C^. Since V is a bounded set in C* x 

C°\ T V is precompact in C^xC* (FRIEDMANN, P.188 [49]) . Hence T„ is a 

0,q 

completely continuous operator. 

Therefore, considering T : V c x — » V c x cP , and by 

SCHAUDER’s fixed point theorem (FRIEDMANN, p. 189 [49]), we conclude that T 

«,q 

has a fixed point in cP x Cp. | 

REMARK : 

The existence of a limit of the iterated sequence -|f| k) j- is not 

guaranteed by the existence theorem obtained above. In the next chapter we 

shall prove a centred contraction mapping Property of the T map under some 
— 

additional conditions. This will be sufficient to show the convergence of the 
iterates if^l and hence establishes the convergence of the iterative scheme 


outlined in this chapter. 



4.6 CONCLUSIONS : 


The question of existence of the solution of a nonlinear source function 
problem of a fixed point of a Heat equation in a bounded spatial domain is 
answered in this chapter. The problem is shown to be equivalent to the 
existence of a fixed point of a certain nonlinear map. An iterative scheme is 
outlined via an operator what is shown to be self mapping on a bounded set in 
the Holder space. A large number of domain estimates have been computed 
explicitly to arrive at this result. Subsequent use of the SCHAUDER’ s fixed 
point theorem ensures the existence of the original source problem. 



CHAPTER V 


UNIQUENESS OF THE SOLUTION OF THE HEAT CONDUCTION SOURCE 
PROBLEM AND CONVERGENCE OF AN ITERATIVE SCHEME 


5.1 INTRODUCTION : 

In this chapter the question of uniqueness and the convergence of 
iteration scheme for the inverse Heat conduction source problem is addressed. 
This Is done by showing a centered contraction property for the map T . 

~ o,q 

DEFINITION : Let X be a Banach space and T is a map from U c X to V c X. Then 
T is called centered contraction map about u q e U if 

fTu o - Tu| s A |u Q - u| 

with V u e U and A Is a constant independent of u, s. t. 0 £ A < 1. 

Suppose that (f ,f ) is a (0,q) -fixed point of T Q 
i. e., [ f l( e(t», f 2 Cq<t» ] - T 0>1 [ q, f, ](« (5.1.1) 

and let (g ,g ) and (g ,g ) are related by 
1 2 1 2 

[i/ett)) , i 2 (q(t)) ] -I 8i< [g,.*,] (t) 


(5.1.2) 
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We show below that T is a centered contraction map in a neighborhood 

© > <3 

of the fixed point (f ,fj. Therefore, we need to estimate f $ - ^ and f - g 
in terms of f - g and f - g , and obtain the following estimate : 

It 2 2 

I L " I T 0,q (f r f z J - T e,, t V*i ) l. 

s A f (f t . f 2 ) - (g s .g 2 ) | w (5.1.3) 

we note that 

I " T e„ (g 1 ’ 8 . )| « = 1 ( W ' (i .' i 2 ) I- 

= ( I f . - i. I« * I f 2 - 8 z !« ) (5 - 14) 

where 


I T e., (f ,>V - = I WW ' WW' 1 


00 


+ 1 T o., (f .- f 2 ) - T e., <g .’*2 )l « 


(5.1.5) 


In section 2 the main theorem is stated. The sup and the semi-norm for 
the operator are calculated in sections 3 and 4 respectively to prove the map 
to be contraction and hence the convergence of the iteration scheme obtained 
in the precious chapter. This chapter ends with a conclusion in section 5. 


5.2 THE THEOREM : 


In this section we state some additional hypotheses which need to be 
satisfied in order that the main theorem regarding centered contraction 
property can be proved. 

Note that if f is a uniformly Lipschitz function and g e c" then, f ° g 6 
C* and 



Ill 

I f ° S L * I f L + i f f, Igl* (5.2.1) 

If f © g vanishes at some point then also 

I ' * * I. s C Ifl, Igl. = C |f I gl a (5.2.2) 

for bounded domains f € Lip class preserves C* regularity under composition. 
This is not sufficient to show uniqueness. Hence one extra property is needed. 

PROPERTY S : A function f € Lip class has the "property S" if given functions 
u, v € C°\ the mapping u — » f(u) is C°\ that is there exits C < » and 

S 

J f(u) - f(v) I* * I u - v | a . (5.2.3) 

Below it is assumed that f^ will satisfy the property S. 

THEOREM : If the over posed-boundary value problem (4.2. 1) — (4.2.5) defined 

in chapter IV possesses a (6,q)-flxed point (f ,f ) with property S, then it 
is unique. Furthermore, the operator T is a centered contraction on C* x 

t),q 

C“ for « < 1/2 under the assumptions (B1)-(B3). 

5.3 SUP-NORM ESTIMATION : 

Let ufx.tjf^f ) and v(x,t;g i> g 2 ) be the respective solutions of the 
following partial differential equations. 


u (x,t) - u (x,t) = (1 - x)f (u(x, t ) ) + xf (u(x,t)) + y(x,t); (5.3.1a) 

t XX 12 

X € (0,1), t > 0; 

u (0,t) =0, u Cl, t) = 0 ; t > 0 ; (5.3.1b) 

X X 

u(x,0) = u (x) ; x € [0,1] ; (5.3.1c) 

o 

u(0,t) = 0(t) ; u(l,t) = q(t) ; t > 0; (5.3.1d) 
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v t (x,t) 


- v (x, t ) = (1 - x)g CuCx.t)) + xg (v(x,U) + 7 (x, t ) ; (5.3.2a) 

XX 1 * 


v (0,t) =0, v (l,t) = 0; 

X X 

v(x,0) = u q (x) ; 


x e (0.1), t > 0; 
1 > 0; 
x € [0, 1 ] ; 


v(0,t) = 0(t) ; v(l,t) = q(t) ; t > 0 ; 


(5.3.2b) 
(5.3.2c) 
(5. 3. 2d) 


The difference between the solutions of the nonhomogeneous equations 
(5.3.1) and (5.3.2) is given by 

u - v = f f K(x,y,t-T) F(1 - y ) (f ^ (u) - g^v)) + y(f 2 (u) - g 2 (v)) dy dr 
•*0 ^ o 

(5.3.3) 

Where K(x,y,t) is the Green’s function for the Neumann problem. Now 
differentiating equation (5.3.3) twice with respect to x one obtains. 


u - v 

XX XX 


J J K xx (x,y,t-X) [ ( 1-y ) •( fj(u) - gi (v) 


f “ 

+ y f (u) - g (v) dy dr 
2 2 / . 


(5.3.4) 


According to the formula (4.3.6) in chapter IV, equations (5.1.1) and 


(5.1.2) can be written as 




6' (t) - u (0, t; f , f ) - y(0,t), 

XX 12 


q/ U) " U xx (1 ’ t:f l’ f 2 ) " 


(5.3.5) 


[*• • «i] =T e„[*. ■ «.] - 


[ e ' (t) " V xx ( °’ t;g l’ g 2 ) " 

q,(U ~ V xx (1 ’ t:g 1 * g 2 ) " ] 


(5.3.6) 


Therefore , 
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e,q 


V f 2 


e,q 


gj ► g 2 


' . f . ' *1 ' { 2 - h 


* [ v „ (0 > ti V< ! 2 ) - \ x <°’ t;f vV' 

V xx ,1,t;g i' g 2 ) - U xx <1,t;f l' f 2 ) _ 


Substituting (5.3.4) into (5.3.7), we get 


- S, * f 2 " i 2 


[ f I ' 1 

L v o v Q 


K (0,y,t-x) 

XX 


f f K (l,y, t-x) 

•*0 •'o xx 


(1-y) | gj(v(y,T)) - f^uCy.x)) j 

+ y | g 2 fv(y.T)) - f 2 (u(y,x)) j j dy dx, 
(1-y) | gj fv(y,x) ) - f^ufy.x)) j 

+ y | g 2 (y(y,x)) - f 2 (u(y,x)) )]] dy dx 


We first of all obtain an estimate of the sup-norm of f - g 5 
that f - g 2 can be estimated analogously. 

fj ~ g x = J J K xx (0,y,t-x) | g^vCy.x)) - f^uly.x)) j dy dx 


u K^CO.y.t-x) y | | g 2 (v(y,x)) - f 2 (u(y,x)) j 


^ g^vCy.x)) - fj (u(y,x) ) j | dy dx 


= 1+1 
1 2 


Adding and subtracting some suitable terms, I as given below can be 


as a sum of I and I 

1,1 1,2 


(5.3.7) 


(5.3.8) 

and note 


(5.3.9) 

rewritten 
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I = r f K x (0,y,t-T) f gj (v(y,T) ) - fMuCy.t)) 
1 h J o xx ^ 


dy dr 


= f f K . (0,y,t-x) I" gj (v(y,r) ) - Mvjy.x)) 

•*0 xx 


dy dx 


n 1 r l 

K (0,y, t-x) f (v (y , x ) ) - f (u(y,x)) dy dx 

XX 1 * 

0 L J 


V + *1.2 


(5.3.10) 


Consider ^ and e C . Making use of the identity property, 1^ ^ is 
estimated in the following manner. 


I - f f K (0,y, t-x) 

i.i J 0 J 0 « 


g 1 (v(y,x) ) - fjfvfy.x)) - 


g 1 (v(0,x) ) + f^vtO.xl) 


dy dx 


Uj jl - J | |K xx (0,y,t-r) i |f j - i v (y , x ) - v(0,x)i w dy dx 

s |f i -■«,!« dvJJ* j‘ J IK„«0,y,t-r)| (y z /2)“ dy dr 


a lf l ‘ B i‘« (|V xxU“ 2"'“*' — C 3 (0,0.2«) t“ 


(5.3.11) 


Where C g ( is defined in Appendix B (B6). Noting that f satisfies the 
property S with constant and that 1^ ^ is equivalent to 

1 i,2 ~ ~ ^ fjluCy.x) ) - f^vCy.x)) j dy dx 

r 1 r 1 r 

■ J o J o K xx (0,y,t * T) [ - fJufO.x)) j dy dx 
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we have 


‘Vs 1 " 


J J !K xx (0,y,t_T) I ( y lf,(u(.,T)) - f^vL.x))^ | dy dr 
C . J J tK xx C0,y,t ~ T)i ( y " v{. ,t) ^ J dy dr 


< C f f |K (0,y,t-r)| 

s 1 J XX 
J 0 J Q 


y I u ( . , t ) - v (. ,t)| 

X X to 


dy dr 


sC | u - v | C (0,0,1) |t| 

s x x CO 3 


( 1 / 2 ) 


(5.3.12) 


To compute the first part of (5.3.8), we need to estimate I 2 , where 


[ = f f K (0,y, t-x) y 

2 J o J o xx 


\ 

g 2 (v(y,r)) - f 2 (u(y,r)) dy dx 


L J>. (0, y, t-x) y | ^ (u(y,x) ) - g j (v(y,x)) j dy dx 


o 'o 


= I + I 

2,1 2,2 


(5.3.13) 


Once we have an estimate of I , the estimate for I will follow 

2,1 2,2 

similarly. Adding and subtracting yf (v(y,x)) inside the integral sign, I 

2 d, 1 

written as sum of two terms. 


* 2,1 = ~ J J K xx (0,y,t-x) y | f 2 (u(y,x)) - f 2 (v(y,x)) j dy dx 
-u K xx (°, y ,t-x) y | f 2 (v(y,x)) - g 2 (v(y,x)) j dy dx 


= I + I 

2 , 1,1 2 , 1,2 


( 5 . 3 . 14 ) 
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Taking modulus on both the sides of 5 2 have. 


|I | a f f |K {0,y, t-T ) | |y| If ( v (y, x ) ) - g (v(y.T))| dy dt 

2.1,2 J Q J 0 xx 2 2 

nl 

IK CO, y, t-x ) I |yl dy dt 

« XX 


* if - gj C (0,0,1) t 

2 2 oo 3 


a/ 2 ) 


Let f € Lip class with constant C, then 


(5.3.15) 


I . .1 s f f IK (0,y, t-x) | |yj |f (u(y,T)) - f (v(y,x))| dy dr 


S C C (0,0,1) | u — v | t 

3 oo 


(1/2) 


(5.3.16) 


Estimating in a similar manner and taking Into account that f satisfies 

1 

both property S and is in Lip class with constant D, we have the following 
estimate. 


1 1 2,2 I S |f i " Sjoo C 3 (0,0,1) t (1/2) + D C 3 (0,0,1) | u - vl^ t Cl/2) (5.3.17) 

Now putting together all the estimates (5.3.11), (5.3.12), (5.3.15), 

(5.3.16) and (5.3.17), the corresponding estimate for f - g is given by 

i t 


|f i " SiL s |f , - g,L Civ I )“ 2 

1 1 00 1 1 <X XX CO 


~(«+l) 


C 3 (0,0,2«) t 


cx 


+ c s |u x ' V X L c 3 (0 »°>n t (1/2> 
+ |f 2 “ 82 L C 3 (0,0,1) t (1/2) 

+ c C (0,0,1) lu - v| t (1/2) 

J CO 
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+ 1 f i - 8,1. ^( 0 , 0 . 1 ) t 


(1/2) 


+ D C (0,0,1) | u “ v | t 

3 co 


( 1 / 2 ) 


Analogously, 


Y - i*'» = l f , - *,'« — c 3 (‘.o.w f 


+ C |u - V } C (1,0,1) t 

s x x oo 3 


(1/2) 


+ if ? - C (1,0,1) t 

2 2 oo 3 


(1/2) 


+ C C (1,0,1) |u-v| t 

3 oo 


(1/2) 


♦ lf , - *,L C 3 (1.0,l) t 


( 1 / 2 ) 


+ D C (1,0,1) | u — v | t 

3 oo 


(1/2) 


(5.3.18) 


(5.3.19) 


To find the sup norm of T , one can add up (5.3.18) and (5.3.19). Hence 

b,q 


l T e, q ( f,’ r 2 ) - Y„ ( VS 2 >L " l f , - *,l. * lf 2 - «*'. 


* J C 3 (0,0,2o<) + C 3 (1,0,2oc) 


{ 

{ 


} if, - g . 


i (iv f r-i— 

oc XX oo -OC+1 

<x2 


+ i C (0,0,1) + C (1,0,1) C |u - V | t 

* 3 3 j s x x co 


{ 


,0,1) J. 

, 0 , 1 ) | 

+ | C 3 (0,0,1) + C 3 (1,0,1) | C lu - vl w t 

+ | C 3 (0,0,1) + C 3 (1,0,1) | If i - g 1 l (0 t 

+ | C 3 (0,0,1) + C 3 (1,0,1) J. D |u - v| M t 


( 1 / 2 ) 


+ -( C 3 (0,0,1) + C 3 (1,0,1) V |f 2 - g 2 l tt t 


(1/2) 


( 1 / 2 ) 


( 1 / 2 ) 


(5.3.20) 
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In order to get the required result, one has to show, I - v } 

{ | v | ) w , and |u - v| are small quantities for small values of 

xx eo ® 

"t". Differentiating the equation (5.3,3) with respect to x, one finds 


11/ \ 
ufc.t) - v x (x,t) = J J KJx.y.t-x) | f j (u(y, x ) ) - g^vty.x)) j dy dx 


r 0 *0 


n l 

K (x.y.t-x) y f 2 (u(y,x}) - g^vly.x) ) 

n V 


'o “o 


f (u(y,T) ) + g. (v (y , x ) } dy dx 
1 1 


= 1+1 
3 4 


(5.3,21) 


I 

3 


can be expanded as 




I = f f K (x.y.t-x) 

J n J A 


F 0 •'O 




f K (x.y.t-x) 


0 v 0 


= I + I 

3,1 3,2 


^ f i (u(y,x)) - f } (v(y,x ) } j dy dx 
f f i {v(y,r)) - g^vfy.T)} | dy dT 


For f e Lip class with constant D, I 

l 3,1 

U ,.! 1 s D lu - vl o j‘ |K x (x.y,t-r)| 

£ D |u - v| C (x,0,0) t 1/2 
» 2 

and I 3 2 can be estimated directly as 

S l f i " f f |K (x,y,t-T)| 

J o J o 


is estimated as 


dy dx 


dy dT 


(5,3.22) 


(5.3.23) 


s l f r g il«o c 2 (x '°* 0) tl/2 


(5.3.24) 
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Hence combining (5.3.23) and (5.3.24), the bound for I is given by 
I I 3 I 5 ( D lu - vl. * If, - 8,1. ) C 2 (x.0.0) t I/2 


We break I into two parts as follows 

4 


W f K^Cx.y.t-x) y | f 2 (u(y,x)) - g 2 (v(y,x))j dy dr 

* 

+ f f K^fx, y, t-x) y 


'o o 


I + I 

4,1 4,2 


gjCvfy.x)) - f (u(y,x) ) 


dy dx 


We estimate I and I in the similar as way it is done for I, 

4 , 14,2 2 


■ I 1 K^tx.y.t-x) y | f 2 (u(y,x)) - f 2 (v(y,x))j dy dx 


.t ,i 


u KMx,y, t-x) y J f 2 (v(y,T)) ~ g 2 (v(y,x))J 


I 4- I 

4 , 1,1 4 , 1,2 


dy dx 


For f € Lip class with constant C, the following estimate for I 

2 1 


I I 4 , 1 , 1 l - c l u - V L J J I K x (x,y,t-r)| | y | dy dr 


sC lu - v| 


7 

0 0 

7 

0 J o 


r t 1 

J J j K x (x,y,t-x)| |y - x| dy dx + 


O i 

| K x (x,y,t-r)| |xj dy dx 

u 0 X 


s C (u - v| w | C 2 (x,0,l) t 1/2 + C 2 (x,0,0) |x| j 


. 1/2 


(5.3.25) 


(5.3.26) 

Hence , 


(5.3.27) 

is given 


(5.3.28) 
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and similarly the estimate for I 41 2 is 8 ivcn b V 

|i, ili2 l - |f 2 -8 2 l. fj'j i'’! 

* \ / ■;> 

< If - g | C (x,0* 1 ) t + c (x f 0.0) 1 x j t " (5.3.29) 

* 2 2'to ^ 2 2 j 

So in view of (5.3.28) and (5.3.29), I is estimated as 


I 1 , ,! s ( c,(*.0.D t ,/2 ♦ c 2 (x.o,o> |x| J ( C |U - v| m - |f 2 - g 2 | ai j t 1/2 

(5.3.30) 

Similarly, the estimate for I is given by 

* * 2 

|I 4(2 I * ( C 2 (x,0,l) t 1/2 + C 2 (x,0,0) |x| j ( D j u — vj^ + | f t - gi | ro ] t V2 

(5.3.31) 

Adding (5.3.30) and (5.3.31), the estimated value for 1^ is computed as 
|IJ =£ [ C 2 (x,0,l) t 1/2 + C 2 (x,0,0) jx| ] 



Hence combining (5.3.25) and (5.3.32) one obtains the required estimate 
for lu - v I 

1 X X 1 00 

|u - v | 

1 X X 00 

- [ [ c 2 (x,0,l) t 1/2 + C 2 (x,0,0) |x| | (C + D) + D C 2 (x,0,0) J |u - v| w t 1/2 
+ ( C 2 {x,0,1) t 1/2 + C 2 (x,0,0) jx | + C 2 (x,0,0) J jf i - gj^ t 1/2 
+ ( V*' 0 ’ 15 t 1/2 + C 2 (x,0,0) |x| ] |f 2 - t 1/2 
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A ju - v| t 1/2 + A [f - g ! t 1/2 + A If - a I + 1/2 

1 ' 2 1 1 *1'® 3 1 2 g 2 'co 1 


1/2 


SA , I ” ' v l . * a 2 


( l f . -8.1. * l f 2 - Sj. ] 


, 1/2 


(5.3.33) 


Now we obtain an estimate for u - v . Indeed 

00 


u(x,t) - v(x,t) = f K(x,y,t-x) f f (u(y,x)) - g (v(y,x)) ) dy dx 

J o J o ' 1 1 J 


.t ,i 


f f K(x,y,t-x) y [ f (u(y,x)) - g (v(y,x)) - 
J o J o v 2 2 

f j (u(y,x) ) + g i (v(y,x)) j dy dx 


= 1+1 
S 6 


(5.3.34) 


Since the very structure of u - v is similar to u - v , instead of the 

X X 

kernel which is of first derivative of x, the estimate for (5.3.34) given 
similar to (5.3.33) except some changes in the domain constants using the 
Appendix B (Bl). 

,t „i 


n |K(x,y, t-x) | dy dx s C (x,0,0) t , and 

v o * 

f f jK(x,y, t-x) | jy - x| dy dx £ C (x,0,l) t 

J r \ Jn 


0 0 


Hence 


I \i — v I s A |u - v| t + A If - g I t+A If - g I t 

1 *00 4 * *co 5*1 & l‘co 6 1 2 2 1 oo 

S A , I" - V L ‘ + A S ( l f l * g .l. * l f 2 ‘ g J. ) ‘ 

’ |U ' V| ” ^ 1 -l t i |f ' * S,l “ + ^ ' g2, “ ) 1 


= \ ( i f , -8.L* |f 2 -8 2 l.) ‘ 


(5.3.35) 
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Substituting (5.3.35) in (5.3.33), one obtains 


u - v s 

1 X x 1 00 


A A t + A 

17 2 


8, „ + i f 2 


8 


. 1/2 


2 oo 


(5.3.36) 


Finally we estimate ( |v y* ■ Utilizing the representation, one can write 


,t ,1 


= i (i + f f K (x.y.t-x) g (v(y,x)) dy dT 

« L Jo ** 5 


'0 “0 
.t „1 


J J K „ (x - y ' t - T) y 

* A » A 


'0 "0 


g 2 < v (y.T)) - gj ( V (y , T ) ) 


dy dt (5.3.37) 


Where \p is the comparison function. Similar procedure can be adopted to bound 
(5.3.37) as it is done for ( |ii y* in chapter IV to obtain analogous 


estimate. 


5.4 SEMINORM ESTIMATION : 

In order to complete the «-norm, in this section we obtain an estimate 
of the «-seminorm of the map T . According to the definition it can be 

t),q 

written as 





lf >- 


g, I + If. 


1 oc 


Sol 
2 » 


(5.4.1) 


We estimate If - g I , the estimate for If - g | will follow similarly. 

1 1 2 2 oc 

Without loss of generality assume t > t . Then 

i 2 


If - g | = 

l 


sup 

eu )*eu > 
1 2 


[ yect^) - g 1 (0(t i )) J - £ yety ) - g^ety)] 


lo(t i ) - e(t 2 )r 


( 5 . 4 . 2 ) 
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It is given that 




’ f 2 ' 8 2 


f f K^(0,y, t-x) I ( 1 - y) 
-o ■'o xx 


| gjfvfy.T)) - fjuly.x)) j 

4 y | g 2 (v(y,x)) - f 2 (u(y,T)) J. J dy dx , 

f f K ( 1 ,y, t-x ) (1 - y) - g (v(y,x)) - f (u(y,x)) } 

J o J o xx L 1 1 J 

4 y | g_,(v(y,x)) - f 2 (u(y,x)) j dy dx j 


(5.4.3) 


In order to get the { f ^ - g 1^, our aim is to compute the numerator of 

(5.4.2). The limits of the integrations are divided into two parts, from 0 to 

t and from t to t . Thus, 

2 2 1 


i 7 = f j(®(t 1 )) " ijOttj)) 


f (e(t )) - g (e(t )) 
12 1 2 


= f 1 f K (O.y.t - x) 
J o J o xx 1 


g^vCy.x)) - f i (u(y,x) ) dy dx 


J^ 1 J^KjO.y,^- x) y | g 2 (v(y,x)) - f 2 (u(y,r)) j - 

* gj(v(y,x)) - f i (u(y,x)) j j dy dx 

J 2 J K^O.y, t 2 ~ x) £ gj (v(y,x) ) - f^uty.x)) j dy dx 


' I o 2 J o K xx to - y -V T) y [ { 8 2 (v(y - T)) - f 2 (u(y - T)) } - 

| gjfvfy.x)) - f^ufy.r)) J dy dx 
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Now the R. H. S. of the above expression can be broken into two parts with 
respect to the limits of Integrations. Indeed. 


r t r i r 

= I 1 I K (0,y, t - t) g (v(y,r ) ) - f iuiy.r)) dy dx 
7 J J XX 1 1 1 j 


f 2 f f K (0,y,t - x) - K (O.y.t ~ t) 

J XX 1 XX 2 

J 0 J 0 L 


gjfvCy.r)) - f jfuCy.r) J dy dT 
f 1 f K (O.y.t - x) y f g (v(y,x) 5 - f (u(y,x)) 

J t J 0 xx 1 L 2 2 


gj (v(y,x) ) + f i (u(y,x ) ) dy dr 


i s r [ 

J o J o L 


K (O.y.t - r) - K (O.y.t - t) y x 

xx 1 xx 2 


g_(v(y,T) ) - f (u(y,x) ) - g (v(y.r) ) + f (u(y.r) ) dy dx 
& * # 


(5.4.4) 


These four parts can be represented as I g , 1^ , I and 1^ respectively. Now 
we proceed to estimate each of these terms : 


I 8 ~ J 1 J K xx (0,y,t r _ fj (u(y,T) ) dy dx 

*2 ° 

= | J K xx* 0,y,t i~ [ gj( v (y»T)) - f i (v(y,x)) dy dx 
+ J 1 J K xx (0,y,t r T) [ f ! (v(y ’ T)) - fjtufy.x)) 1 dy dx 




( 5 . 4 . 5 ) 
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Let f , g € C . By the application of the identity property I can be 

1 8,1 

rewritten as : 




I = f 5 f K (O.y.t - t) 
8,1 J, ** 1 


t 0 
2 


| gjCvCy.T)) - (v(y, t) ) | - 

| gjfvCO.x)) - fjvfo.x)) | J dy dx 


This will imply, 


'W S If ! ’ g l'« "'Jj" ” 2 ' { “ +n 'S - 'J 


(5.4.6) 


To estimate I use the fact that f satisfies the property S. Making use of 

8,2 1 

the identity property, I Is first of all written as 

8,2 




l s,2 = I 1 J K xx (0,y,t r T} fj(v(y,x)) - f^uCy.T)) - 


f (v(0',r)) + f^uCO.x)) dy dr 


Following the arguments utilized for the estimation of I and Appendix B 

1,2 


(B5), we estimate I as 

8,2 


II J * 
8,2 


| 1 J |K^(0,y, t^- t)| | y If^uC-.x)) - f^vf-.x))! J- dy dx 


„i 


s C f 1 f |K (0,y, t - x ) | | u( . , x) - v(.,x)| lyl dy dx 

«J, J„ ** i 1 


t ■'o 
2 




C f 1 f IK (O.y.t - x ) | | u (. ,x) - v (. ,x) I lyl dy dx 

sj j xx 1 X x 00 


t *0 
2 


* C lu - v | C (0,0,1) It - t \ 

8 X x CD 3 12 


1/2 


(5.4.7) 


From (5.4.6) and (5.4.7), the estimate for I is 

8 
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Ills If - gl ( tv I )“ -L 2- 1 ”-" C (0,0, S.) It - II 

1 o ' 1 l CK XX ® <* J I * 


4 c 111 - V i C CO, 0,1) it - t I 

m X x CO 3 1 2 


1/2 


(5.4.8) 


I can be estimated In a similar fashion, as it Is done for I Indeed, 


10 


Jt. .1 


a 

I = I K (0,y,t - T ) y 

io J, J„ xx 1 


t "o 
2 


g (v(y,T ) } - f (u(y,T ) ) 

m* 


dy dx 


„t „i 


- r j K „ ( °’ y ’V i > y 

** *. ** n 


t -'o 
2 


gj ( v (y , x 5 ) - f j (u(y,x ) ) 


dy dx 


= I + I 

10,1 10,2 


(5.4.9) 


Since we can not apply the Identity property at x = 0, I ^ and 1^ ^ , are 
estimated directly. 


10,1 


- [ 1 f K (0,y, t,- t) y 

J t J o 


XX 1 


2 


- [ 1 f K (0,y, t - t) y 

J , K x * 1 


t *0 
2 


f ( v ( y , T } ) - g ( v ( y , T ) ) 
2 2 


f (u(y,x) ) - f (v(y,x)) 
2 2 


dy dt 


dy dx 


= I + I 

10 , 1,1 10 , 1,2 


(5.4.10) 


I is estimated as, 
10* 1,1 


ft „1 


1 1 s lf a ■ J 1 I IK lot (0 ’ y ’V t) 1 lyl dy dT 


* lf 2 - 8 2 l„ c 3 t°.0.1) |t, - t 2 l 


1/2 


(5.4.11) 


^2 e ^Ip class with constant C. Then the bound for I is given by 


10 , 1,2 


II 


10,1,2' 


■ r * r 
\ ° 


M l f 2 (u(y»T)) - f (v(y,x))| dy dx 
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From (5 


II 


10,1 


s C 



IK (0,y, t - r) J 

XX 1 


lyl 


lu(y.t) 


2 

s C C (0,0,1) |u - v| |t - t | 1/2 
3 eo l 2 

.4.11) and (5.4.12), the estimated value for 

' = If 2 - IS - S'" 2 

4 c c (o.o. n iu - vi it - t i 1/2 

3 oo 1 2 


- v(y,r)| dy dr 


I 10( , is given by 


(5.4.12) 


(5.4.13) 


For fj e Lip class with Lipschltz constant D, 1^ ^ is estimated in the same 
way as it is done for I J0 ^ So the result can be written as 


|I | * |f - g ( C (0,0,1) |t - t | 

10,2 1 00 3 1 2 


1/2 


+ D C (0,0,1) lu - v| |t - t i 
3 to 1 2 


1/2 


(5.4.14) 


Hence from the results that are obtained in (5.4.13) and (5.4.14), the 
estimate for I is given by 


S lf 2 - g 2 l„C 3 (0.°.l) It, - t 2 l ,/2 
* C C (0,0,1) |u - vl It, - t l w 

3 CO 1 2 

+ D C 3 (0,0,1) | u - vl It, - t 3 l ,/2 (5.4.15) 

The expressions for I and I , which are almost similar can be estimated by 

9 11 

the application of Mean value theorem. I is written as 
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I « 

9 


I 


t .1 


K (O.y.t - O - K (O.y.t - x 


{v(y,T ) ) - (ufy.x ) ) 


dy dr 


-t ,t _ ,i 


- J 1 I 2 I ’w' 0 -*-*- 1 ’ 


t * 0 *0 

2 


( v (y, T ) ) - f } (ufy, t ) ) 


dy dr ds 


,t ,t „i 


= f 1 f 2 f K (O.y.s-t) g^vCy.T}) - f^vfy.x}) 

J + n J n *• 


0 *0 

2 


dy dr ds 




rri K „t <o -^ s - T) 

J . J A A 


't rf 0 *0 
2 


(v{y,T ) ) - f ^ (u(y,T ) ) 


dy dx ds 


= I + I 

9,1 9,2 


£5.4.16) 


Let f , g i e C , then application of the identity property will give us the 
following equivalent expression for I 

9, 1 


*t 


I ».> = I. , L 2 L K -t (o - y - s ' T) 


t 'O ■'0 
2 


gj £v(y,x) ) - f s £ v(y , x ) ) 


- g (v (0 , T ) ) + f (v(0,t)) 
1 t 


dy dx ds 


Taking modulus on both the sides and estimating, we have 


_ i*l 


~ J ^ 1 J 2 J lK xxt (0,y,s-T)l Ifj - |v£y,x) - vfO.TlT dy dx 


ds 


£ 2 |f i " S,L ( > v I >* 

1 1 CC XX CD 


rS r^2 r 1 

O 0 


|K xxt (0,y,s-T)| lyf 2 * dy dr ds 


s 2 




|f i " g i'oc UvJJ* C 4 (0,0,2 oc) J 1 J 2 Is - tT - 2 dr ds 

* " (from (B7)) 


t -o 
2 
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iy 5 'y, 1 * iI 9 .2 l 


. c* 4 


C (0, 0, 2» ) 


s ^ ifj - 8j t a uvjj & ( i - oi i 2 ' 


it. - t. 


+ 4C | u - v | C (0,0,1) it - t J 

g x X ® ^ 


1/2 


(5.4.19) 


Similarly, using the Mean value theorem, the equivalent form for I i s 

expressed as 


„t „t „i 


1 = f 1 [ 2 f K (O.V.s-x) y g 2 ( v ( y . t ) ) 

11 J t j o J o xx 
2 


f (u(y,x) ) 
z 


dy dx ds 


„t „t A 


- j ‘ I 2 | K x «‘ 0 ' y ’ S - T! * 

J *. J r \ J rv 


t *'0 "o 
2 


g £v(y,x) 5 


f } (u(y,x) 3 


dy dx ds 


= I + I 

11,1 11,2 


(5.4.20) 


I is estimated first. The result obtained is on same line as I . With 

11,1 11,2 

suitable choice of terms, 1^ j is written in an equivalent form. 


„t „t ,i 


y, ■ - J. 1 f 2 f y 


t ~0 "0 

2 


fJv{y,T>) - g CvCy*x) ) 

<c 2 


dy dx ds 


„t „t a 


I t ‘ I/ J„ K ^ (0 -5 r ' s - t) y 


f 2 (u{y,x) ) - f 2 (v(y,x)) 


dy dx ds 


= I + I 

11 , 1,1 11 , 1,2 


(5, 4.21*3 


I is estimated as 
1* * * * * 


nt -t „1 


^n, 1,1* £ J o J * K xxt £°»y* s “ T )l lyl lf 2 (v(y,T)) - g 2 (v{y,T))| dy dx ds 
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£ |f “ g I 

2 2 ® 


0 * 0 


m * 

fK **t C0,y>s ' T)I ly i d y dT ds 

* n 


3 4 if - g I c (0,0,3) it - t | 

2 °2 ® 4 1 2 ‘ 


1/2 


(5.4.22) 


If f € Lip class, then the bound for I is given by 

2 11 , 1,2 


f t 

5 j j j 1 K x X t ( °’ y,S " r)1 !yl !f 2 (u(y * T ^ - f 2 (v(y,T))| dy dx ds 
l 2 ° ° 


3 C 




0 ' 0 


ft ft 

I K xxt ( 0 . y . s-T 5 1 lu(y.T) - v(y,T)| |yl dy dx ds 

J 


A 


3 C (u - v! I* 1 2 f |K (0,y,s-x)j |yl dy dx ds 

00 J * J J « * Xt 


t "*0 "0 

z 


& 4C |u - Vi C ( 0 , 0 , 1 } It - t 

03 4 1 2 


1/2 


(5.4.23) 


Now we sum up Inequalities (5.4.22) and (5.4.23) to obtain the following bound 
for I 

11,1 


II ,1 a II , J + II , J 

11,1 11 , 1,1 11 , 1,2 


£ 4 


If - g_L + C iu - v| 

2 2 00 co 


C 4 (0,0,1) |t 4 - t 2 l 


1/2 


(5.4.24) 


Now 


taking f € Lip class with Lipschitz constant D, the estimate for 1^ 2 is 
given by 


II I £ 4 

11,2 


If - g | + D |u - vi 

1 1 00 


] C (0,0,1) It - t I 
00 J 4 12 


1/2 


(5.4.25) 


Thus from (5.4.24) and (5.4.25) , one obtains 
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I 1 ,,' 5 I".,.. 1 '"u.t 1 


s 4 |f - g I * C u - V 

2 2 co » 


C (0. 

4 

f 


0,11 it - t i 

i 


. If l - 


+ D iu - v| C (0,0,1) It - t J 

<c 4 1 « 


(5.4.26) 


Hence from (5.4.8), (5.4.15), (5.4.19) and (5.4.26) one obtains the estimate 
for I 7 as given below. 


Ii 7 l = fjtettj)) - 


f , toe t )) - g (0(t )) 

12 5 2 


li I + III + li J + li 

8 9 10 11 


If l " g i'« {|V xx ! « )0< ~^r 2 ’ <<X+U c 3 (0,0. 2 «) It, - t 2 r 


+ C lu - v | C (0,0,1) |t - t 1 

s x x CO 3 12 


„ C (0,0,2*) 

+ 2 |f - g | (lv | ) -JU, — It - t I 

1 1 <* xx co oc ( 1 - <x ) % 2 


* ic . iu , - \i„ c j°’ o -v i‘, - y 

* S 10,0,11 [ lf 2 - y„ + c lu - v >„ ] it, - s 

♦ C 3 (0,0,1) |f, - g,l m + D |u - v| 1 1 1 3 - t 2 l 


V°.°.D [ If, - 8,1. ♦ C lu-vl.] 


It, - t 2 l 


* *C t (0,0,l) [ If, - g,| m * D lu - vl^ J |t, - t, 


(5.4.27) 


Denoting C (0,0,«) = CI - «j c (O.O.ZoO ♦ 2““ C (0.0.2«> 


C b (0,0,l) = C s C 3 (0,0,1) + 4 (^(0,0,1) 


(5,4.28) 
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We obtain 


C £0,0,0.) 

iu s “ A n r-~irr if i ~ •*. - u* 


l ot xx co 1 2 


+ C (0,0, 1 ) |u - v | it - t 

b x x co l 2 


1/2 


+ C (0,0,13 

b 


if - g,L + C |u - v| 

2 2 co oo 


u, - t 2 i 


1/2 


+ c ( 0 , 0 , 1 ) 

b 


If, - g, l„ ♦ D |u - v|_ 


it, - y 


1/2 


Now we proceed to evaluate the oc-semi norm of }f ^ - g |. For 

(5.4.29) is divided by left } - 6(t ) I*. 

1 2 


Therefore, 



f (e(t )) - i (e(t )) 

Jr „ 


ii i i 



f(e(t 2 )) - gj ( e ( ) ) 


lettj) - e£t 2 ) | v 


C (0,0, «) 

ft 




« XX CO 


+ C (0,0,1) |u - V | |t - t | 

b x x » 1 2 


(1/2) -« 


+ C (0,0,1) 

D 


If - gj + C |u - v| 

2 2 oo 


l U 1 ' 

CO J 1 


t 2 l 


( 1 / 2 ) -« 


C b (0.0,l) [ If, 8,l tt + D|u - v l(o] ' V 


( 1 / 2 ) -oc 


It * 


J le(t i ) 


(5.4.29) 

that purpose 

]| 


e(t 2 )l“ 


Hence 



134 


( iv 1 ) 

If, - g.L * 

1 1 “ Unfle'i) 


C (O.o.c) 

® { 5 - a ) ’ "! # 


I*. ■ B. 


1 


(Inf f0' i 3 


|u - v I * 
x * a 


If " P,l 4 £' ~ v| 

4 * €D < 


l f . “ g,l 4 D |u - v| 
1 1 so 


C (0.0.1) it - t 

b 12 




(5.4.30) 

Equation (5.4.3) expresses the exact expression for f g - g^, One can see that 

an estimate similar to the expression (5.4.30 3 will be obtained for if - Z | 

2 6 2 « 

as well. Here the only difference Is that the kernel is evaluated at X = 1. 
Therefore, the estimate has the form 


If 2 ’ g 2 '« ^ 


Iv | 1 

ex 

C <1,0, «) 

XX CO 


A 

infle' | J 


« { i - <x)’ 


|u - V | + 

x x ce 


(inf |0' I ) 

+ [ |f i " g iL + D i u - v l 
Equation (5.4.1) implies the following estimate : 


If " 8J 4 C |u - v! 

2 2 co co 


C (1,0,1) it - t (5.4.31) 

D 12 


IT e., (f i • V - T e. q (g i • g 2 >l„ = ( If, - i,l„ * lf„ - i 2 l„ ] 


I Iv II* 

XX CO 


[ infle 


i 1 

— I C (0,«) If - g I ] 

' | L c i *v«J 
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1 


(Inf IS' |) 


‘ V, - { l f 2 - 8,1. + C lu - vl^ 


{ 


l “i to ~ ” } C c * 0 - 1 ' I tj - t 2 l (1/2) “ (5.4.32) 


* < If. - 8.L ♦ D lu - v | 


Where 


C (0,«) 


c ( 0 , 0 , «) c ( 1 . 0 , oc) 

fi ft 

« n ~ «) « (i - oc) 


C e (0,l) = C b (0,0,l) + C b (l,0,l) (5.4.33) 

The estimates for |u^ - vj^ and ju - vl^ can be given by (5.3.36) and (5.3. 
35) respectively, l.e.. 


lu - v 1 

X X » 


f \ { \ 

A A t + A Jf-g| + |f - g I 

1 7 1 2 1 *1 00 2 *2 00 J 


1/2 


(5.4.34) 


|u - vl^ a A 7 [ If, - 8,1. * If, - 8,1. 


(5.4.35) 


Now we are in a position to write the complete expression for the <x-norm of 

V ,- 

■ y - T e., ( *. • yi« - 'v,(f,. y - • 8 ,)i. 

* IT e,, (f i’ y • 8 2 )l « (5 - 4 - 36) 

Substituting (5.3.35) and (5.3.36) in (5.3.20) and abbreviating the constants, 
the 1st term of R. H. S. of the above expression is repeated for convenience, 

IT e,, ( y y - T e.,y • 


* A e |f i ■ *.'« 1 


<x/2 
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♦ \ A , A 7 S * A * i "l ‘ "l". * lf .’ * J 1 


♦ A |f - g I * If, - R,i 1 
10 5 °> » 2 * * J 


+ A A, If, * g I ♦ If ~ K J K «' 

11 7 1 1 05 « * * 


(5.4.37) 


Similarly, substituting (5.4.34) and (5.4.35) in (5.4, X> i the second term of 
R. H. S. of (5.4.36) becomes 


1T e, q (f r V " T e,q {g i ’ h )l <* * h iz if $ ~ 8 i ! c 


+ A A A t + A 

13 17 1 2 


4A .3 ( '5 lf 2 -* 2 '.) 15 - 51 " ,2 "" 

4A i« 5 ( '5 ' 5 1 .* '5 ' 55 ) l ! /2 I5"5 


If - g | + if - g I t l/2 j t - t j £5/2)_c * 

1 °1 CO 2 r 2 oo 1 ' 1 2> 


( 1/2) -» 


(5.4.38) 


Considering individual norm to be less than the «-norm, (5.4.37) and (5.4.38) 
together will give us 


- T 8,q <g l 

. j 

A t W2 * 

A 

8 

12 

f* (a 

A 1 

l » l ' 

7 

> 

> 

t 

13 ^ 1 7 

1 

f A.. t 1 ' 2 

+ J 


, " 8,1* + l f 2 ~ « 2 1 W 


t + A t 
2 


(l/2}-oc 


* 15 - tj"^) 


(5.4.39) 


(5.4.40) 


and therefore, the map T 0 is a contraction if A <1. 

IS 
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Let us repeat A ( , 1 = 1 15 for convenience. 


f i /2 'I 

A 1 = c 2 ( x > 0 , 1 ) t + C^x.O.Q) I x j I (C + D) + D tyx.O.O) 

’. 0 , 0 ) J 


A = 


C 2 (x,0,l) t"‘ + C 2 (x,0,0) |x| + C 2 (x, 


C (x,0, 1 ) t 1/Z + C (x,0,0) x 
2 2 


A = 


A = 


| C^x.0,1) t 1/2 + C^x.0.0) j x | j (C + D) + D CJx.0,0) 

^(x.0,0) J 

= [ C^x.O,: 


C (x.O.l) t I/2 + C^x.0,0) | x | + C 


1/2 


1) t + C (x,0,0) | x 


1 - A t 

4 


A = f C (0,0,2«) + C (1,0,2a) ] (|v | ) a 2~ (<x+1) _L 

8 1 3 3 J xx oo <x 


A = C C (0,0,1) + C (1 

9 s 3 3 


,0,1) | 


A 10 = | C 3 (0,0,1) + C 3 (1,0,1) j 


( C (0,0,1) + C (1,0,1) ] 

A 

= C (0,«) 

|v | ) 

1 XX 1 CO 

(3 3 J 

12 

c 

. l«'l . 


A 3 - C c (0,«) 


(Inf |e' I )* 


A - (C + D) C (0,1) ■ 

14 c (inf |0' | ) 


that 


In order to obtain T to be a center contraction map we need to show 
' |v ' 

is bounded. This can be achieved by utilizing the 


XX ‘ oo 


l !*'l J 

representation (5.3*37) and the steps (4.5.39)~{4.5.43) of chapter IV. Hence 
using the flat initial value, we can see that T Q is indeed a center 
contraction map. 



Now we define the sequence of iterates 4 p ( [ * r, r i * 1 . t 


defined by 


> (n * 1J {eu)i - t„ Ip ,ni 1 (t) 

i Cl, q i 


since f is a (e.q)-fixed point of T , clearly ty th«- estimate 

i w,<i 


if - p <n+55 I 5 A jf - p M \ 

» i *i *« IS 1 ) ' l •< 


and consequently p 


f ( - — > 0 as n ■ -» ®, Since C is complete, 


f in C . To see that T has a unique fixed point , we proceed by the method 
i y*<j 

of contradiction. Suppose h * (h ,h ) be another fixed point of T , 

t 2 6,q 

different from f but we know that 


» - f l„ - I WV - T e I. * V ! h " f L < l» - f|. 


This is a contradiction. This shows that f * h. 


This competes the basic assertions of this chapter . | 


5.5 CONCLUSIONS : 

In this chapter we have shown the uniqueness of a solution which is 
established by the method of center contraction. With the assumptions of the 
flat initial values the result is shown to be true for small time interval. 



APPENDIX-A 


COMPUTATION OF DOMAIN CONSTANTS 


The terms of the following form are needed to be estimated. 


- V J 


t r 00 I a J 


t J o 1 dx 
2 


- KCx.y.^-x) 


I x— y | dy dr , 


for 


j = 0, 1, 2, ; 0<«<1. 


In the case j = 0, (x,y are positive) 


t 00 


I Q = J 1 J I K(x,y, t j -x ) | Ix-yT dy dr 


t o 
2 




j s f [ 


t "0 
2 


-<x-y) 2 /4(t -T) -(x+y) 2 /4(t -T) 

* i- + i— ] ix- y r 


Vt~~T 

1 


Vt-r 

l 


± V 1 


t. ^ -(x-y) /4(t -T) 

e — | x-y I dy dx 


ViT J t J 0 Vt -T 
2 1 


let <r = 


y - x 


2Vt - x 

i 


for t > x, then 


t JCD 


I S 
0 




VST 


— 0- 0< _OC , . , we. . _ 

e (T 2 |t -t) acr dx 

X 1 


, oc/2 


Vt ~X 
1 


dy dx 
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i m 0 

2“ f 1 ,, i® /2 r 

s j It -t| r e 

* V M 


4<r rft 




m 


v'tT 


HX €* _ 

e $ 


ft - t 
1 # 


{ €* ' # > * 1 


]Mr» V 


- CM Cy t 2 ) 


«*/2>»S 


(Al) 


2 «« 

r 

r ® 2 

-4T r* 

where C {«) * — 


e cr dcr 

1 C2+«)v / x~ 

■ m 4 



In the case x = 0, 


,<x 


I i 


ViT 


[ 


,00 2 

ht a 
e <r dcr 


{t - tj 
a 2 


t€»/2Ul 


l * (<*/2 ) 


= (1/2) C (oc) (t - t ) 

1 12 


(*x/ 2 ) + l 


(A2) 


When J = 1, 


K (x,y,t) = 


-A 


2VT t l 2VT“ 


(x-y) , 2 (x+y) 2 \ 

„-*ix~y> /4t , J (x-*y > /4t 

— — — e + e 


2vT" 


Since |x - y| < |x + y| 


I 

1 


-r*r 

IMx.y.t^ 

-x)| | x-y l" dy dx 


* r* r 

1 

-<x-y) 2 /4(t -X) 

“ 1 

lx-y| 

*2 ° 

2t/h~ 

1 t i -x 1 

2Vl x -x 

*r* r 

1 e' 

•(x+y) /itt^-X) 

lx+yl 

V° 

2Vn~ 

It -x| 

2VT~^ 


I x-y I “ dy dx 


lx+y|“ dy dx 
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Therefore, 


i s i rS r“ e - (x - y) 2 /<u r T) !x ~y |14tx 

1 * 2 VT K ^ -CO 


I tj-T I 


2vt^T 

1 


dy dx 


y-x 


with the substitution cr = 


2vt-x 

i 


for t > x, 


oc _t ^oo 2 


i a 


1 2 VT 


J.’J 


t -00 
2 


-cr l+oc . {cx-D/2 , , 

e c r 1 -t | dcr dx 


2 

vr 


I 


-cr l+oc 

e cr dcr 


(t r V 


(OC+D/2 


1 + OC 


C (<*) (t - t ) 
2 12 


«x+l)/2 


oc co 2 

where C (oc) = — j cr 1 *** dcr 

2 VtT (1+oc) -oo 


if x = 0, 


I £ __L_ C (oc) (t - t ) (0C+I>/2 
1 _ 2 12 


When j = 2, 


K (x,y, t) 

XX 


zViT t { 


-<x-y) / 4 t 


(x-y)‘ 


2t 4 t 


(x+y) 


2t 4 t 


2 T 


-(x+y) /4t 


(A3) 


(A4) 



I *k£* 


I « f 1 f IK (x.y, t -t ) I Ix-y I" rtv dt 

2 j J xx * 
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-(*-y» 2 /«U -t! 

e i 


Ix -y s' ’ 

• 

1 
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\ 
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"*0 

2 

It -T| J " 

\ 

% 

* 

4 1 1 - T I 

j J 

i 

rr 

S ° 

2 

-(x»y) /4lt -T) 
e I 

1 

(x* y r' ’ 

# 

2^r 

. . , 3/2 

t -T 

1 

2 

411 -t | 

1 J 


y I dy dt 


lx yl dy dt 


Since I x - y J a 1 x + y| 


— 1 1 f 

y^T j* j 


t -m -(x-y) /ill -TS 
e i 


2v^r It r xi 3/2 


1 

2 


Ix-y I 


4 1 i } -t 


|x - y| dy dT 


2“ f° -cr 2 f l a 

J e — + * 

— ** 


v'r 


■r” » f ' It, - T I IT 
l 2 


c (oc) 


(t - 1 r 

1 2 


(A5) 


oc to 2 

Where C («) = _ f g -0 " I 1 

3 Vn~ J-co 


♦ <r 


« . 

0 * d<r 


if x = 0, then 


I 5 
2 


c 3 («) 

2<x 


(t - t J* 72 
1 2 


(A6) 
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The estimate for 


K (x,y, t ) * 

XXt 


1 

f 2 

-(x-y) / 41 

|tv * 

6 

(x-y) 2 

■f 

(x-y) 4 ( 3 j 

?./n 

\ 

8 

t ( 7 / 2 ) 

16 t <9/2 > ’ 3 t s/2 


2 j 

- {x+y ) / 4 t 
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+ - + 0 


¥ C 


8 t (7/2) 
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-J 
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1 


t ** a 
2 
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xxt 1 


<x -oo 2 r 

2 f e HX 


Vn~ J o 


** r 


2cr 4 + 6cr 2 + (3/2) 


1 : >v 


ti ( “ / 2) - 2 dr 


(A7) 


(t - t 

s C («) - 

4 1 - («/ 2 ) 


(A8) 


where C («) = 

4 VT J o 


I, 


.03 2 

-cr oc 

e cr 


2<r 4 + 6<r 2 + (3/2) dcr 


(For complete calculation please see Appendix B (B8)) 


Let the value of I is estimated, where 

4 


,t ,00 
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^ J * J r\ 


t o 
2 


C (x.y.t) = — f 

1 t 
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Now 
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1 


lx 
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dy 
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4 2 VT 
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1 2 1 
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J e 
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cr dcr 
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-« 


J 
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C (oc) 


(t - 1 r 


oc 12 


c 5 <«) = 


vr 


2 r 



f HT J 
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J e 

"2- + c 

O' d<r 

" —CO t 



V 


j 


yl 0 * dy dx 


yl** dy dx 


yl* dy dx 


dx 


(A9) 


Where 



APPENDIX-B 


COMPUTATION OF DOMAIN CONSTANTS 


We will need to estimate terms of the form 


t -! 


I = 
J 


J. J 


, dx 
O ■'0 1 J 


K(x,y,t-r) 


lx - 2n - yl dy dr ; 


for j = 0, 1, 2, and 0 < « s 1 


r 1 

Where K(x,y,t) « ^ 


nr - 


« -/h n t 


2 2 -i 

• £x-y-2n) /4t . -(x+y-2n) /it 

+ e 


In case J * 0 (x, y are positive) 


t 


I = f f | K(x,y, t-r) ] jx - 2n - y|“ dy dr 
® "0 "^0 


= rf I — r 

J 0 J 0 — *rT L 


-(x-y*-2n) / 4 (t-T) 


n=~oo /4ir(t - tT 


+ e 


(x+y-2n) /4(t-T) 


J | |x - y - 2n| 


I* J' 2 [ ^ [ e' (x 

J 0 J 0 n=0 ^471 1 1 - tT *• 


:-y-2n) /4(t-T) 


+ e "(x+y- 2 n) 2 / 4 (t-T) | x - y - 2n| 


dy dT 


dy dx 
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since |x-y-2n| * jx + y - 2nj , n*=0 , 


then 


-{x-y-2n) 2 /4«t-T) . - <x*y-2n> 2 /4 ( t -T 1 

e * e 


so 




1 

Vir j t — x| 


2 

- <x-y-£n) /4U-T) , 

e x 


y 


?nj dy dx 


y - (x - 2n) 

Let us consider, cr = — 

2v^t - x 


, t > x then 


dcr = 



dy 


Hence 


1 -(x-2n) 


= 0 


4 « Vt-T 


'.*/ 7= I I 

J 0 V7T n = 0 J (x~ 


2 

~CT « „ A 

e <r dcr x 2 j t 


x dcr dx 


2n) 


2V^ 


0 


a+2 , ® ft 2 


7n~ 


f [ r 2 e- s * i f | t -xr /2 dx 

^ n=0 J e > J o 


oc+2 co ft 2 


I J 

=n J / 


VrT n=0 J 0 


2 hT <x 

e cr dcr x 


,«/2+l 


(«/2) + 1 


= C^Cx.n.oc) t 


OC/2+1 


(Bl) 


Where C^x.n.oc) 


.oc+3 


(a + 2} VtT" *- n = 0 ''O 


co z%0 2 

2 -O' oc , 

e cr dcr 


’ i r 

- « = A J C 


(B2) 


When j = 1 ; 


I 

l 


= «f 0 lo l K x (x,y,t " T) l l x - 2n - y|“ dy dx 



where 
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K (x,y, t ) * Y 

V / ■* fm** 

n--co 1/4 Tt t 


(x - y - 2n) 2 

-(x-y-2n) /4t 


x e 


2 t 


(x + y - 2n) 


2 t 


x e 


- (x+y-2n) /4t 
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= j 
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A 
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0 J 
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* 

1. • 

0 
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2 (t - t) 
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(x + y - 2n) 2 
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2 

00 

1 
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x e 
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- 1 x-2n-y I /4 I t-T | - | x-2n+y | /4 1 t-T | 

6 & 0 


Therefore, one can write 1^ as 
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r l r 1 ? 2 

r s ^ ( n 

x - 2n ~ y , ^ 2 , . 

5 7 {x-y- 2h) /4 | t ~ T 1 

I 2 V 

— . — - * — f* 

•'o J o n=0 /4rt!t - tT 

2 It - T| 


dy dx 


y - ( x - 2n 5 
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2 /r^"x 
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, dcr dx 
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o“ +3 » r G 
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Vn - (ot + 1 ) n=0 J 8 
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Now for J = 2, one gets 


r i: 


2 J o J o I a x 2 


K(x,y,t~x ) |x - 2n - y| a dy dr 


Integrating the Neumann function twice with respect to x, it becomes 
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1 00 
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<x 
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Finally, we consider I where 
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, CO r 2 
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Now we estimate where 
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